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Abstract 


In this second technical report for the investigation of the effect 
of atmospheric dust on the estimation of total ozone from the earth’s 
ultraviolet reflectivity along a local nadir direction, we will describe 
results of a detailed analysis of the simulated measurements for the BUV 
(Nimbus-4) configuration by using a total ozone estimation procedure 
similar to the one reported by Matoor, Heath, Kruegar (1971). Based 
on this detailed analysis, we will discuss the following two fcitems: 

(I) A set of recommedations for increasing the accuracy and confidence 
level of the total ozone values estimated from the measurements of the 
earth's ultraviolet reflectivity at five different wavelengths (BUV con- 
figuration); and (2) a tentative procedure for the estimation of total 
ozone from measurements of reflectivity at six different wavelengths 
specified in. the SBUV/TOMS (Nimbus'tG) configuration. - 

The above-menti noed analysis .is based on simulated measurements for 
pseudo-spherical models of the earth's atmosphere resting on a Lambert 
surface and having arbitrary vertical profiles of ozone, but otherwise 
free of aerosols and/or water droplets. 



I. INTRODUCTION 

In the first technical report with the same title as this one, we 
described a set of computer programs (FORTRAN IV language) for evaluating 
the azimuth- independent component of the intensity of scattered radiation 
emerging at the top of a pseudo-spherical atmospheric model with arbi- 
trary height-distributions of ozone, and/or aerosols. The model is 
assumed to rest on a Lambert ground of reflectivity R. By pseudo- 
spherical , we mean that the basic atmospheric model corresponds to that 
of a plane-parallel atmosphere of infinite extent along the horizontal 
directions, but of finite optical extent along the vertical; the 
sphericity is only partly accounted for by computing attenuation suffered 
by the incoming solar radiation (neglecting refraction) for the appropriate 
spherical case. Furthermore, the nadir angle of observation is also duly 
corrected for the satellite’s altitude above the earth's surface. 

These computer programs were primarily developed for assisting the 
NASA, Ultraviolet Ozone Team (Dr. D. F. Heath and others) in their total- 
ozone estimation work by performing analysis of the simulated measure- 
ments for various realistic models of the earth’ s atmosphere. Estimation 
of total ozone in an atmospheric column from multi-spectral measurements 
of the radiation backscattered by the earth-atmosphere system, is not a 
trivial problem for several reasons. First, a relatively noise-free set 
of spectrally pure measurements is required with an accurate knowledge; 
of the instrument's position and attitude in space and time. Second, 
there are a good number of atmospheric unknowns which can modulate the 
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signal in an undeterminable manner. Besides the total ozone to be 
estimated, some of the other unknowns are its height distribution, 
reflectivity of the underlying surface and its wavelength dependence, 
surface pressure, aerosols, water droplets, ice crystals, and nonhomo- 
geneities in the field of view due to one or more of these factors. 

Some aspects of this complex problem were first investigated by 
Dave and Mateer [1967] with simulated measurements restricted to the 
Rayleigh models of the earth-atmosphere system. Based on this prior 
study, Mateer, Heath, and Krueger [1971] developed a procedure for esti- 
mating total ozone, in an atmospheric column underneath the NIMBUS-IV 
satellite, from measurements of the back-scattered ultraviolet radiation 
in five spectral regions (viz., 0.3125, 0.3175, 0.3312, 0.3398, and 
0.3800 pm). They then used this procedure for estimating total ozone 
for a sample of 320 cases of approximate coincidences (in space as well 
as in time) of the NIMBUS-IV data, . and ground-truth measurements of total 
ozone with the Dobson Spectrophotometer. They find that, on the average, 
the satellite values were lower than the corresponding Dobson values by 
about 0.025 atm-cm. Some cases are also reported where the difference 
between these two values is as large as 0.060 atm-cm. They believed 
that a part of this error is due to a lack of a perfect coincidence be- 
tween the Dobson, and the satellite data. Based on their experience with 
the analysis of the NIMBUS-IV BUV data, this team of investigators have 
recommended an additional measurement at 0.3600 pm for the SBUV/T0MS ex- 
periments to be performed aboard the forthcoming NIMBUS-G satellite. 
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The purpose of this report is to discuss the results of a detailed 
analysis of the simulated measurements for the five-wavelength, BUV con- 
figuration. These simulated measurements are obtained by making use of 
the computer programs mentioned at the beginning of this section. This 
analysis is for models of the earth-atmosphere system resting on a Lambert 
surface whose reflectivity R may or may not vary with the wavelength, but 
the atmosphere free of any particulate matter comparable in size to, or 
larger than, the wavelengths of incident radiations under investigation. 

A set of five-wavelength simulated measurements is first analyzed in 
terms of the best ozone estimate , and the effective albedo by making use 
of a procedure identical to the one discussed by Mateer, Heath, and 
Krueger [197.1]. Further analysis is then performed by imposing addi- 
tional, valid and meaningful constraints for acceptance of the best 
ozone estimate. This analysis is followed by a set of recommendations 
which may be incorporated in the final, BUV, total-ozone estimation pro- 
cedure for increasing the accuracy and confidence level of the satellite, 
total-ozone values. Experience gained from the aforementioned BUV simu- 
lation study, is then applied for developing and testing a total-ozone 
estimation procedure for the six-wavelength, SBUV/TOMS configuration. 


II. TOTAL-OZONE ESTIMATION PROCEDURE 

2.1 General: In §1, it was mentioned that a procedure for estimation 

of total ozone from five spectral measurements of the radiation back- 
scattered by the earth-atmosphere system along the local nadir direction. 
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was developed and used by Mateer, Heath, and Krueger [1971]. Since a 
detailed written description of this procedure was not available at the 
beginning of this project, a tentative description was developed by the 
author after several telephone conversations with Mr. Larry Novak who is 
very closely associated with the BUV, total-ozone-analysis project at 
NASA/GSFC. This tentative description was sent to Dr. Robert S. Fraser 
(Consultant to the Technical Officer for this Contract) for a verifica- 
tion on August 12, 1976. An additional point to this tentative descrip- 
tion was brought out by Mr. Novak in the material which Dr. Fraser en- 
closed with his letter to the author, and dated August 25, 1976. In the 
meantime, Mr. Novak also sent, to the author, a set of subroutines used 
in estimation of total *ozone, . and a copy of his preliminary dpcumentation 
of these subroutines. 

Some basic information required in following the total-ozone esti- 
mation procedure described in §2.4, is given in §2.2 and §2.3. Descrip- 
tion of the procedure given in §2.4 is primarily for the purpose of 
bringing out details of the procedure as used in the interpretation of 
the simulated, BUV data. However, this description is kept general 
enough for use in modification to the original total-ozone estimation 
procedure (one being used at NASA/GSFC by Mr. Novak and others), and 
also for use in analysis of the simulated observations for the SBUV/TOMS 
configuration. 

Because of a rather limited scope of this simulation experiment and 
contro lability of several factors, some minor differences occur between 
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the original total-ozone estimation procedure (used at NASA/GSFC) , and 
the one used for the first analysis of the simulated, BUV observations. 
These differences are brought out in §2.5 for the sake of completeness. 

2.2 Basic Tables: The backbone of the total- ozone estimation procedure 

is several sets of tables of computed quantities I c (X, Poj R^ = 0.0, 

6 0 ), T^CX, ft, P 0 , 6 o) , and S c (X, ft, Po) obtained with the help of computer 
programs described in the Technical Report I having the same title as 
this one. Each table in these sets is for a pseudo-spherical atmospheric 
model with a surface pressure of Pq mb, and a total-ozone amount of ft 
atm-cm. These models resting on a perfectly absorbing surface (R^ = 0.0), 
are free of any .aerosols, water droplets, ice crystals, or horizontal 
nonhomogeneities. The remaining two parameters, viz X and 6q represent 
the wavelength (in pm, 1 p = IQ -4 cm) , and the local solar zenith angle 
(in degrees) for the incident radiation under investigation. For all 
wavelengths, these quantities are computed after assuming that a ttcosSq 
units of solar energy is incident normal to a unit horizontal area loca- 
ted at the top of the model, in each spectral band. 

The quantity I C (X, ft, Pq, R^ = 0.0, 0' o ) is the intensity of the 
scattered radiation emerging at the top of the atmospheric model along 
the local nadir direction when all four parameters in the parentheses 
are specified. The quantity T c (X, ft, P 0 , 0o) represents the total 
(direct + diffuse) transmission along the local nadir direction (in the 
same model) when the atmospheric model is illuminated isotropically, 
from below, by a Lambertian reflection of all direct solar energy arid 
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the diffuse skylight energy (due to illumination of the model from top 
only) emerging at the bottom. The third quantity £> c (X, £2, P 0 ) repre- 
sents the diffuse flux reflectivity for an isotropic illumination of the 
model from below. 


The following expression is then used for computations of 1^(1, £2, 
Pq, R-, , 9 q ) when a given model rests on a ground obeying Lambert’s law of 
reflection, and having a reflectivity R^ at the wavelength X: 

I C (X, £2, P 0 , R a , e 0 ) = I c (X, £2, P 0 , \ = 0.0, 0 O ) 


1 - R x S c (X, £2, P 0 ) c 


T (X, £2, P n , 9n) * (1) 


In all, 20 sets of tables are prepared; the first ten sets are for 
Pq = 1,000 mb, and the remaining ones are for P 0 = 400 mb. The total 
ozone amount for an atmospheric model with Pq — 1,000 mb is given by 
one of the following values of the parameter £2: £2 = 0.200, 0.250, 

, l 

0.300, 0.350, 0.400, 0.450, 0.500, 0.550, 0.600, and 0.650 atm-cm. (The 
vertical distribution of ozone for each of these amounts is for a mid- 
latitude atmosphere. The actual distribution used is based on informa- 
tion provided by Dr. C. L. Mateer.) 


The atmospheric models for Pq = 400 mb are obtained after termin- 
ating the corresponding 1,000 mb models at a level which is 7 km above 
the mean sea- level. Thus, if the 1,000 mb has a total ozone amount of 
0.200 atm-cm, the corresponding 400 mb model has an actual total ozone 
amount of 0.189 atm-cm. For clarity, we will denote the actual total 


7 



ozone amount in a model by the symbol ft, and the total ozone amount in 

the model when its lower boundary is extended to the mean sea-level, by 

the symbol ft. . Thus, for the above-mentioned case, ft = ft. = 0.200 
m in 

atm-cm for the Pn = 1,000 mb model: but ft = 0.189 atm-cm, and ft. = 
u m 

0.200 atm-cm for the Pn = 400 mb model. Values of ft and ft. for 

u in 

the remaining models are given in §3.1. For convenience, all ten sets 
of tables for a given value of Pq are arranged in the increasing order 
of the parameter ft. 

In order to confirm with the original, BTJV total-ozone estimation 
procedure, we have used the value of ft in the place of that of ft 
during the analysis of simulated observations. This ingenious step per- 
mits an estimation of total ozone in the entire (right down to the mean 
sea- level) atmospheric column underneath the satellite even when an op- 
tically thick cloud layer is located between the satellite and the 
. ground. Presence of such a cloud layer is not expected to affect the 
amount of ozone within, or under, the cloud layer. 

A given set of tables contains values of I , T , and S at six 

c c c 

wavelengths (viz, 0.3125, 0.3175, 0.3312, 0.3398, 0.3600, and 0.3800 pm), 
and those of I and at ten values of 9g (viz, 0°, 45°, -60°, 70°, 
75.6°, 79.6°, 82.5°, 84.7°, 86.7°, and 90°). The 0.3600 ym table is for 
use in the analysis of the simulated measurements for the SBUV/T0MS con- 
figuration. 

If the conditions for which these basic tables are prepared are 
• satisfied within acceptable limits, and if values of the parameters Pq, 
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0q» and are available when the intensity measurements [ 1^ ( A , 6 q)] 
along a local nadir direction are taken, a value of ft for the atmo- 
spheric column underneath the satellite at the time of observation may 
be estimated by making use of secondary tables of I (A, ft, Pq, R , 0q) 
vs. Pq, vs. ft (or ft^ ) , and vs. 0 q. Some interpolation between tabu- 
lated values of various variables will be necessary for this purpose. 

Out of the several values of ft so obtained (one corresponding to each 
wavelength of measurement), one with the highest sensitivity of 1^ to 
changes in ft, in the range of Immediate interest, may be selected as 
the best ozone estimate. [Only two parameters are listed for the mea- 
sured intensity, I (A, 0 q). This is because the third parameter ft is 
to be determined. Furthermore, as it will be discussed later, values of 
the remaining two parameters (viz, Pq and R ) are 'also unknown within 
any meaningful accuracy for either the BUV, or SBUV/TOMS configuration. 3 

Before closing this subsection, we would like to add that the above- 
mentioned basic tables can also be used to obtain a value of after 

rewriting Eq. (1) as follows: 

f (A, ft, P 0 , R = 0.0, 9 0 ) 

R „ 1 (2) 

^ T c (a, ft>?0> 6o) + S c (A, ft, Po) f(A, ft, P 0 , R^ = 0.G, 0 q) 

where 

f (A, ft,P 0 , R x = 0.0, 6 0 ) = I m (A, 6 0 ) - I c (A,ft, P 0 ,R x = O.O, 9 0 ) . 

(3) 

The ft dependence of the terms on the righthand side of Eq. (2) vanishes 
when the wavelength for which the measurement is available, is located 
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outside the bands of absorption by ozone (e.g., 0.3800 ym) . 


2.3 The Quantity N: It is a well-known fact that the atmospheric con- 

ditions for which the basic tables (§ 2 . 2 ) are prepared, are rarely satis- 
fied in practice. Furthermore, values of 1^(1, may have some biased 
error due to calibration problems and/or deterioration. It is therefore 
customary to use the ratio of the intensity at two wavelengths with 
significantly different amounts of absorption by ozone, but located 
sufficiently close enough to minimize the effects of other uncertainties. 
A more convenient quantity used in the ozone estimation field is given 
by the following equation: 

N c (Xi, X j ; Rx j » ^in 5 ®o) = 100 

X c (^i> ^in s p 0> ®-Xj[> 

* 10810 M x j. “in' p 0> 0o)J ’ 

(4) 


|VV 9 0>] 

W 0 O ). - 100 X 1 O g l0 • 


For a given set of the remaining parameters, the quantity N c generally 

increases with increase in (see §IV). Therefore, a value of Q 

in 

can be estimated from a measured value of N (X., X., 0n) and appropriate 

m r 3 

sets of tables of vs. Pq, vs. 0 q, and vs. (some interpolation 

would be required) provided information about the parameters R^, ^X j 5 
Pq, and 0 q is available at the time and the place of observation. For 
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the BUV (or SBUV/TOMS) configuration, the criterion for selecting the 
wavelength pair mentioned at the beginning of this subsection is satis- 
fied by two pairs, viz, 0.3312 - 0.3125 pm, and 0.3398-0.3175 pm. The 
values for the third pair, viz, 0.3398 - 0.3312 pm, are generally too 
insensitive to changes in ft^. 

Another quantity of interest to us is the slope of the vs. ft _.^ 
curve in the region surrounding the point where an ft value is esti- 
mated. It will be given by 


lZ - vv 

' x ,y - Vl - °k 


where 


ft . £ ft ^ ft. . 
k x,y k+1 


In Eq. (6), we have omitted the parameters A-;, A-j , R^., Ra-> ^0 an< ^ ®0 

j x j • 

for the quantity N in the numerator, for brevity. The quantities ft. 

c 

and ft, ,, are two successive values of ft. for which basic tables are 
k+1 m 

available. 


The subscript x (to ft or S) stands for the wavelength pair used; 
x = 1 for the first pair, viz, 0.3312-0.3125 pm; and x = 2 for the 
second pair, viz, 0.3398- 0.3175 pm. The subscript y stands for the 
value of P 0 identifying the set of basic tables. The superscript Z 
identifies the nature of Rj^, ^-Aj information (see §2.4). 
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2; 4 Basic "Procedure: For estimation of total ozone in a simulated BUV 

configuration, simultaneous measurements of the intensity of back- 
scattered radiation along the local nadir direction are available at 
five wavelengths; i.e., we have values of 1^(1, 0o) f° r ^ = 0.3125, 
0.3175, 0.3312, 0.3398, and 0.3800 ym. Value of Qq for the place, and 
the time of observation is also known. We then proceed as follows to 
obtain values of the best ozone estimate (Q^) , and the effective albedo 
<£). 

Step I: Use the measured value of I m (0.3800, 9 0 ) for evaluation 

of the righthand side of Eq. (2) by making use of the 1,000 mb tables 
for the appropriate 0 q. (We do not have information on Pq 5 the 1,000 
mb table is used arbitrarily.) Furthermore, there may be other unknown 
factors, e.g., absorption and/or back-scattering by aerosols, or the 
Lambert’s law may not be' valid. We have therefore not computed Rq ^goo 
in this step. The computed quantity will therefore be referred to as 
the course effective albedo 3 and it will be represented by the symbol 

R C 

0.3800,1000’ 

Step II: Assume that the coarse effective albedo so measured is 

independent of X. Therefore, 

R 0 _ R 0 ,-jK 

1,1000 0.3800,1000 v J 

for X = 0.3125, 0.3175, 0.3312, and 0.3398 ym. 

C 

Step III: Use values of the quantity R^ ^qqq ^- n t ^ ie pl ace of re- 
flectivities Rj^ and R^ in Eq. (4). Then obtain values of the quan- 
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- Q - - Q — 

titles £2 , nnn and S , for x = 1, and 2 by making use of the 

x,1000 XjlOOO 

P 0 = 1,000 mb tables, information content of §2.3, and values of N m (A_^, 
X , 9 0 ) obtained from appropriate values of 1^(1,' 0 O ). 

C C 

Note: (1) The quantity £2 , AnA is set to 0.0, and S 1AAA is set to 

n x,1000 x,1000 

— 100.0 when the quantity is found to be greater than 

all values available for the given combination of vari- 

ous parameters. 

(2) As mentioned earlier, the quantity generally increases 

with £2. . Thus, if the N value under scrutiny is less 

m m 

C 

than H for £2. = 0.200 atm-cm, values of £2 A „ A and 

c xn. x,1000 

Q 

S x ^000 are °^ ta ^ ne ^ carrying out a backward extrapola- 
tion -with values of N at £2. = 0.200 and 0.250 atm-cm. 

c • in 

c 

Step IV: Select one of the two £2 , AAA values for a use in Step 

r • x,1000 

V with the help of the following criterion: Set x to 2. Reset x to 1 

if, and only if, s ljl0 00 > S 2,1000' 

c c 

If £2 1 1000 = £2 2 1000 = 0.0 (i.e. , it was not possible to estimate 
total ozone for either of the pairs in Step III) , initialize the follow- 
ing quantities as shown for clarification, and omit Steps V and VI. 

^x, 1000 = °*°’ S x,1000 = ~ 100 - 0 ’ for x = 1, 2; and R Aji000 = R X,1000 
for the first four wavelengths. 


Step V: Use the value of £2^ j qqq selected in Step IV to compute 

a value of the improved effective albedo 3 Rq 33 gg as follows: 

Determine two successive values of £2. for which the basic tables are 

m 
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available so that £2 < £2 C . < £2, - . Obtain two values of 

k XjlOOO k+1 

I 

Rq 339 g 1000 ( one with £2^ and the other with £2^ + ^) , by making use of 
I m (0.3398, 0 O ), 1,000 mb tables for the appropriate £2^ and 9 q, and 
Eq. (2), These two values of R^ 3398 1000 are t * ien usec ^ to obtain a 

c 

final value of this quantity at £2^ ^qqq after a linear interpolation 

between £2, and £2, , , . 

k k+1 

Then, set R Xjl000 = R o. 3398, 1000 f ° r X = 0>3125 ’ 0-3175, and 
0.3312 pm. (Improved effective albedo is also assumed to be independent 
of wavelength.) 

Step VI: Repeat Step III after using R^ ^qqq in the place of 

C 1 XX 

R^ j,000" w iii then have values of the quantities £2^ ^qqo anc ^ 1000 

for x = 1, and 2. - . - 

Step VII: Repeat Steps I through VI above by using the 400 mb tab- 


les instead of the 1,000 mb tables. We will then have values of R 

,1 


for all five values of X, those of R 


X,400 


X , 400 

for the first four values 
I 


C C I 

of X, and of £2 S , £2 as well as S'"" , for x = 1 

x,400 x,400’ x,400 x,400 


and 2. 


Step VIII: Determine the effective albedo (R) as follows: 


R = 0.5 


R 1 +R 1 

[_ X,1000 x,400j 


( 8 ) 


using any value of the parameter X as is assumed to be indepen- 

dent of wavelength. 
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Step IX: This final step is for obtaining a value for the quantity 

' X 

ft , the best ozone estimate , from the computed values of ft for 
e’ x,y 

x = 1, and 2, and for y = 1,000 and 400 mb. The following paths are 
taken depending upon the availability of two, or all of values 

x ,y 

I I 

(a zero value of ft implies that that particular ft is unavail- 

x,y * r x,y 

able) : 


(i) All four values of ft are available. Select the wave- 

x,y 

length pair (a value for x) as follows: Set x to 2. Reset 


XX X 

. X to 1 if, and only if, > s 2j10 qq> and S l,400 

j $2 ^qq. Then proceed to (v) below. 

(ii) -Only values of ft^ ^000 and ^1 400 are ava:i -l a ^l e * Set x to 

. 1, and proceed to (v) below. 

I I 

(iii) Only values of ^qqq an ^ ^ 2 400 are ava ^--'- a ^-*- e * Set x to 

2, and proceed to (v) below. 

(iv) None of (i) to (iii) above. It is not possible to^obtain a 
value of ft^. Set it to zero, and omit Step (v) below. 


(v) 


ft = ft innr . 
e x, 1000 


^e ^x, 400 


for R < 0. 2 , 

for R > 0.8 , 


(9) 

( 10 ) 


and for 0.2 < R < 0.8, the ft £ is given by 


Q 

e 0.6 x, 1000 


R - 0.2 I 
0.6 x,400 


(ID 


15 



The foregoing description of the procedure is rather mechanical. 

We will now attempt to provide some rationale behind various steps of 
this basic procedure. As mentioned earlier, a value of P 0 for the 
place and time of observation is not available. It has therefore been 
thought appropriate to perform computations, at all times, for two pres- 
sure levels, the 1,000 mb value corresponding to an average sea— level 
atmospheric pressure, and the 400 mb value corresponding (hopefully) to 
the top of high clouds encountered in the earth's atmosphere most of the 
times. 

Chances of encountering high reflectivities at high surface pres- 
sures are very small unless there is snow on the ground, or there are 
very low-level, optically thick clouds. Thus, a low value of R can be 
associated with a high value of Pq, and vice versa. Hence, Eqs. (9) to 
(11) represent- a heuristic but ingenious attempt to relate R to Pq. 

(In the original BUV, total-ozone .estimation procedure, it is recognized 
that higher values of R can frequently occur for higher values of P 0 , 
especially at high latitudes in the winter months. For such cases, an 
additional special step is set up in that program.) 

An assumption is made that the computed effective albedo is inde- 
pendent of wavelength. However, we are not certain of this assumption 
under all possible conditions. It is therefore prudent to measure 
effective albedo as close to the wavelength region of interest (0.3125- 
0.3398 ym) as possible. However, we have no information on ozone-amount 
and hence the wavelength for measuring effective albedo has to be 
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located outside the ozone absorption band, i.e., beyond 0.3650 pm. The 
nearest measurement for an effective albedo measurement is at 0.3800 pm 
which cannot be considered very close to the wavelength region of inter- 
est. Therefore, it is advisable to call the 0.3800 pm effective albedo 

measurement as the coarse effective albedo j use it to obtain a first 

G C 

guess of the total ozone value (Q ) , and use this to compute 

x,y x,y 

(hopefully) a better effective albedo value closest to the spectral re- 


gion of interest. This new effective albedo value called the improved 

effective albedo R* be better calculated for the region of the least 

A.,y 

absorption by ozone to minimize errors due to uncertainties in the value 


of . 

x,y 


Hence, the choice of the 0.3398 pm measurement for this 


purpose. 


A computer- program for estimating total ozone using the simulated 
BUV measurements and steps outlined above, is described in Appendix A. 


2.5 Differences from the Original Procedure: Some minor differences 

occurring between the original total-ozone estimation procedure (see 
§2.1), and the one used for estimating total ozone from the simulated 
BUV data (see §2.4), are described below for the sake of completeness. 

The basic tables described in §2.2 are prepared in the original 
version of the procedure for an atmosphere with 1,013 mb surface pres- 
sure, and for the same atmosphere chopped off at a level where the pres- 
sure is four-tenths of 1,013 mb. These two sets are then respectively 
referred to as 1.0 atmosphere table, and 0.4 atmosphere table. This 
difference is’ of no significant importance in our analysis. 


17 



Because of changes in the amounts and profile characteristics of 
ozone encountered from the equator to pole, the original procedure calls 
for three subsets of tables for each Pq value. The subset used depends 
upon the geographic location of the area under observation; i.e., low 
latitude, midlatitude, or high latitude. Two additional subsets are 
also provided for use under special circumstances. A migration is gen- 
erally permitted from one subset to another whenever the £2^-range of 
the subset selected at the first time, is found to be insufficient. The 
subset to which a migration is made depends upon several factors such as 

the nature of the original subset, the £2. —direction along which more 

xn 

tables are required, and the exact- location of the site of observation. 

C X 

If a value of £2 and/or £2 cannot be determined after one such 
x,y x,y 

migration from 'one subset to another (a special case can arise where 

2 

even one migration cannot be permitted) , the corresponding £2 is - 
labeled as undeterminable. We have not incorporated this sophisticated 
maneuver in our procedure as our interest is primarily restricted to the 
analysis of problems arising from several ingenious but heuristic paths 
taken to hopefully minimize the effects due to unavailability of values 
of the basic parameters Pq and R , and those due to the breaking down 
of various assumptions (e.g., independent of X; models free of aero- 

sols and water droplets; horizontal homogeneity; and Lambert law of 
reflection) in real life. 

The original procedure provides a subprocedure for an interpolation 
of tables between two successive values of Gq for which the basic 
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tables are available. This is because real observations can exist for 
any arbitrary value of 0q. Our simulated observations are only made 
for those values of 0 q for which the basic tables are available. 

The specific BUV configuration results in measurements of I (X, ©o) 

at X = 0.3125, 0.3175, 0.3312, and 0.3398 pm, not simultaneously, but at 

a time interval of about 2 sec. In order to correct for changes in the 

effective albedo due to the translation of the satellite, 1^(0.3800, 0 O ) 

is measured four times, one coinciding with measurement of each of the 

four wavelengths in the region of ozone absorption. This additional 

C I 

complication results in different values of R-. and R at different 

A,y 

wavelengths. The effective albedo - R given by Eq. (8) is then redefined 
to include weighting for R^ values for the pair selected in Step IX. 
It should be pointed out that this use of different coarse or improved 
effective albedoes at different wavelengths is in no way connected with 
a real possible dependence of on X. 

In the original procedure, a quadratic interpolation procedure is 
used for computation of Rq 3398 y in 

Finally, the basic tables in the original procedure are prepared 
after taking into account polarization of the scattered radiation when 
a given atmospheric model is illuminated by Fq-^ units of solar energy 
passing normally, in unit time, through a unit area located at the top 
of the atmosphere. Thus, as the incident energy is wavelength depen- 
dent, the term inside the large parentheses in Eq. (4) or (5) has to be 
multiplied by the term Fq^/Fo-^. The atmospheric models in the orig- 
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inal' procedure extend toa height of 80 km, while those used by us ex- 
tend to a height of 70 km, only. There is also some difference in the 
manner in which the basic models are subdivided into layers. 


III. INPUT DATA 

3.1 Atmospheric Models: The atmospheric models with a surface pres- 

sure (P 0 ) of 1,000 mb are divided into 32 layers, and those with Po = 
400 mb are divided into 25 layers. Values of the geometric thickness in 
km, pressure thickness in mb, and ozone amounts in atm-cm for each of 
the 32 layers of the 1,000 mb models are given in Tables I and II. Only 
the ozone amount of a layer changes from one model to the other. Ozone 
amounts in different layers of the first four models (viz, models with 
a total ozone amount £2. given by £2. = 0.200, 0.250, 0.300, and 

0.350 atm-cm) are given in Table I. Table II contains values of h£2 

for each of the 32 layers for modeis with £2. = 0.400, 0.450,' 0.500, 

xn 

0.550, 0.600, and 0.650 atm-cm. 

Total ozone amounts mentioned in the preceding paragraph are for 

the 1,000 mb models. As mentioned earlier (§2.2), the 400 mb models are 

obtained after deleting the bottom-most seven layers of the 1,000 mb 

models. Thus, £2 = £2. for the models with Po = 1,000 mb, but £2 < £2. 

in • u xn 

for the models with Pq < 1,000 mb. Table III contains information 
about P q , £2, and £2^ for 22 different models for which computations 
are performed. 
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Table I. Values of the geometric thickness in km, pressure thickness in 
mb, and ozone amount in atm-cm for various layers of the model 
No. 2, 3, 4, and 5 with 0.200, 0.250, 0.300, and 0.350 atm-cm 
of total ozone, respectively. Geometric and pressure thick- 
ness remain unchanged from one model to the other. 


Layer 

Thickness 

Amount of ozone (atm-cm) in the layer 
for the model 

number 

Geometric 

(km) 

Pressure 

(mb) 

No. 2 

No. 3 

No. 4 

No. 5 

1 

10.0 

0.17 

0.00002 

0.00002 

0.00002 

0.00002 

2 

10.0 

0.74 

0.00016 

0.00016 

0.00016 

0.00016 

3 

5.0 

0.80 

0.00290 

0.00290 

0.00290 

0.00290 

4 

5.0 

1.53 

0.00722 

0.00722 

0.00725 

0.00732 

5 

5.0 

3.16 

0.01735 

0.01750 

0.01769 

0.01768 

6 

5.0 

6.55 

0.03630 

0.03340 

0.03792 

0.03925 

7 

5.0 

13.75 

0.05376 

0.06650 

0.06900 

0.07050 

8 

1.0 

4.30 

0.01100 

0.01500 

0.01750 

0.01850 

9 

1.0 

5.10 

0.01040 

0.01460 

0.01830 

0.01940 

10 

1.0 

6.00 

0.00975 

0.01360 

0.01710 

0.01990 

11 

1.0 

7.20 

0.00890 

0.01220 

0.01590 

0.01940 

12 

1.0 

8,30 

0.00800 

0.01120 

0.01472 

' 

0.01820 

13 

1.0 

9.60 

0.00690 

0.00920 

0.01330 

0.01650 

14 

1.0 

11.20 

0.00510 

0.00710 

0.01120 

0.01390 

15 

1.0 

12.60 

0.00360 

0.00510 

0.00845 

0.01180 

16 

1.0 

14.50 

0.00195 

0.00400 

0.00640 

0.00920 


. . .to be continued 
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Table Ijcont'd.' 


Layer 

Thickness 

Amount of ozone (atm-cm) in the layer 
for the model 

number 

Geometric 

(km) 

Pressure 

(mb) 

No. 2 

No. 3 

No. 4 

No. 5 

17 

1.0 

17.50 

0.00105 

0.00240 

0.00471 

0.00802 

- 18 

1.0 

21.50 

0.00090 

0.00120 

0.00380 

0.00720 

19 

1.0 

24.50 

0.00070 

0.00100 

0.00340 

0.00640 

20 

TTo 

27.50 

0.00050 

0.00080 

0.-00260 

0.00590 

21 

1.0 

31.50 

0.00040 

0.00070 

0.00210 

■ 0.00510 ‘ 

22 

1.0 

35.50 

0.00031 

0.00065 

0.00195 

0.00440 

23 

l.o- • 

40.50 

0.00035 

‘ 0.00065 

0.00186 

0.00340 

' . 24 

1.0 

45.00 

0. G0056 

0.00080 

0.00180 

0.00250 

25 

. i.o 

51.00 

0.00075 

0.00125 ' 

0.00181 

0.00200 

. 26 _ 

1.0 

61.00 

0.00105 

0.00150 

0.00194 

0.00220 

27 

1..0 

69.00 

0.00125 

0.00190 

0.00225 

0.00260 

28 

1.0 

75.00 

0.00167 

0.00230 

0.00248 

0.00300 

29 

1.0 

87.00 

0.00180 

0. 00245 

0.00289 

0.00310 

30 

1.0 

97.00 

0.00190 

0.00265 

0.00294 

0.00320 

31 

1.0 

100.00 

0.00180 

0.00265 • 

0.00286 

0.00320 

32 

1.0 

111.00 

0.00170 

0.00240 

0.00280 

0.00315 
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Table II. Values of the ozone amount in atm-cm for various layers of 

the model no. 6, 7, 8, 9, 10, and 11 with Q = 0.400, 0.450, 
0.500, 0.550, 0.600, and 0.650 atm-cm of total ozone, respec- 
tively. 


Layer 

Amount of ozone 

(atm-cm) in the layer for the model 

number 

No. 6 

No. 7 

No. 8 

No. 9 

No. 10 

No. 11 ' 

1 

0.00002 

0.00002 

0.00002 

0.00002 

0.00002 

0.00002 

2 

0.00016 

0.00016 

0.00016 

0.00016 

0.00016 

0.00016 

3 

0.00290 

0.00290 

0.00290 

0.00290 

0.00290 

0.00290 

4 

0.00732 

0.00742 

0.00742 

0.00742 

0.00742 

0.00742 

5 

0.01810 

0.01780 

0,01790 

0.01790 

0.01790 

0.01790 

6 ' 

0.03905 

0.03950 

0.03950 

0.03950 

0.03950 

0.03950 

7 

0.07290 

0.07380 

0.07390 

0.07390 

0.07390 

0.07390 

8 

0.01980 

0.01990 

0.02000 

0.02000 

0.02000 

0.02000 

. 9 

0.02170 

0.02220 

0.02280 

0.02340 

0.02400 

0.02460 

10 

0.02340 

0.02440 

0.02570 

0.02700 

0.02830 

0.02960 

11 

0.02120 

0.02380 

0.02700 

0.03020 

0.03340 

0.03660 

12 

0.02020 

0.02290 

0.02500 

0.02710 

0.02920 

0.03130 

13 

0.01860 

0.02110 

0.02400 

0.02690 

0.02980 

0.03270 

14 

0.01720 

0.01880 

0.02200 

0^02520 

0.02840 

0.03160 

15 

0.01440 

0.01700 

0.01900 

0.02100 

0.02300 

0.02500 

16 

0.01220 

0.01550 

0.01800 

0.02050 

0.02300 

0.02550 


.. . .to be continued 
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Table II 3 coni' d. 


Layer 

Amount of ozone 

(atm-cm) in the layer for the model 

number 

No. 6 

No. 7 

00 

• 

o 

J3 

No. 9 

No. 10 

No. 11 

17 ‘ 

0.01150 

0.01490 

0.01800 

0.02110 

0.02420 

0.02730 

18 

0.01110 

0.01480 

0.01900 

0.02320 

0.02740 

0.03160 . 

19 

0.01070 

0.01500 

0.02000 

0.02500 

0.03000 

0.03500 

20 

0.01030 

0.01540 

0.02100 

0.02660 

0.03220 

0.03780 

21 

0.00950 

0.01430 

0.01900 

0.02370 

0.02840 

' 0.03310 

22 

0.00710 

0.01080 

0.01400 

0.01720 

0.02040 

0.02360 

- 23 

0.00520 

0.00710 

0.00900 

0.01090 

0.01280 

0.01470 

24 

0.00310 

0.00410 

0.00500 

0.00600 

0.00700 

0.00800 

-25 

0.00245 

0.00320 

0.00410 

0.00500 

0.00590 

0.00680 

26 

' 0.00185 

0.00300 

0.00380 

0.00460 

■ 0.00540 

0.00620 . 

27 

0.00195 

6.00290 

O'. 00350 

0.00410 

0.00470 

0.00530 

28 

0.00250 

0.00310 

0.00350 

0.00390 

0.00430 

■ 0.00470 

29 • 

0.00310 ■ 

0.00330 

0.00350 

0.00370 

0.00390 

0.00410 

30 

0.00360 

0.00360 

0.00370 

0.00390 

0.00410 

0.00430 ' 

31 

0.00345 

0.00370 

0.00390 

0.00410 

0.00430 

0.00450 

32 

0.00345 

0.00360 

0.00370 

0.00390 

0.00410 

0.00430 
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Table III. Surface pressure (Pg in mb) , actual total ozone 

amount (ft in atm-cm) , and the total ozone amount 

when the model is extended to the mean sea- level 

(ft. in atm-cm) for various atmospheric models. 
m 


Model number 

Pg in mb 

ft in atm-cm 

ft. in atm-cm 
xn 

1 


0.0 

0.0 

2 



0.200 

3 


0.250 

0.250 

4 



0.300 

5 

1 BE 

0.350 

0.350 

' - 6 

1 

0.400 

0.400 

7 


’ 0.450 

0.450 

8 

Be 

0.500 

0.500 

9 

1,000 

0.550 

0.550 

10 

1,000 

0.600 

0.600 

11 

1,000 

0.650 

0.650 

12 

400 

0.0 

0.0 

13 • 

400 

0.189 

0.200 

14 

400 

0.234 

0.250 

15 

400 

0.282 

0.300 

16 

400 

0.330 

0.350 

■ 17 

400 

0.380 

0.400 

18 

400 

0.427 

0.450 

19 

400 

0.474 

0.500 

20 

400 

0.522 

0.550 

21 

400 

0.569 

0.600 

22 

400 

0.617 

0.650 
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3.2 Optical Bata: Table IV contains information on the normal Rayleigh 

scattering optical thickness of the standard atmosphere and 

absorption coefficient of one atm-cm of ozone to the base e for all 
six wavelengths at which computations are carried out. 


Table IV. Normal Rayleigh scattering optical thick- 
ness [T b (s ’ R) l of the standard atmosphere, 
and absorption coefficient (a) of one 
atm-cm of ozone to- the base e at six 


different wavelengths for which computa- 
tions are performed. 


Wavelength (ym) 

(s,R) 

b 

a 

0.3125 

1.0200 

1.6700 

0.3175 

, 0.9570 

0.9100 

0.3312 

0.8000 

0.1750 

0.3398 

0.7180 

0.0482 

0.3600 

0.5634 

0.0012 

0.3800 

0.4494 

0.0000 


IV. SOME CHARACTERISTICS OF THE QUANTITY N c 

4.1 Variations of N with : In §2.3, it was mentioned that the 

J c -tn 

quantity N (X^, Xjj R^., R^.; Pqj generally increases with an 
e l j 

increase in the total ozone content of the atmospheric model when 

other parameters appearing within the parentheses are held constant.- 
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curve then, forms the basis for 


This characteristic of the N vs . ft . 

c xn 

estimation of total ozone from the measured value of this quantity, 

i.e. , N (X., X., 0 q) given by Eq. (5). We shall therefore examine this 
m x 3 

property of the vs. ft^ n curves, in some detail, in this subsection. 

Variations of N as a function of the total ozone content (ft. ) 
c xn 

in the atmospheric -model are shown in Figs. 1 and 2 for the first 
(0.3312-0.3125 ym) , and the second (0.3398- 0.3175 ym) pair, respective- 
ly for p 0 = 1,000 mb, and R^. = . =0.0. Different curves are for 

different values of the parameter G 0 . For the first pair, the quantity 

N increases with an increase in ft. for values of 0 O up to 75.6°. 

c m 

The 0 n ■=. 79.6° curve of N vs. ft. of Fig. 1 shows only a very small 
u , c xn 

increase in N as ft. is increased from 0.550 to 0.600 atm-cm (the 
c ■ in , 

corresponding values of are 100.33 and 100.70 units, respectively). 

A further increase in ft. from 0.600 to 0.650 atm-cm results in a very 

in 

small decrease in N from 100.70 to 100.68 units. For still higher 

c 

values of 0 n» the N vs. ft. curve depicts a clear maximum which shifts 
u c in 

from 0.500, to 0.400, to 0.350 atm-cm as 0g is increased from 82.5°, 

to 84.7°, to 86.7°. The 8 0 = 90° curve shows a continuous decrease in - 

N as ft. is increased from 0.200 atm— cm to 0.650 atm-cm. For the 
c m 

second pair (Fig. 2), the quantity N c increases with an increase in 

ft. at all ten 0 n s-. However, the curves for 0 q = 86.7° and 90.0° 
in 

are rather flat for high values of ft. . As for example, N (90°, 0.600 

. m c 

atm-cm) = 77.11 units, and 1^(90°, 0.650 atm-cm) = 77.22 units. The 

curves of N vs. ft. for the first pair are generally more sensitive 
cm 




0.200 0.300 0.400 0.500 0.600 

Ozone 'amount (atm-cm) 


Fig. 1 Variations of the computed quantity [see Eq. (4)] 

as a function of the total ozone amount (Q. ) in the 

m 

atmospheric model. Different curves are for different 
values of the parameter 0 q, solar zenith angle.. 

X^ = 0.3312 pm, X = 0.3125 pm, R^ = = 0.0, and 

P 0 = 1,000 mb. 
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: ioo 



Fig. 2 Variations of the computed quantity [see Eq. (4)] 

as a function of the total ozone amount (ft. ) in the 

m 

atmospheric model. Different curves are for different 
values of the parameter 0 qj solar zenith angle. 

A. = 0.3398 ym, A. = 0.3175 ym, Rx- = % . =0.0, and 
Pq = 1,000 mb. 
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to a change in for the small values of 9 q, than the corresponding 

curves for the second pair. As for example, an increase in total ozone 

from 0.200 atm-cm to 0.250 atm-cm results in an increase in N (0.3312, 

c 

0.3125, 0°) from 16.80 to 22.72 units, but an increase in N^CO.3398, 
0.3175, 0°) from 4.57 to 8.10 units, only. 

The higher sensitivity of the vs. curves (to changes in 

) for the first pair at small 9 q values compared to that of the 
corresponding curves for the second pair, is due to a relatively stronger 
absorption by ozone for the spectral regions of the first pair compared 
to that of the second pair. The loss of sensitivity of the first pair 
at large 0q to changes in ozone content, is primarily due to a rise in 
the height of the effective scattering layer through the ozone layer. 
Further information on this particular aspect can.be found in a paper by 
Dave and Mateer [1967]. 


Results similar to those presented in Figs. 1 and 2 but for f 

R^ . — 1.0 were also analyzed for both the wavelength pairs. However, 
no graphical results are presented for these cases in this report. An 


increase in the surface reflectivity resulted in a very significant 


change in the slopes of the N vs. Q. curves at large 9g for the 

c xn 

first pair. The 0 q = 79.6° N vs. £2. curve for R-\ . = R\. = 1.0 
r u c xn 13 . 

did not show any decrease in N with increase in £2. . The maxima of 

c in 

the N vs. 0. curves for 0 O = 82.5°, 84.7°, and 86.7° shifted to the 
c xn 

right for the first pair as the reflectivity of the underlying surface 
was increased from 0.0 to 1.0. The 0 q = 90° curve for the first pair 
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with . = 0.0 (Fig. 1) shows a continuous decrease in N with 

i J c 

increase in ft, , while the same curve but for Rt, . = R-; = 1.0 showed 
xn A i J 

a weak maximum around ft^ n = 0.500 atm-cm. Similar improvements in the 

sensitivity of N to changes in Ji. at large 0ns due to an increase 
c xn 

in reflectivity were also seen for the second pair after comparing re- 
sults in Figs. 2, with those for R^ . = R ( \. - 1.0, and Pq = 1,000 mb. 

From the preceding discussion on the trends of N vs. ft. varia- 

tions, it is clear that the first wavelength pair loses its sensitivity 

completely to changes in ft^, f° r at 0o = 82.5° when R^ = 

R^. = 1.0, but at 0 Q = 79.6° when R^. = R^. = 0.0. The arbitrary 

method of determining the coarse, or the ■improved, effective albedo 

described in §2.4, results in a requirement of generation of the N c vs. 

ft curves with negative values of R^. To this - effect, it should be 

pointed out that for R^ . = R^. = -0.1, the quantity N (0.3312, 0.3125, 

— 3 

ft. , 1,000 mb, 75.6°) is equal to 92.33, 94.70, 96.36, and 97.46 units, 
xn 

for ft^ = 0.500, 0.550, 0.600, and 0.650 atm-cm, respectively. 

Trends in variations of N with ft. for the Pn = 400 mb models 

cm u 

are, in general, very similar to those discussed in the preceding para- 
graphs for corresponding models with Pq - 1,000 mb. Representative 
results for the first and the second pairs are shown in Figs. 3 and 4, 
respectively for Pg = 400 mb, and R^ = R^ = 0.0. 

In Fig. 5, we have plotted values of 11^(0.3398, 0.3175) and 
13^(0.3398, 0.3312) as a function of ft_^, for = = 0*0, = 

86.7° and 90°, and P 0 = 1,000 and 400 mb. The results for the new pair 
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0.200 0.300 0.400 0.500 0.600 

Ozone amount (atm-cm) 


Fig. 3 Variations of the computed quantity N c [see Eq. (4)] 

as a function of the total ozone amount (fl. ) in the 

xn 

atmospheric model. Different curves are for different 

values of the parameter 6q, solar zenith angle. - 

X. = 0.3312 pm, X. = 0.3125 pm, Rt . = R>. =0.0, and 
x - j A x 3 

Pg = 400 mb. 
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0.200 0.300 0.400 0.500 0.600 

Ozone amount (atm-cm) 


Eig. 4 Variations, of the computed quantity [see Eq. (4)] 

as a function of the total ozone amount (Q. ) in the 

m 

atmospheric model. Different curves are for different 
values of the parameter Qq, solar zenith angle. ' . 

= 0.3398 pm, = 0.3175 ym, R;^ = R^ = 0.0, and 
P 0 = 400 mb. 
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'Fig. 5 Variations of the computed quantity [see Eq. (4)] 

as a function of the total ozone amount (£2. ) in the 

.in 

atmospheric model. Solid curves are for the 0.3398-- 
0.3175 ym pair, and the broken curves are for the 
0.3398- 0.3312 ym pair. R A± = = 0.0. Only the 

. curves for 0q = 86.7° and 90° are shown. Pq = 1,000 
mb as well as 400 mb. 
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(viz, third pair; 0.3398- 0.3312 ym) are shown, by broken curves, while 

those for the second pair (viz, 0.3398- 0.3175 pm) are shown by solid 

curves. A comparison of the curve number m with the curve number 

m + 4 suggests a small but significant increase in the slope of the 

curve as we transpose from the second, to the third pair. 

Thus, in principle, there is a small advantage in use of the third pair 

when Q . and 0n assume large values. However, there are several 
m u 

additional factors such as the rapid changes in brightness within the 
field of view and the reliability of the 0g value at the place and 
time of observation, which may control the desire of estimating total 
ozone in such extreme cases. 


Values of 11^(0.3398, 0.3175, 84.7°) are plotted in Fig. 6 as a func- 
tion of the total ozone (^ n in- atm-cm) for Pg = l,000mb (solid curves), 
and Po=4QOmb (broken curves). There are three sets of curves, one set 
for each of the following values of the ground reflectivity (R^-R^. = R) : 

0.0,. 0.5, and 1.0. This diagram is provided primarily for realizing 

Z ' 

the spread between &2 y f° r a given measurement, N m « As for example, 

Z 

for R^ . = R\. - 0.0, a value of N of 60 units yields £2- i nnn = 
i j m z j i uuu 

_ 2 

0.238 atm-cm, and = 0.272 atm-cm. On the other hand, for = 

z z 

R, . = 0.5, we have , nnn = 0.212 atm-cm, and = 0.220 atm-cm. 

Aj 2,1000 2,400 

If we enter this diagram with = 80 units and R^ = R^_. = 0.0, we 

z z 

will come out with a value of 0.423 atm-cm 'for but /r ._ 

2,1000 2,400 

will be undeterminable. It may be necessary to construct more diagrams 
like this one for a better understanding of the results obtained with 
the total-ozone estimation procedure described in §11. 
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Fig. 6 Variations of the computed quantity [see. Eq. (4)] 

as a function of the total ozone amount (Q. ) in the 

in 

atmospheric model for a solar zenith angle of 84.7°. 
X. = 0.3398 ym, and X. = 0.3175 ym. The solid 
(broken) curves are for the models with 1,000 mb 
(400 mb) surface pressure. Values of Rj^ = = R 

are as shown in the diagram. 






4.2 Values of as Obtained with Different Procedures: In §2.5, we 

pointed out several differences between computational procedures used 
for evaluation of the quantity N^A^ ; R^, R^_. ; ^ in , Pq, e o) f° r the 
original NASA/GSFC total-ozone estimation procedure for the analysis of 
the actual BUV data, and for that used by us in the analysis of the sim- 
ulated observations. It was therefore considered desirable to compare 
values of this quantity as obtained from these two procedures. 

Values of N c (A^, A j ; R^, r a j » %n» R 0> f° r fc h e following param- 

eters were computed using the NASA/GSFC procedure by Dr. Ashok Kaveeshwar 
associated with the BUV ozone-estimation project, and using the procedure 
for the analysis of the simulated BUV data (as described in this report) 
by the author: 0.3312- 0.3125, 0.3398- 0.3175, and 0.3398 - 0.3312 pm 

wavelength pairs; R^. = R^ . = R given by R = - 0. 1 (0. 1) 1, 0; ten values 

i J 

of 0 q listed earlier; Pq = 1,000 mb and 400 mb for the author's calcu- 
lations, i.e., 1.0 and 0.4 atmospheres for the NASA/GSFC work (see the 
second paragraph of §2.5); and = 0.200 (0.050) 0.450 atm-cm. These 

sets of tables of N^ were then exchanged for independent intercompari- 
sons. 

The results of these independent intercomparisons showed that values 
of the quantity N^ as obtained from these two independent procedures 
described above, agree within ±0.5 units on the average for cases with 
the solar zenith angle (0 q) not exceeding about 70°. For still larger 
values of 0 q, this difference increases to about ± 2 units. Consider- 
ing several differences between two procedures, it was decided by us 
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(Drs. Mateer, Novak, and Kaveeshwar, as well as the author) that the 
agreement between these two sets of N £ is very good, indeed. Further- 
more, it should be added that values of I , T , etc. , for the analysis 

c c 1 

of the simulated observations are to be computed using the same proce- 
dure. Thus, the aforementioned procedure dependence of cannot be 

expected to have any effect on the results of investigations carried 
out under this contract. 

The above-mentioned procedure dependence of can be due to one, 

or both, of the following causes: 

(1) Difference in vertical profiles of ozone as used by two procedures. 
The author has based ozone profiles for his work on the data used 
in the NASA/GSFC work. (These data were supplied to the author, 

in the numerical form, by Dr. C. L. Mateer.) However, because of 
the manner in which an atmospheric model is divided into layers by 
two procedures, it is possible that some minor differences have 
' popped in. 

(2) Neglect of polarization of the atmospheric radiation: Values 1^ 

for N for the NASA/GSFC procedure were obtained after taking into 
c 

account polarization of the scattered radiation. However, as men- 
tioned earlier, this aspect of atmospheric scattering was neglected 
in the investigation of the analysis of the simulated data. Dr. R. 
S. Fraser has shown that a significant portion of the above- 
mentioned difference is due to the neglect of polarization of the 
scattered radiation. 
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V: SIMULATION FOR THE BUV CONFIGURATION 


5.1 Results for the R-Independent-of-X Cases: Values of the effec- 
tive albedo [R given by Eq. (8)] and the best ozone estimate given 

by Eq. (9), or (10), or (11)] were evaluated for the atmospheric models 
no. 2 to 11, and 13 to 22 (see Table III) for ten different values of 
the solar zenith angle (0q — 0°, 45°, 60°, 70°, 75.6°, 79.6°, 82.5°, 
84.7°, 86.7°, and 90°), and for eight different values of the reflectiv- 
ity (R = 0.0, 0.1, 0.2, -0.3, 0.4, 0.6, 0.8, and 1.0) of the Lambert sur- 
face underlying the model. Results presented in this section (§5.1) as 
well as in the next section (§5.2) are for the cases where the quantity 
R is assumed to remain constant with changes in the wavelength (A) . 

Four- sets of programs were written for this purpose. They are re- 
ferred to as SITC1, SITC2, SITC3, and SITC4, or— for brevity— as Cl, C2, 
C3, and C4. (It may be noted that programs SITAA and SITBB described in 
the first technical report are for computing values of I c , T^, etc., 
for a- given model.) 

The program SITC1 faithfully follows the steps described in §2.4. 

As mentioned earlier, this basic procedure of §2.4 is expected to be 
identical to the one being used at NASA/GSFC for the estimation of total 
ozone from the BUV data. Some minor differences between the original 
NASA/GSFC procedure and the one used by us for the analysis of the sim- 
ulated observations, are listed in §2.5. 
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a 


The program SITC2 is a copy of the SITC1 program modified to set 


1,400’ ^1,1000 of Step 111 ’ and ^1,400 and ^1,1000 


of Step VI, to 


zero whenever 9 0 exceeds 79.6°. This is to avoid the possibility of 

- -- - • ■ - ' 2 " ' ■ " ’ 

obtaining erroneous values of ft , arising from peculiar properties 

x,y 

of vs. £2 curves discussed in §4.1. In an oral presentation of the 

results discussed in this §V to Dr. Mateer on September 20, 1976, we 

found that the original NASA/GSFC procedure contains an additional step 

which can be interpreted to carry out the same function (i.e., of setting 
2 

£2., to zero when 6n is large) in an indirect manner. We therefore 

i,y . 

feel that, based on this additional information, the SITC2 program more 
closely resembles the original NASA/GSFC procedure than the SITCl pro- 
gram. 


The SITC3 program is a copy of the SITC2 program modified to set 


£2 to zero whenever 
x,y 


ft 1 - £2° 
x,y x,y 


2: 0.030 atm-cm. This criterion is 

C 


based on the following argument: As mentioned in §2.4, £2^ is 

obtained by making use of the coarse effective albedo Rq yj and 

£2^ is obtained by making use of 'the improved effective albedo 
x,y 

I ■ C ‘ I 

R n rjno • A large absolute difference between £2 and £2 implies 

x 9 y x,y 

either a large difference between two computed values of the albedos, 

or alternately, a strong sensitivity of the vs. £2 curve to small 

changes in R. In either case, there is sufficient ground to discard 

the corresponding value of £2* . 

x,y 

The SITC4 program is a copy of the SITC3 program modified to accept 

-j* . _ 

the appropriate value of £2^ as £2 & provided R < 0.2 [see Eq. (9)], 
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and the appropriate value of ft* ^ 0 ^ as ft £ provided R s 6.8 [see Eq. 
(10)]. This modification is introduced to free the procedure of an un- 
necessarily strict acceptance criterion. When R '< 0.2 (R > 0.8), a 
value of ft* (ft* igOO^ * s not re 9 u ^ rec * computations of ft^; see 
Step IX of §2.4. 

The purpose of this simulation is to compare the value of ft g with 
the corresponding total ozone amount in the atmospheric model when its 
lower boundary is extended right down to the mean sea-level, i. e., the 
quantity ft^,. and to compare R with the corresponding value of R. 
Results of this latter comparison will be presented in §5.2. It is 
clear that this comparison task is by no means a trivial one (4 proce- 
dures x 20 models x 10 values of 0 q x 8 values of R) , and hence a 
good degree of discretion is used in the presentation of results. For 
convenience, we will refer to the quantity 1,000 x (ft - ft^) as the 
deviation.; A. A positive (negative) value of the quantity A implies 
an overestimation (underestimation) of total ozone by the procedure. 

Values of the quantity A are tabulated in Tables V, VI, and VII 
for the atmospheric models no. 3, 6 , and 9, respectively. These atmo- 
spheric models are for a surface pressure (Pq) of 1,000 mb. Their re- 
spective 'total ozone contents (fi^ ) are 0.250, 0.400, and 0.550 atm-cm. 
For each case, results are presented for all ten values of 0 q, but for 
only four values of R, viz, 0.0, 0.2, 0.6, and 1.0. Similar results 
for the atmospheric models no. 14, 17, and 20 [respective ozone amounts 
(ft. ) 0.250, 0.400, and 0.550 atm-cm when corrected for the ozone con- 
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Table V. Values of the deviation [A = 1,000 x (fi -ft. )] for the atmo- 

e xn 

spheric model no. 3; ft. = 0.250 atm-cm; P 0 = 1,000 mb. 


Q q in 
deg. 


Deviation A 

obtained with 

the procedure 


Cl 

C2 

C3 

C4 

1 

Cl 

C2 

C3 

C4 

1 


R = 

0.0 



R = 

0.2 


0.0 

0 

0 

0 

0 


B 

-1 

-1 

B 

45.0 

0 

0 

0 

0 



-1 

-1 

im 

" 60.0 

0 

0 

0“ 

0 


IS 

-1 

-1 

-1 

70.0 

0 

0 

0 

0 


0 

0 

0 

0 

75.6 

0 

0 

0 

0 


0 

0 

0 

0 

79.6 

0 

0 

-o 

0 


0 

0 

0 

0 

82. 5 

0 

0 

0 

0 


0 

0 

0 

0 

84.7 

0 

0 

0 

0 


0 

0 

0 

0 

86.7 

0 

0. 

0 

0 


0 . 

0 

0 

0 

90.0 

0 

0 

- 

0 


-1 

-1 

-1 

-1 



. R = 

0.6 



R » 

1.0 


0.0 

11 

11 

11 - 

11 


28 


28 

28 - 

45.0 

.9 

9 

.9 

9 


. 25 


. 25 

25 

60.0 

8 

8 

8 

8 


21 

21 

21 

21 

70.0 ' 

6 ' ‘ 

6 

' 6 

6 


16 

16 

16 

16 

75.6 

5 

5 

5 

5 


12 

12 

12 

12 

79.6 

4 

4 

4 

4 


10 

10 

10 

10 

82.5 

4 

3 

3 

3 


10 

7 

7 

7 

84.7 

3 

3 

3 

3 


• 8 

8 

8 

8 

86.7 

3 

3 

3 

3 


9 

9 

9 

9 

90.0 

5 

'5 

5 

5 


18 

18 

18 

18 


























Table VI. Values of the deviation [A = 1,000 x (ft -ft. )] for the atmo- 

g in 

spheric model no. 6; ft. = 0.400 atm-cm; Pq- = 1,000 mb. 


9 0 in 
deg. 


Deviation A 

obtained with 

the procedure 


Cl 

C2 

C3. 

C4 


Cl 

C2 

C3 

C4 



R = 

0,0 



R = 

0.2 











0.0 

0 

0 

0 

0 


0 

0 

0 

0 

45.0 . 

0 

0 

0 

0 


0 

0 

0 

0 

60.0 

0 

0 

0 

0 


0 

0 

0 

0 

70.0 

0 

0 

0 

0 


0 

o 

0 

0 

75.6 

0 

0 

0 

0 


-1 

-1 

- 1 . 

-1 

79.6 

0 

0 

0 

0 


0 

0 

0 

0 

82.5 

0 

0 

0 

0 


0 

0 

0 

0 

" 84.7 

0 

0 

0 * 

0 


0 

' 0 

0 

0 

86.7 

0 

0 

- 

0 


0 

0 

0 

0 

90.0 

- 

-■ 

- 

0 


-4 

-4 

- 

- 



R - 

0.6 



R = 

1.0 


- 0.0 

15 

15 

15 

15 


33 

33 

33 

33 

45.0 

13 

13 

13. ' 

13 


30 ■ 

30 

30 

' 30 ' 

60.0 

11 ' 

11 

11 

11 


- 26 

26 ' 

• 26’ 

26 

70.0 

9 

9 

9 

9 


20 

20 

20 

20 

75.6 

6 

6 

6 

6 


16 

16 

16 

16 

79.6 

5 

5 

5 

5 


12 

' 12 

12 

12 

82.5 

5 

5 

5 

5 


11 

11 

11 

11 

84.7 

6 

6 

6 

6 


13 

13 

13 

13 

86.7 

8 

8 

8 

8 


18 

18 

18 

18 

90.0 

8 

8 

8 

8 


24 

24 

24 

24 
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Table VII. Values of the deviation [A = 1,000 * (ft -ft. )3 for the atmo- 

e m 

spheric model no. 9; ft. = 0.550 atm-cm; Pq = 1,000 mb. 


6 0 in 

Deviation 

deg. | 

1 C1 

C2 

C3 


R =0.0 


A obtained with the procedure 



Cl 

C2 

; C3 

C4 


. R = 

0.2 



0.0 

0 

0 

- 

0 


0 

0 

0 

0 

45.0 - 

0 

0 

- 

0 


1 

mm 

1 

1 ' 

60.0 

- 0 

0 

- 

0 


1 


1 

| 

70.0 

0 

0 

- 

0 


-1 

H 

-1 

-1 

75.6 

0 

0 

- 

0 


0 

0 

0 

0 ' 

79.6 

0 

0 

- 

0 


0 

0 

0 

0 

82.5 

0 

0 

- 

0 


0 

0 

0 

0 

84.7 

91 

0 

- 

0 ' 


67 

0 

- 


86.7 

-251 

- 

- 

0 


- 

-3 

- 

- 

90.0 

-181 

- 

- 

0 


-183 


- 

- 



R = 

0.6 



R = 

1.0 


0.0 

21 

21 

21 

21 




43 

43 

45.0 

19 

19 

19 

‘ 19 




40 

40 

60.0 

16 

16 

16 

16 




' 34 

34 

70. a 

12 _ 

12 

. 12 

12 


26 


26 

26 

75.6 

10 

10 

10 

10 


22 

22 

22 

22 

79.6 

9 

9 

9 

9 


17 

17 

17 

17 

82.5 

9 

9 

9 

9 


17 

17 

17 

17 

84.7 

11 

11 

11 

11 


21 

21 

21 

21 

86.7 


13 

13 

13 


28 

28 

28 

28 

90.0 

3 

3 

- 

- 


33 

33 

33 

33 
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tents of the deleted layers; see’ §3.1] with a surface pressure of 400 
mb are presented in Tables VIII, IX, and X, respectively. Presence of 
the symbol (dash) in the place of- a numerical value for A implies 

that it was not possible to obtain a value of £2^ for that particular 
case. 


From the values of A given in Table V for the atmospheric model 
no. 3 with R = 0.0, we find that the total ozone amount in this parti- 
cular case can generally be estimated very accurately [A = 1,000 (Q - 
£2^ n ) = 0] for all ten positions of the sun by making use of any of the 
four procedures (Cl through C4) described above. The only exception is- 
the subcase of 0 q = 90° with the procedure C3 where it is not possible 

to estimate £2 as £2^ could not be determined. We find that 

e 2,400 

R 0. 3800, 400 = ^2,400 = °* 270 atm-cm, r q. 3398, 400 = °* 221, and 


£2„ = 0.231 atm-cm. Thus, 

2,400 


X c 

^2,400 ~ ^2,400 


hence , £2 


2,400 


> 0.030 atm-cm, and 
is set to zero. However, the effective albedo R is < 


0.2 for this subcase,' and hence £22 ^qqq = 0*250 atm-cm is returned as 
a value 'of £2^ by the C4 procedure. For the atmospheric model no. 3 
resting on a Lambertian surface with R = 0.2, the deviation A is 
either 0, or -1, for all procedures and all values of 0 q. Hence, we 
can state the total— ozone estimation procedure is capable of returning 
Very accurate (A = ± 2) values of £2^ provided £2^ is small (0.200- 
0.300 atm-cm), Pq is very high (very close to 1,000 mb), and R £ 0.2. 
[Evidently, this last statement (and also similar generalized statements 
appearing later in this report) is made after examining results in 


45 



Table VIII. Values of the deviation [A = 1,000* (£2 -£2. )] for the atmo- 

e in 

spheric model no. 14; £2. = 0.250 atm-cm; Pq = 400 mb. 
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Table IX. Values of the deviation [A = 1,000 x (52 -52. )] for the atmo- 

e m 

spheric model no. 17; 52. =» 0.400 atm-cm; Pq = 400 mb. 
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Table X. Values of the deviation [A = 1,000 x (ft -ft. )] for the atmo- 

e xn 

spheric model no. 20; ft. = 0.550 atm-cm; Pq = 400 mb. 


9 0 in 
deg. 


Deviation A 

obtained with 

the procedure 


Cl 

C2 

C3 

C4 


Cl 

C2 

C3 

C4 

„ , 


R = 

0.0 




R = 

0.2 


0.0 

44 

44 

- 

- 


9 

9 

9 

9 

45.0 

25 

25 

- 

- 


6 

6 

- 

- 

60.0 

3 

3 

- 

- 


i 

1 

- 

- 

- 70.0 

94 

94 

- 

- 


30 

30 

- 


75.6 

64 

64 

- 

- 


17 

17 

- 

- 

.79.6 

53 

53 


- 


10 

10 

- 

- 

82.5 

83 

39 


- 


9 

9 

’ - 

- 

" 84.7 

-208 

82 

- 

- 



-20 

- 


86.7 

-315 

- 

- 

- 


.-280 

93 

- 

- 

90.0 

-225 

- 

- 

- 


-207 

- 

- 

-- 



R = 

0.6 



R = 

1.0 


BB 

-10 

-10 

-10. 

-10 


0 

0 

0 

0 

B 

-9 

-9 

-9 ' 

-9 


0 

0 

0 

o - 

m o 

' —7 

-7 

—7 

-7 


- 0 

o • 

0 

0 

- 7 °-°. 

-5 

-5 

-5 

-5 


0 

0 

0 

0 , 

75.6 

-4 ’ 

-4 

' -4 

-4 


0 

0 

0 

0 

79.6 

-4 

-4 

-4 

-4 


. 0 

0 

0 

0 

82.5 

-4 

-4 

-4 

-4 


0 

0 

0 

0 • 

84.7 

-5 

-5 

. -5 

-5 


0 

0 

0 

0 

86.7- 

- 

-6 

-6 

-6 


' 0 

0 

0 

0 

90.0 

-2 

-2 

- 

- 


0 

0 

0 

0 
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appropriate intermediate ranges . ] If we moderate our definition of 
very accurate from A = + 2 to A = ±4 units, the range can be ex- 

tended to about 0.350 atm-cm, and the R range to about 0.4. For 
R h 0.6, all total-ozone estimation procedures yield a significantly 
overestimated value of The degree of overestimation is independent 

of the procedure used but, it increases with an increase in R, and is 
very significant indeed, at small (Qq £ 45°) and large (0g ~ 90°) values 
of 0 g . This is because the estimation procedure is forced to accept 
^x 400 instead of 10 oo ^ see Eq * at high values of R due to . 

our lack of knowledge of the quantity Pg. It may be noted' that for the 


model no. 3 with 0 O = 0° and R = 1.0, all estimation procedures come up 

I .1 

with ft, , nnr , - 0.250 atm-cm, but = 0.278 atm-cm. Thus, we see 

1,1000 1,400 

a clear need of having some information about the parameter P 0 at the 


time of observation for increasing the accuracy of Sl^, especially if 
the sun is very close to the local zenith. 


• Numerical results presented in Tables VI and VII for the atmo- 
spheric models no. 6 and 9, respectively, support the analysis given in 
the preceding paragraph, in a general sense. For the model no. 6 with 
00 = 90° and R = 0.0, a value of could be obtained with the C4 

procedure. On the other hand, the restriction imposed by the C3 proce- 

. 1 

dure results in an inability to accept an otherwise good value of 

y 

for the 0 q = 90°, R = 0.2 subcase. For the model no. 9 (Table VII) 
with 0.550 atm-cm total ozone, the Cl procedure yields very erroneous 


values for ft when the sun is near the horizon, and the ground reflec- 
e 
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tivity is small. A value of 91 for the R = 0.0, 0 O = 84.7° subcase 

occurs because ^.3398,1000 and R 0. 3398, 400 are evaluated to be - 0. 091 

I I 

and 0.248, respectively. Corresponding values of i000 and ^2 400 
then turn out to be 0.641 and 0.535 atm-cm, respectively. A large nega- 
tive value for A for the R = 0.0, 8 q = 86.7° subcase (and also for 
the subcases R = 0.0, 0q = 90°, and R = 0.2, 0q = 90°) is due to the 
presence of a max imum in the vs. curve discussed in §4.1, asso- 

ciated with the existence of the condition of ^qqq being greater 

X 

than S 2 10 qq. in the region of immediate interest. Thus, we have sever- 
al examples stressing the need of using only the second wavelength pair 
for large values of 0q. As mentioned earlier, the NASA/GSFC procedure 
for estimation of total ozone from the BUV measurements, does contain a 
step to perform the same function in an indirect manner.' 

From the results presented in Tables VIII to X for atmospheric 
models with a surface pressure of 400 mb, we find that the deviation A 
vanishes for all subcases with R = 1.0. Results for the same models 
but with R = 0.8 for which no numerical data are presented in this 
report, also show A to be zero for all subcases. This is because the 
estimation procedure is forced. to select values of £2^ ^qq as S for 
R > 0.8,' and the surface pressure for these models also happens to be 
400 mb. For R = 0.6, the quantity A carries a negative value which 
increases with an increase in 6q, and also with a decrease in ^n‘ 

The maximum underestimation in total ozone is by about 2%. For still 
lower values of R (R = 0.2, and 0.0), we find the total ozone amount 
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to be overestimated by as much as 31% in some extreme cases (e.g., A = 

77 units for the R = 0.0, 0o = 0° subcase for the model no. 14) by the 
Cl and C2 procedures. Many of such bad cases can be discarded by impos- 
ing the extra constraint added in arriving at SITC3 from SITC2. It 
should be pointed out that, once in a while,, this extra constraint does 
result in discarding of a good case also; as for example, a use of Cl or 
C2 procedure yields a value of zero for A for the subcase R = 0.0, 

0q = 15 . 6 ° of the model no. 17, but a use of C3 or C4 procedure leads 
to the inability in estimation of £2 g . 

Somewhat detailed discussion in the preceding paragraphs can be 
summarized as follows: 

After exercising due precautions for not using the N £ vs. £2 curves 

exhibiting a maximum, two problem areas have developed due to a need 

for estimating total ozone in the absence of any information about 

■ the surface pressure. The first problem area consists of cases 

associated with values of the parameters Pq and R in the upper 

parts of their respective ranges. For such cases, we find that the 

estimated ozone amount (£2 ) can be greater than the input ozone 

e 

amount by as much as 0.045 atm-cm. The exact degree of overestima- 
tion depends upon values of the various parameters such as P 0 , R, 

0n, and £2. . 

. u xn 

•The second problem area consists of cases associated with values 
of the parameters Po and R in the lower part of their respective 
ranges. For such cases, we again find that the total ozone amount 
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can be overestimated by as much as 0.075 atm-cm. Some of these 


cases can be discarded by setting.- & to zero whenever 

x,y 


a 

x,y 


x,y 


2: 0.030 atm-cm for a given x,y combination (see discussion 


of SITC3 procedure). 


5.2 Effective Albedo: The total-ozone estimation procedure also com- 

putes another quantity called the effective albedo which is denoted by 
the symbol R [see Eq. (8)].* From the procedure outlined in §2.4 for 
the evaluation of R, it is clear that R is not related to the Lambert 
reflectivity (-R) of the surface underneath the atmosphere, only. Because 

of the use of tables for arbitrarily chosen surface pressures of 1,000 

C 

and 400 mb, and also because of the use of Q in evaluation of' 

x,y 

Rq 339Q R can also depend upon the surface pressure (Pq) of the 
model under investigation, and upon the solar zenith angle (6 q). How- 
ever, because of our knowledge that the Lambert's law is not valid for 
realistic surfaces, there seems to be a considerable temptation for 
relating R- to the optical characteristics of the surface only. The 
following discussion is provided to assist the reader in seeing some 
pseudoproperties of R which are in no way connected to the actual 
optical characteristics of the surface underlying the atmospheric model. 

c c 

Values of the coarse effective albedo (R Q ^goo 100 q and R Q 33 QQ ^qq) > 
improved effective albedo (R g 33gg 1Q00 and r q. 3398, 40(P’ as wel1 as of 
the effective albedo (R) as obtained during various stages of calcula- 
tions in the total-ozone estimation procedure are given in Table XI, for 
two arbitrarily selected (Nos. 4 and 18) atmospheric models with R = 0.0 
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Table XI. Values of albedo obtained during various stages of 
calculations in the total-ozone estimation proce- 
dure, for two arbitrarily chosen models. 


Model 

No. 

R 

e o 

R c 

0.3800,1000 

R I 

*0. 3398, 1000 

R c 

0.3800,400 

R I 

0.3398,400 

R 

4 

0.0 

45.0 

0.000 

0.000 

0.120 

0.188 

0.094 



70.0 

0.000 

0.000 

0.167 

0.243 

0.122 



90.0 

0.000 

0.000 

0.147 

0.226 

0.113 


1.0 

45.0 

1.000 

1.000 

1.006 

1.003 

1.002 



70. Q 

1.000 

1.000 

0.966 

0.970 

0.985 



90.0 

1.000 

1.000 

0.985 

0.997 

0.999 

18 

0.0 

45.0 

-0.168 

-0.333 

0.000 

0.000 

-0.167 



70.0 

-0.268 

-0.517 

0.000 

0.000 

-0.259 



90.0 

-0.223 

-0.487 

0.000 

0.000 

-0.244 


1.0 

45.0 

0.995 

0.997 

1.000 

1.000 

0.999 



70.0 

1.035 

1.031 

1.000 

1.000 

1.016 



90.0 

1.015 

1.004 

1.000 

1.000 

1.002 


and 1.0, and 0 O = 45°, 70°, and. 90° . For the model No. 4 with 1,000 mb 

C I 

surface pressure, values of Rq 3 qqq iqqq and r q .3398,1000 

the corresponding value of R as the surface pressures of the model and 

tables used are exactly the same. Such is also the case for values of 

c x 

R~ and R_ for the atmospheric model No. 18. However, 

0.3800,400 0.3398,400 

for the determination of Ro > 3800j400 and **.3398,400 for the model No. 
4, we enter the 400 mb tables with a relatively large value of 
I (0.3800, 0 n ) and I (0.3398, 0 n ), respectively. This step then results 

m u m v 
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in an over estimation of the improved effective albedoes. The degree of 

overestimation is naturally controlled by the slope of I (X , . , Pg, R. , 

c xn A 

9q) vs. R^ curves, and the separation of such curves for Pg = 1,000 

and 400 mb, in the region of immediate interest. Hence, we have pseudo- 
C X — 

dependence of R 0 . 3800,400’ R 0. 3398,400’ and R ° n the P arameter 9 0» for 
the atmospheric model No. 4. The above-mentioned arguments can also be 

used after appropriate changes to explain under-estimation and 9 0 - 

dependence of r£. 3800 , 1000 > R o. 3398, 1000 and R ValueS f ° r the atm °“ 
spheric model No. 18. 

A detailed analysis of R values for all 20 models confirmed the 
presence of the above-mentioned pseudo- dependence of R on 0 q, and on 
R, in all cases. Besides, a small, but significant, dependence of R 
on was also noted. We will not go into this aspect -any further. 

However, we will consider the variations of mean R (mean over 10 values 
of 0 g , and over 10 atmospheric models with the same surface pressure, 
but for a given value of R) as a function of R, for models with 1,000 
mb, and 400 mb surface pressures! For Pg = 1,G00 mb, this quantity 
mean R is found to assume a value of 0.108, 0.189, .0.270, 0.353, 0.438, 
0.613, 0.797, and 0.993 for R = 0.0, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, and 
1.0, respectively. Thus, the range of mean R is smaller than that of 
R even though values of mean R for small values of R are greater 
than the corresponding values of R. For ‘Pg = 400, mean R assumes a 
value of -0.235, -0.077, 0.071, 0.210, 0.340, 0.582, 0.803, and 1.007 
for the values of R given by R = 0.0, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 
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and 1.0, respectively. Thus, the current total-ozone estimation proce- 
dure tends to expand (shrink) the range of mean R compared to that of 
R if the models being analyzed by it have low (high) surface pressures. 

From the results presented in the preceding paragraphs, it is 
clear that the negative values of the effective albedo is an artifact of 
the total-ozone estimation procedure. The analysis of the actual BTJV 
measurements has provided some negative values of R in the range -0.05 
to 0.0. Absence of very large negative values (~ -0.1 or smaller) in 
the results of the analysis of actual measurements, can be partly attri- 
buted to the lack of perfectly absorbing (R = 0.0) natural surfaces. 
Several values of the effective albedo measured over the polar regions 
by the BUV configuration, are found to fall in the range 1.1 -1.6. It 
is possible that some of these large values of R are due to a specular 
reflection, or some other physical reason, or due to a combination of 
circumstances. However, one more aspect of the estimation procedure has 
to be considered before assigning these large values of R to natural 
causes. ' This aspect is related to the possible errors in estimation" of 
the parameter 0 O at the place and time of observation due to some 
biased errors as a’result of the nonparallelness of the axis of the 
satellite and the optical axis of the cone of observation, or due to 
some random errors resulting from changes in the orbital parameters. 

From the plots of 1^(0.3800, 8= 0.0, 1,000 "mb, R^, 9q) vs. 0 o shown in 
Fig. 7 for four different values of R^ (= R) , we can see that an error 
of about 1° in 0g for large (~ 84° or greater) values of 0q can re- 
sult in errors of ±0.2 in the estimation of R. 
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Intensity along the nadir direction (arbitrary units) 


1 



Eig., 7 Variations of the intensity of the scattered radiation 
emerging along the local nadir direction [l c (X, 

- R, , e 0 )] as a function .of the solar zenith angle (0 q) 

A __ 

for X = 0.3800 pm, and P 0 = 1,000 mb (atmospheric, 
model Nos. 2 through 11). Different curves are for 
different values of the reflectivity (R^) of the sur- 
face' underlying the atmosphere. 
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5.3 Results foT the Cases with R Dependent up to A; In §5.1, we 

I I C 

justified the constraint of setting £2 to zero if £2 -£2 S 

x,y x,y x,y 

0.030 atm-cm in the SITC3 (or C4) procedure on the ground of a possible 
large difference between values of Rq 3800 y an£ * ^0 3398 y* or on tiie 
ground of a large sensitivity of £2^ to small changes in albedo. 

This constraint is on a very sound basis if we have an a pviovi infor- 
mation that R is independent of X. If we do have information on 
A 

variations of R, with A for various surfaces, we can modify the 
A 

total-ozone estimation procedure to take advantage of this additional 
information. The question then arises as to how much error is intro- 


duced in £2 e as a result of a given R^ vs. A variation, and 'to what 
extent the criterion of the SITC3 procedure assists in improving values 


of £2 so obtained, 
e 


Very little information is available in the ‘open literature about 
the spectral dependence of various types of natural surfaces in the 
region of immediate interest, viz, 0.3100- 0.3800 pm. We have therefore 
attempted to obtain a preliminary answer to this question by using sim- 
ulated Lambert reflectivity characteristics shown in Fig. 8. In fact, 
we have used three distinct types of characteristics. The first type 
of characteristics represents surfaces whose Lambert reflectivity in- 
creases with an increase in wavelength (curves marked A, B, and C) . 

Three different slopes of R vs. A are considered, viz, AR = 0.01, 
0.02, and 0.05 for a 0.01 ym change in wavelength. The second type 
of characteristics represents surfaces whose Lambert reflectivity 
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Simulated Lambert reflectivity 



Fig. 8 Variations of the simulated Lambert reflectivity (R^) 
as a function of the wavelength (A) . 
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decreases with an increase in wavelength (curves marked D, E, F). Again, 
the same three slopes are considered, but with an opposite sign. The 
third type of characteristics would represent surfaces whose vs. X 
curves would have either a maximum (curve G) , or a minimum (curve H) , 
in the middle of the spectral region of interest. 

The SITC2 and SITC4 procedures (see §5.1) were modified to handle 
this generalized input, and renamed SITC6 and SITC8, respectively. The 
new procedures are capable of handling cases where the Lambert reflec- 
tivity (R, ) of the surface underlying the atmospheric model, can be a 
function of wavelength (1). Atmospheric models Nos. 2 through 11, and 
13 through 22 were again used to obtain values of ft^ and R for ten 
different values of 0 q, and eight different R^ vs. X characteristics 
using the SITC6 and SITC8 procedures. For brevity, we will restrict our 

discussion to a mean value of ft where this mean is taken over a maxi- 

e 

mum of 20 subcases for a given value of the input ozone amount (ft^ ) , 

but for 10 values of 8 q and 2 values of the surface pressure (Pq ) « 

Number of subcases for which the best ozone estimate (ft ) could be eval- 

e 

uated for a given value of ft and the corresponding mean value of ft^ 
as obtained with the SITC6 and SITC8 procedures are given in Tables XII, 
XIII, and XIV for the three different types of reflectance characteris- 
tics, respectively, described in the preceding paragraph. 

From the results presented in Table XII for the models resting on 
a Lambert surface whose reflectivity increases with an increase in wave- 
length, it can be seen that the constraint of the SITC8 procedure dis- 
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Table XII. Number of subcases for which the best ozone estimate (ft ) 

e 

could be evaluated out of a maximum of 20 subcases (10 

values. of .6 0 > and 2 values of P 0 ) for. a given input ozone 

amount (^ n )> and the corresponding mean ft g for these 
subcases as obtained with the SITC6 and SITC8 procedures. 

Ground Reflectivity: R. is assumed to increase linearly 

A 

with an increase of X (curves A, B, and C of Fig. 8). . 


Number of subcases for which ft could be evaluated. 

e - - 


ft. 

Ref. 

Curve A 

Ref. Curve B 

. Ref. Curve C 

. ln 
atm-cm 

SITC6 

SITC8 

SITC6 

SITC8 

SITC6 

SITC8 


20 

17 

20 

17 

20 

15 

0.250 

20 

17 

20 

15 

. 20 

13 


20 

15 

20 

12 

19 • 

8 

0.350 

19 

10 

19 

9- 

. 19 

8 

0.400 

19 

9 

19 

8 

19 

5 

0.450 

19 

9 

19 

8 

17 

0 

0.500 

18 

8 

17 

7 

17 

0 

0.550 

17 

7 

17 

6 

15 

0 

0.600 

15 

' 7 

- 14 

5 

11 

0 " 


Corresponding value of mean ft^ in atm-cm. 



0.215- 

n 


wmm 

0.233 

0.221 

0.250 

0.270 

Wmm 

l&lfl 

mSmm 

0.292 

0.272 

■ 

0.324 

WEm 

B 

■ES 

0.337 . 

0.315 

0.350 

0.366 

0.362 

0.372 

WEM 

0.390 • 

0.367 

0.400 

0.420 

0.405 

0.428 

mmm 

0.455 

0.417 

0.450 

0.477 

0.455 

0.485 

0.458 

0.495 

- 

0.500.' 

0.527 

0.505 

0.531 

0.509 

0.554 

- 

0.550 

0.575 

0.555 

0.585 

0.559 

0.592 

_ 

0.600 

0.606 

0.606 

0.605 

0.608 

0.601 
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Table XIII. Number of subcases for which the best ozone estimate (Si^) 
could be evaluated out of a maximum of 20 subcases (10 
values of 0q, and 2 values of Pg) for a given input ozone 
amount (ft ) , and the corresponding mean fl for these 
subcases as obtained with the SITC6 and SITC8 procedures. 

G round Reflectivity: R^ is assumed to decrease linearly 

with an increase of X (curves D, E, and F of Fig. 8). 


Number of subcases for which SI could be evaluated. 

e 


S^in 

atm- cm 

Ref. Curve D 

Ref . Curve E 

Ref. Curve F 

SITC6 

SITC8 

SITC6 

SITC8 

SITC6 

SITC8 

0.200 

20 

20 

20 


• 20 

18 , 

0.250 

20 

20 

20 


20 

18 

0.300 

20 

20 

20 


20 

17 

0.350 

20 

19 

20 

18 

. 20 

13 

0.400 

20 

18 

20 

16 

18 

4 

0.450 

20 

18 

20 

14 

17 

0 

0.500 

20 

16 

19 

12 

16 

0 

0.550 

19 

14 

17 

12 

13 

0 

0.600 

16 

13 

- 13 

10 

8 

0 - 


Corresponding value of mean Sl^ in atm-cm. 


mm 

0.196 

0.196 

0.192 

nn 

0.178 

0.180 

■SS 

0.247 

0.247 

0.243 

WBESM 

0.230 

0.231 


0.297 

0.297 

0.291 

0.293. 

0.276 

0.278 

eh 

0.346 

0.347 

0.341 

0.343 

0.323 

0.324 

0.400 

0.396 

0.397 

0.389 

0.392 

0.376 

0.381 

0.450 

0. 445 

0.446 

0.436 

• 0.442 

0.421 

- 

0,500 

0.493 

0.497 

0.485 

0.491 

0.470 

- 

0.550 

0.543 

0.547 

0.537 

0.541 

0.519 

- 

0.600 

0.595 

0.597 

0.591 

0.591 

0.572 

- 
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Tabte XIV. Number of subcases for which, the best ozone estimate (fl ) 

e 

could be evaluated out of a maximum of 20 subcases (10 
values of 0 O , and 2 values of Pq) for a given input ozone 
amount (&^ n ) > and the corresponding mean for these 

subcases as obtained with the SITC6 and SITG8 procedures. 

Ground Reflectivity: Variations of R vs. X used for 

A 

these results are as shown by the curves G, and H of Fig. 8. 


Number of subcases for which 


a 


e 


could be evaluated. 


a ln 
atm- cm 

Ref. Curve G 

Ref . Curve H 

SITC6 

SITC8 

SITC6 

SITC8 

0.200 

19 

15 

20 

15 

0.250 

18 

7 

20 

14 

0.300 

18 

5 

20 

9 

0.350 

17 

0 

20 

6 

0.400 

15 

0 

20" . 

4 

0.450 

9 

0 

19 

0 

0.500 

5 

0 

18 

0 

0.550 

0 

0 

18 

0 

0.600 

0 

0 

16 ' 

0 


Corresponding value of mean in atm-c’m. 


0.200 

0.241 

wgm 

SSE9 

0.173 

0.250 

0.292 


BB 

0.226 

0.300 

0.342 

■n 

US . IS 

0.275 

0.350 

0.394 

1 

0.315 

0.329 

0.400 

0.448 

1 

0.369 

0.376 

0.450 

0.499 


0.410 

- 

0.500 

0.537 

aBMB 

0.458 

- 

0.550 

- 

I 

0.508 

- 

0.600 

- 

1 

0.556 

— 
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cards a large number of subcases at large values of the parameter ft ^ . 
This is especially. true when the slope of the R vs. X curve is also 

A 

the steepest of all ■ (curve C). Use of the SITC6 (or SITC8) procedure 
for atmospheric models resting on a Lambert surface with its reflectiv- 
ity increasing with an increase- in wavelength, results in an overestima- 
tion of mean ft (in fact, most of the individual values of ft were 
e e 

found to be greater than the corresponding value of ft ). In general, 

values of the mean ft obtained with the SITC8 procedure are signifi- 

6 

cantly closer to the corresponding ft^ values than those obtained with 
the SITC6 procedure. 

Results presented in Table XIII for the models resting on a Lambert 
surface whose reflectivity decreases with an increase in wavelength, 
also support the observations for the models with the first surface 
characteristics reported in the preceding paragraph except that the 
quantity [(mean ft & ) “ ^ ] now carries a negative sign. Furthermore, 
the magnitude of (mean ft g ) - ft^ is also generally smaller for this 
latter characteristic than the preceding one. 

When the spectral reflectance curves of the. Lambert surface under- 
lying an atmospheric model exhibit a strong maximum or a minimum in the 
middle of the spectral region of interest, 'use of the SITC8 procedure 
results in discarding of a large number of subcases (Table XIV) . Fur- 
thermore, the quantity [(mean ft ) - ft . ] is positive- (negative) when the 
reflectance characteristics of the underlying surface are represented by 
the curve G (H) of Fig. 8. It may also be noted that the values of 
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~| (mean -h [ for the results presented in Table XIV are generally 
greater than those in Table XII, or XIII. 

Other additional aspects which should be looked into in this type 
of analysis are the maximum, mean, and the R.M.S. deviations of 
around its mean value. From a cursory analysis of the detailed output, 
we feel that the use of the SITC8 procedure would also result in a sig- 
nificant improvement in this area. 

5.4 Effect of Raising, the Ground: Results presented in §5.1, §5.2, and 

§51 3 are restricted to dust-free atmospheric models with a surface pres- 
sure (Pq) of 1,000 or 400 mb. In fact, these models are the same as 
those used for the preparation of basic tables for the total-ozone esti- 
mation procedure (see §2.2). In order to investigate the effect of 
raising the lower boundary of the atmospheric model above the mean sea- 
level, the atmospheric model No. 3 and 7 with = 0.250 and 0.450 

atm-cm were used to generate several additional models whose lower 
boundaries are located at h km above the mean sea- level. The following 
values are used for the parameter h: 0, 2, 4, 6, 7, -9,. 11, 13, and 15 

km. It may be noted that the model with its lower boundary located at 
7 -km. above the mean sea- level has a surface pressure (Pq) of 400 mb. 
Information about the surface pressure (Pq in mb), and the. total ozone 
content (fi in atm-cm) of these additional atmospheric models along 
with their corresponding serial numbers as used in generation of the 
output, is given in Table XV. 
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Table XV. Values of the height (h in km) of the lower boundary of the 
atmospheric model above the mean sea- level, surface pressure 
(P 0 in mb) at this lower boundary, and the actual total 
ozone amount [S2 in atm-cm] for several additional models. 
Information about models No. 3, 7, 14, and 18 is reproduced 
here from Table III for assuring completeness of presentation 


Model Number 

Height 

(km) 

Surface Pressure 
(mb) 

Ozone Amount 
(atm-cm) 

3 

0 

1,000 

0.250 

31 

2 

789 

0.245 

32 

4 

605 

0.240 

33 

6 

461 

0.236 

14 ' 

"7 " 

' 400 

0.234 

34 

9- 

304 

0.232 

35 

'll 

228 

0.231 

36 

13 

169 

0.229 

' 37 

' 15 

‘ 123 

0.227 

7 

0 

1,000 

0.450 

38 

2 

789 

0.443 

' 39 

4 - 

605 

0.436 

40 

6 

461 

0.430 

' 18 

7 

400 

0.427 

41 

9 

304 

0.419 

42 

11 

228 

0.402 

43 

13 

169 

0.372 

44 

15 

123 

0.342 
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Values of the deviation A [= 1,000 x (ft^-ft j] are given in Table 
XVI for all atmospheric models with = 0.250 atm-cra, and for all ten 

values of 0 q, but for the wavelength- independent reflectivity of 0.1, 
and 0.2 only. (The SITC4 procedure was used for this purpose.) Similar 
results but for R - 0.8 and 1.0 are given in Table XVII. For R = 0.0 
and 0.1, the deviation A increases with an increase in the height of 
the lower boundary of the model above the mean sea-level. Large values 
of A are generally seen at small values of 6q, and at high values of 
h. For 8 0 - 0°, 45°, and 60°, it is not possible to estimate a value 

of ft for values of h > 9 km when R = 0.1. For R = 0.8 and 1.0, 

e 

the deviation A (which is positive for h = 0 km) decreases with an in- 
crease in the value of the parameter h, passes through a zero value at 
h = 0 km, and assumes negative values thereafter. 

From the values of the quantity A given in Table XVIII for all 
atmospheric models with ft^ = 0.450 atm- cm, and with low values (R - 
0.1 and 0.2) for the Lambert reflectivity of the underlying surface, we 
find that it is not possible to obtain a value of in many subcases. 

This is not the case for large values of R (0.8 and 1.0) for which 
numerical results are presented in Table XIX. Comparison of results 
presented in Tables XVII and XIX shows a very significant increase in 
the magnitude of A for given R, 8 q, h combinations as ft is in- 
creased from 0.250 to 0.450 atm-cm. 

The following analysis was performed for a part of the output for 
the atmospheric models listed in Table XV at the request of Dr. A. J. 
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Table XVI. Effect of raising, the lower' boundary by h km on the devi- 


ation [A = 1,000 x for an atmospheric model with 


ft. = 0.250 atm-cm. 

in 


e o 



Deviation for 

the model with h, = 



0 km 

2 km 

4 km 

•6 km 

7 km 

9 km 

11 km 

13 km 

15 km 




R = 

0.1; Independent of 

X. 



■ 0.0 

0 

9 

22 

36 

44 

- 

- 


- 

45.0 

0 

7 

18 

30 

37 

- 

- 


- 

60.0 

. 0 

5 

. 13 

23 

28 

- 

- 

- 

- 

70.0 

0 

3 

8 

14 

18 

27 

- 

- 

- 

75.6 

• 0 

2 

4 

8 

11 

17 

- 

- 

- 

79.6 

0 

1 

• 3 

5 

6 ' 

9 

13 

- 

- 

'82.5 

0 

3 

* 8 " 

12 

' 15 

' 23 



— ” ’ 

84.7 

0 

4 

8 

12 

14 

19 

24 

- 

- 

86.7 

0 

4 

-10 

15 . 

18 

22 

25 

- 

- 

90.0 

0 

17 

- 

- 

- 

- 


- 

- 

- 



R = 

0.2; Independent of 

X. 



- o.o‘ 

' -1 

' 5 

12 

19 ' 

23 

31 

-38 

43 

47' 

45.0 

-1 

4 

10 

17 

' 21 

' 28 ' 

34' 

39 

42 

60.0 

. .-1 

3 

8 

13 

17 

23 

29 

33 

36 

70.0 

0 

2 

5 

8 

11 

16 • 

21 

25 

28 

75.6 

. 0 

1 

3 

5 

6 

10 

14 

17 

19 ' 

79.6 

0 

1 

1 

3 

3 

5 

7 

10 

11 

82.5. 

0 

2 • 

5 

8 

10 

14 

19 

25 

- 

84.7 

0 

2 

5 

.8 

9 

11 

14 

18 

22 

86.7 

0 

3 

7 

10 

11 

13 

14 

14 

14 

90.0 

-1 

11 

26 


- 

- 

- 

- 

- 
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Table XVII. Effect of raising the lower boundary by h km on the devia- 


tion [A = 1,000 x (fi -£2^)'] for an atmospheric model with 
£2. = 0.250 atm-cm. 

in . 


e 0 



Deviation 

Hi 

o 

H 

the model with h = 



0 km 

2 km 

4 km 

6 

km 

7 km 

9 km 

11 km 

13 km 

15 km 

■ 



R = 

0. 

8; Independent of 

X. 



Q 

24 

16 

8 


2 

0 

-3 

-5 

-7 

-9 


20 

14 

8 


2 

0 

-3 • 

-5 

-7 

-9 

H 

16 

12 

6 


2 • 

‘ 0 

-‘-3 

’ -4 

-6 

-9 

70.0 

12 

9 

5 


1 

0 

-2 

-3 

-5 

-8 

75.6 

10 

7 

4 


1 

0 

-1 

-2 

-4 

-6 

79.6 

8 

6 

3 


1 

0 

-1 

-1 

-2 

-3 

82.5 

6 

4 

2‘ 


0" 

0 ‘ 

0 

0 ~ 

0 

' "-1 

84.7 

6 

4 

2 


1 

0 

0 

0 

0 

1 

86.7 

7 

5 

. 3 


1 

0 

-1 

-1 

-2 

-1 

90.0 

15 

10 

5 


1’ 

0 

-1 

-2 

-6 

-15 



• 

R = 

1. 

0; Independent of 

X. 



o.o- 

. 28 

- 19 

in 


2 • 

0 

-3 

-5 

-6 

-8- 

45.0 

25 

.. 17 

9 


2 

0 

-3 

-5- 

-7 

-9 

60.0 

21 

14 

8 


2 

0 

-3 

-5 

-7 

-10 

70.0 

16 

ii 

6 


2 

0 

-2 

-4 

-6 

-8 

75.6 

12 

8 

5 


1 

0 

-2 

-3 

-4 

-6 

79.6 

10 

7 

4 


1 

0 

-1 

-1 

-2 

-3 

82.5 

7 

5 

2 


0 

0 

0 

0 

0 

0 

84.7 

8 

5 

2 


0 

0 

0 

0 

1 

2 

86.7 

9 

6 

3 


1 

0 

-1 

-1 

0 

1 

90.0 

18 

12 

6 


1 

0 

-1 

-2 

-5 

-14 
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Table XVIII. Effect of raising the lower boundary by h km on the devi- 


ation [A = 1,000 x (ft -ft. )] for an atmospheric model with 

e xn 
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Table XIX. Effect of raising the lower boundary by b km on the devi- 
ation [A = 1,000 x (fl -jj. )] for an atmospheric model with 

e xn 

ft. = 0.450 atm-cm. 

xn 


Deviation for the model with h = 


00 



0 km 

2 km 

4 km 

6 km 

7 km 

9 km 

11 km 

13 km 

15 km 




R = 

* 0.8; 

Independent of X. 



0.0 

32 

22 

12 

4 

0 

-9 

-30 

-63 

-94 

45.0 ' 

29 

20 

11 

4 

0 

-9 

-29 

-62 

-93 

60.0 

24 

17 ! 

10 

1 

4 ; 

0 

-8 

-28 

-60 

-91 

70.0 

19 

14 1 

8 

3 

0 

-7 

-25 

-56 

-86 

75.6 

15 

10 

6 

2 

0 

-5 

-21 

-51 

-79 

79.6 

12 

9 

5 

' 2 

0 

-5 

-19 

-47 

-75 

82.5 

12 

8 ■ 

5 

2 

' 0 

-4 

-15 

-40 

-66 

84.7 

14 

10 

6 

2 

0 

-4 

-15 

-37 

-60 

86.7 

19 

14 

- 8 

3 

0 

-6 

-21 

-46 

-67 

90.0 

24 

16 

9 

3 

0 

-8 

-30 

-67 

— 


0.0 

37 

• 24 

45.0 

34 

23 

60.0 

29 

20 

70.0 

23 

16 

75.6 

18 | 

12 

79.6 

14 

9 

82.5 

14 1 

9 

84.7 

16 

10 

86.7 

22 

14 

90.0 

28 

19 


R = 1.0; 


13 

13 

12 

9 

7 

5 

4 

5 


4 

4 

4 

3 

3 

2 

1 

2 


8 


3 

3 


Independent of X. 


0 

-9 

-30 ! 

-63 

-94 

0 

-9 ■ 

-30 

-63 

-94 

0 

-9 

-29 

-62 

-93 

0 

-7 

-26 

-58 

-88 

0 

-6 

-22 

-53 

-81 

0 

-5 

-20 

-49 

-77 

0 

-3 

-15 

-40 

-67 

0 

-3 

-14 

-36 

-59 

0 

-5 

-19 

-44 

-64 

0 

-7 

-28 

-65 

- 
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Fleig: The total-ozone estimation procedure may provide values of 

I I 

£2 innn an ^ & for x = 1, and 2. Based on the value of 8 n and 

x, 1000 x,400 u 

the slope of the N c vs. ft^ curve in the region of immediate interest, 

we select one of the two wavelength pairs, and then obtain a value of 

ft £ with the help of the procedure discussed in §2.4, and §5.1. If 
I I 

values of ft , nriri and ft . are available for a given subcase and 
x,1000 x,400 

for the wavelength pair used in the determination of ft^, we can then 

obtain a value of ft 1 from the knowledge of the Pq parameter for a 

given model, and a linear interpolation or a linear extrapolation of 

ft'*' vs. P 0 , or vs. log P 0 , or vs. h curves. After several trials, 

x,y 

we found that a linear interpolation or extrapolation in Pq provides 
the best results. 

Values of ft ' were obtained after following the procedure out- 
lined in the preceding paragraph for the atmospheric models No. 31, 32, 
36, and 37 illuminated by the sun at 0°, 45°, 60°, 70°, and 75.6° from 
the -local zenith. For all cases for which it was possible to obtain a 

value of ft ’ , we found the deviation A' [= 1,000 x (ft f - ft. )] to be 
e e rn 

in the range ±4 units. Thus, for the model No. 32 with R = 0.0 and 
9 q - 0°, we have A = 34 units but A* = - 3 units only. For the 

model No. 37 with R = 0.2 and 0 q = 0°, we have A = 47 units but A' = 

0. For the model No. 37 with 0 q = 0°, we have' A = - 9 units and 

A 1 = 0 for R = 0.8, but A = - 8 units and A' = 4 units for R = 1.0. 

These results suggest that a knowledge of P 0 associated with the 
availability of the values of ft'*' 1nnn and ft"*" for a given subcase. 
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helps considerably in increasing the reliability of the best ozone esti- 
mate provided is small. 

For the moderate values of fi. (i.e., 9.. - 0.450 atm-cm), it was 

„ . . . in in . . 

not possible to obtain a value of £2^' in many instances due to the un- 

I I 

availability of both 9 and £2 For subcases for which the 

x,1000 x, 400 

quantity 9^ could be determined, the corresponding A* was found to 
be close to zero only if Pg for that particular subcase was located 
between 1,000 and 400 mb. The results were of rather poor quality when- 
- ever a linear extrapolation was required. As for example, for the model 
No. 41 and 42 with 0 q =0°, and R = 0.4, the quantity A has a value 
of - 12 and -28 units, respectively, while the respective values of the 
quantity A 1 are - 5 and - 21 units. 

Results presented in the last two paragraphs demonstrate the need 
of developing a new total-ozone estimation procedure with basic sets of 
tables for additional values of the parameter Pg if the value of the 
surface pressure at the time of observation is available from independent 
sources. It may suffice to have three additional basic sets, viz, for 
Pg = 800, 600, and 200 mb. 

5.5 Convergence of the Best-0 zone-Estimate Value: The total-ozone esti- 

mation procedure described in §2.4 can be briefly outlined as follows: 

(a) Measured values of the intensity of the scattered radiation emerging 

at the top of the atmosphere along the local nadir direction 

[I (X, 0 O )] are available for 5 different values of X (viz, 0.3125, 
m 
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0.3175, 0.3312, 0.3398, and 0.3800 pm), and a given, value of 9 q, 


(b) Use 1,000 mb tables for the given 0$ to compute Rq ^390 1000 

a 

from I (0.3800, 8n). Assign the computed value of R- nr , nrt 

— — . , d , v - . . . . U- 3oUU, luuvJ 

Q 

to -^qqq for the remaining four wavelengths. 

c c 

(c) Compute 0^ L00Q using values of Ih. 100Q and 1^(1, 9 0 > at the" 

first four wavelengths listed in (a) above. . .... 


C I 

(d) Use SI lnnri to obtain 0<jri0 , nn „ from the value of I (0.3398, 

x, 1000 0.3398,1000 m 

6 0 ), and assign Rj >33g8il000 value to R^ 1000 for the first 
three wavelengths. 

(e) Compute 0^ >1000 using values of 1Q0Q and 1^(1, e 0 ), 

OCX 

(f) Repeat (b) to (e) above to compute R^ R^ and 

o 1 

x,400‘ 

(g) Evaluate using values of ^ jl000 . ^ >400 > Q x, 1000’ and fi x,400 

after selecting the proper wavelength pair, i.e», a value for the 
subscript x. 


These steps of the procedure suggest a possibility’ of obtaining a 

converged value of by setting up an iteration scheme. The value of 

12 12 

at step (g) can be taken as n x "y to evaluate following step 

(d) . Steps- (e) , (f) , and (g) can then be executed to arrive at the 

2 

second value of the best ozone estimate represented by the symbol 0 . 


This iteration procedure can then be repeated to obtain values of 
e e 
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Values of ft^ for n = 1, 2, 3, 4, 5, and 6 were computed for the 
atmospheric model No. 2 through 11, and 13 through 22 by using the iter- 
ation procedure outlined in the preceding paragraph at the suggestion 
and request of Dr. A. J. Fleig. These computations were carried out for 
all 10 values of 0g, and 8 values of R listed earlier. A modified 
version of the SITC4 procedure was used for this purpose. 


Out of a total 800 subcases for the 1,000 mb models (10 models * 10 

values of 0 O x 8 values of R) , it was not possible to evaluate even 

1 1 n 

, o (or ft ) for 110 subcases (see §5.1). Values of ft were found to 
i , “e e - e 

.remain unchanged with an increase in the order of iteration for 615 sub- 


cases. For 59 subcases, the quantity 


ft 6 - ft. 
e an 


for a given R, 0 O , ft in 


combination was found to be greater than the corresponding value of 


ft Q 

e xn 


by 0.001 or 0.002 atm-cm suggesting a deterioration in the 


value of ft with an increase 'in the order of iteration. For 10 other 
e 

subcases, this deterioration was found to be as much as 0.007 atm-cm. 


For two subcases with ft. = 0.650 atm-cm, we could estimate ft , but 

xn e 

2 - - • 

not ft^. Only for four subcases out of a total of 800, the iteration 
procedure assisted in improving the final value of the best-osone- 


for a given R, 0 q, ft combina- 



6“ i 

< ^ 

1 

> j 

ft -ft. 
e xn 


ft - ft. 
e xn 


tion) by -0.002 atm-cm, at the most. 


Similar analysis of the results of iteration for the atmospheric 
models with 400 ‘mb surface pressure, yielded somewhat more encouraging 
information than the one reported in the preceding paragraph for the 
1,000 mb models. Out of a total of 800 subcases, it was not possible to 
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evaluate even -A for 251 subcases. For 419 subcases,, values of A n 
e e 

were found to remain unchanged with an increase in the order of itera- 
tion, n. For 18 subcases with A . = 0.650 atm-cm, we could evaluate 

m 

A , but not A . An oscillatory character of A n vs. n was observed in 
e e J e 

seven cases; e.g., for the subcase defined by A. = 0.400 atm-cm, R = 

m 

0.2, and 0 q = 75.6°, the quantity A^ was found to. change from 0.400, 

to 0.421, to 0.398, to 0.421, to 0.398, and back to 0.421 atm-cm as the 

superscript n is increased, by one, from 1 through 6, respectively. 

For 63 subcases, the final value of the best ozone estimate showed some 

6 1 

small improvement, i.e.. A' -A. - A -A. < 0.002 atm-cm. A small 

e in- e m 

deterioration in the value of A n with an increase in the order of iter- 

e 

ation was noticed in 28 subcases. The remaining 14 subcases were found 
to exhibit either a very significant deterioration (five cases) , or a 

very significant improvement (nine cases). As for example, we found 

1 2 3 4 

A = 0.614 atm-cm, A = 0.591 atm-cm, A = 0.592 atm-cm, A = 0.592 atm- 
e e e e 

5 6 

cm, A = 0.592 atm-cm, and A = 0.592 atm-cm for the subcase defined 
.e ’ e 

by A..- ?* 0.600 atm-cm, R = 0.3, .and 0n = 60°. 

J in u 

From the results presented in the preceding two paragraphs, we con- 
clude that the use "of an iterative procedure to obtain an improved value 
* t . * *. • 
of the best ozone estimate by plugging in the current ' value of the best 

ozone estimate to get a new value of the improved effective albedo 3 is 
ineffective in most, cases. For a, few cases for which some effect is- 
noticeable, the converged value may or may not represent a real improve- 
ment. 
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5.6 Recommendations: An oral presentation of some of the results dis- 

cussed in §5.1 to §5.3 was given by the author to Dr. C. L. Mateer on 
September 20-21, 1976 at the Atmospheric Environment Service, Toronto, 
Canada, and to a group of scientists (Drs. A. J. Fleig, R. S. Fraser, 
and others) on September 22- 24, 1976 at NASA/GSFC, Greenbelt, Maryland. 
A set of tentative recommendations was made by the author at these meet- 
ings, and discussed extensively. Results of further investigations 
based on these discussions are given in §5.4 and §5.5. A final set of 
recommendations which can assist in increasing the reliability and con- 
fidence level of total ozone values derived from the BUV measurements, 
is given below. Some of these recommendations include suggestions by 
various participants at these two meetings, and are not strictly con- 
fined to the primary purpose of this contract. 

Recommendation I: In §4.1, we presented results to show that the 

N^ vs. curve for the first wavelength pair (0.3312- 0.3125 pm), and 

for 0n > 79.6° has a maximum which is shifted to lower values of ft. 

- .in 

with an increase in 0 q. In §5.1, we presented several examples point- * 
ing out the large errors which can be introduced in the hest ozone esti- 
mate (ft £ ) value (in some cases) , due to an inadvertent use of the 
vs. ft. curves with a maximum. As mentioned earlier (U4 of §5.1), Dr. 
Mateer pointed out that the NASA/GSFC total-ozone procedure does contain 
an additional step which can be expected to perform the function of 
avoiding the use of such vs. ft_^ curves. However, use of an expli- 
cit criterion discarding the use of the first pair during analysis of 
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'observations for all values of 6 0 greater than a prespecified upper 
limit, seems very appropriate and it will avoid misunderstanding and 

duplication of effort on the part of future investigators. The specif i- 

... . ... ~ - 2 • - - - Z 

cation of this upper limit would depend upon the 3312 p » R o 3125 P 

(Z = C or I) combinations which can be encountered, and on the range of 
ft for a given set of basic tables. From the results presented in 
§ 5 . 1 , we can set this upper limit around 80°. 

Recommendation II: The process of digitization of the signal 

representing a 'given 1 ^( 1 , 0 o) measurement is known to introduce some 

■error. This information should be used to determine the R.M.S. error in 

a given measurement of the quantity N m [see Eq. (5)]. If the error of 

2 

measurement in N is AN units, and if the quantity ft is to be 
m . m . . . x,y 

Z 

determined within Aft atm-cm, we recommend discarding of a value of ft 

x,y 

when the corresponding slope of the N c vs. ft^ curve in the region of 

z 

interest [the quantity given by Eq. ( 6 )] is less than AN^/Aft 

units' per atm-cm. 

Recommendation III: This recommendation is for the .inclusion of 

the criterion for discarding a value of ft"*" when the corresponding 

x,y 

I C 

ft -ft is greater than, or equal to, some preset limit. This 
x,y x,y 

limit was set, somewhat arbitrarily, at 0.030 atm-cm for our investiga- 
tions reported under the heading of the SITC3 procedure (see H5 of §5.1). 
Some advantages of using this criterion are brought during the discus- 


sion of numerical results presented in Tables V through X. It is de- 

X G X 

sirable to modify this criterion to |ft -ft |/ft being greater 

x,y x,y x,y 
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than or equal to a preset limit. .Several tests with different values 
for the limit for this ratio are in order ., 


Recommendation .IV: This recommendation is related to • the relaxa- 

tion of the unnecessarily strict criterion for accepting a value of ft^ 
when the effective albedo (R) falls outside the range 0.2 -0.8. Further 
information on this aspect can be obtained in the sixth paragraph of 

§5.1. It is further recommended that if R £ 0.2 and if values of 

X X I 

^1 , 1000 and ^2,1000 are available, the value of fl ljl0Q0 should be ac- , 

cepted as ft^ only if ^q 00 > S 2 1000 ^ or c ^ s P art icular condition. 


A similar check must also be inserted for the selection between ft 

and ft^ , nr , as a final value for ft when R S 0.8. 

2,400 e 


1,400 


Recommendation V : In §5.1, we have shown that the analysis of simula- 
ted BUV measurements in the total absence of any information about the sur- 
face pressure at the time and place of observation, leads to a very sig- 


nificant over-estimation of the- total ozone amount if the sun is near the 


local zenith. or local horizon, and if the associated values of the sur- 
face pressure and the ground reflectivity are also near the upper or 'the' 
lower ends of their respective ranges. (Some of these cases- of over- 
estimation can be objectively declared as "unable to determine the best 
ozone estimate" by following the third recommendation.) Further errors 

in the evaluation of R and hence of ft are introduced: due to errors 

e 

in the determination of the value of 8 q at the place and time of ob- 
servation (§5.2) when the sun is near the local horizon. We therefore 
recommend that extra precautions be exercised before accepting the values 
of ft^ for values of 9 q > 85°. It should be further pointed out that 
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the second wavelength pair (viz, 0.3398- 0.3175 ym combination) also 
loses its sensitivity to changes in for large values of and 

00 . Under such circumstances, we recommend that the possibility of 
using the third wavelength pair (viz, 0.3398- 0.3312 pm; see Fig. 5) be 
also looked into. For the future satellite experiments aimed at measur- 
ing total ozone at middle and high latitudes with the sun near the local 
horizon, we suggest an additional measurement of intensity in a narrow 
spectral interval between 0.3175 ym and 0.3312 yin. 

Recommendation VI: In the preceding paragraph, we listed several 

.conditions under which the total ozone in an atmospheric column can be 
overestimated due to a lack of any information about the surface (or 
cloud top) pressure at the place and time of observation. The degree of 
this over estimation can be reduced by a very significant amount by fol- 
lowing the step outlined in §5.4 provided the surface pressure lies in 
between the two values of Pq for which basic tables of the total-ozone 
estimation procedure are available. (If necessary, basic tables for 
additional values of the parameter Pq can be prepared without much 
difficulty.) In principle, it is possible to obtain some indication of 
the surface pressure for each set of the BUV measurements from the analy- 
sis of the infrared measurements taken by another spectrometer on the 
NIMBUS-IV satellite. We therefore recommend an extensive effort to ob- 
tain such information about the surface pressure, and appropriate modifi- 
cations in the total-ozone estimation procedure for taking advantage of 
this additional piece of information. It should be pointed out that we 
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do not require a very accurate measurement of the surface (or cloud top) 
pressure at the place and time of observation; an error of ±50 mb, or 
even ± 100 mb, is acceptable. 

Recommendation VII: Another problem area is the lack of any infor- 

mation on the spectral dependence of reflectivity of the underlying sur- 
face in the 0.3100- 0.3400 pm region. The type of the errors that can 

be introduced in a value of the best ozone estimate (ft ) due to a lack 

e 

of our knowledge on R vs. X, is brought out in §5.3 by carrying out a 

.A 

complete analysis for eight different simulated variations of the Lambert 
reflectivity with wavelength. An extensive project to obtain R^ vs. X 
characteristics of various types of natural surfaces using aircraft 
measurements, satellite measurements, laboratory measurements, and numer- 
ical simulation is in order. 

Very recently. Dr. Fraser provided the author with a copy of the 
rough draft of a manuscript entitled, "The Apparent Spectral Ultraviolet 
Reflectances of Various Natural Surfaces" by P. M. Furukawa and D. F. 
Heath. The observations for this apparent spectral reflectance investi- 
gation were taken with single monochromators mounted on NASA’ s CV-990 
aircraft. We find from this report that the apparent spectral reflec- 
tivity of a cloud layer is only slightly dependent on wavelength in the 
spectral region 0.3312-0.3800 pm, but decreases very rapidly as the 
wavelength is decreased from 0.3312 pm to 0.3100 pm. This sharp discon- 
tinuity in the apparent reflectance vs. X curve for a cloud layer is 
attributed by the authors, to some absorption by liquid water in that 
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part of the spectrum. Further information on this aspect is very essen- 
tial. In the absence of such information, it is desirable to extend the 
total-ozone estimation procedure for simulated (and then actual) measure- 
ments in the following manner: 


Let us assume that a value of is obtained from the value of 

e 

^2 1000 anc ^ ^2 400 using the SITC4 (or its equivalent) procedure. 

For this purpose, we have made use of measurements of 1^(1, at ^ - 
0.3175, 0.3398, and 0.3800 pm. This value of can then be used to 

obtain values of the second improved effective albedo represented by 


A value of R' 


■R. from those of I (X, 0n) at X = 0.3125, 0.3312, and 0.3398 pm. 
X,y m ’ u 

,12 


0.3175,y 


can then be obtained after a linear interpolation. 


12 12 

Values of R_ and R- can then be used to obtain values of 

0.3175,y 0.3398,y 

12 

&2 y> and finally, a new value of If the original £2^ is obtained 

from measurements at the first wavelength pair,’ the procedure would re- 

12 

quire an extrapolation for obtaining Rq ^125 y ^ rom t * ie values of 

R 12 and R 12 

K 0.3175,y 0.3398,y* 

The extension of the SITC4 procedure in a manner outlined in the 
preceding- paragraph must be tested extensively with simulated observa- 
tions-. Its adoption to the analysis of the actual observations should 
await results of this extensive testing. 


Recommendation VIII: It is recommended that selected BUV data for 

which approximate coincidences of the ground-truth measurements of total 
ozone using the Dobson Spectrometer (only the so-called "direct sun mea- 
surements using A-D pairs") are available, be thoroughly analyzed after 
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including as many of the' above-mentioned recommendations, as possible. 
Effect of each of these recommendations can then be tested on actual 
■data. Such a study is bound to assist considerably in increasing the 
reliability and confidence level of the BUV, total-ozone data provided 
the sample is sufficiently large and representative. It can also draw 
our attention to other problem areas. To this effect, it is desirable 
to see if any significant differences between the ground-truth and the 
BUV ozone values are related to changes in the underlying terrain due to 
the motion of the satellite. For this purpose, we can make use of a 

c 

spread amongst four values of R^ ^gQQ y obtained from each of the four 
measurements of I m (0. 3800, 0g) taken through the duration of a given 
set of observations [i.e. , during measurements of I m (X, at ^ - 
0.3125, 0.3175, 0.3312, and 0.3398 pm which are taken at an interval of 
about 2 seconds] . 

VI. SIMULATION FOR THE SBUV/TOMS 'CONFIGURATION 

6.1 Procedure used for Analysis of the SBUV/TOMS Measurements: The main 

difference between the section of the BUV experiment on the NIMBUS-IV 
satellite, and of the SBUV/TOMS experiment on the NIMBUS-G satellite 
aimed at the estimation of total ozone in an atmospheric column, is in 
the number, of wavelengths for which measurements are available for 1 this- 
purpose. For the BUV configuration, measurements are available- -at the 
wavelength (X) given by X = 0.3125, 0.3175, 0.3312, 0.3398, and 0.3800 
pm. For the other configuration, it is planned to have an additional 
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measurement, at >0. 3600 ym. There are some other important differences 
(e.g., observations along the off-nadir directions), but our study does 
not cover them. 


This additional measurement is made in the spectral region with a 
very small absorption by ozone (Table IV) , and hence it can be used for 

c 

evaluation of ' R„ by making use of the Eq. (2) . In this section, 

0.3600,y 

we will outline the procedure used by us for the analysis of the simu- 
lated SBUV/TOMS measurements along the local nadir, for the estimation 
of, the quantities (best ozone estimate) 3 and (effective albedo ). 
.This procedure to be referred to as the SITD1 procedure, is essentially 
a modification of the SITC4 procedure discussed in the sixth paragraph 
of §5.1. In turn, this SITC4 procedure was obtained after inserting 
three important modifications in the SITC1 procedure outlined in §2.4. 

We will therefore ask the reader to refer to the discussion of §2.4 when- 
ever appropriate. 


Step I: Use the measured values- of 1^(1, 0 o) at ^ = 0.3600 and 


0.3800 ym for. .evaluating the coarse effective albedoes. Rq ^qq iq 00 

C 

and R_ oonn respectively by making use of the 1,000 mb tables for 

U • JuUU j IUUU — 

the appropriate 0g. The tables for = 0.400 atm-cm are arbitrarily 

c 

chosen for the- evaluation of the quantity Rq ^qq £000 as 310 informa- 
tion on the total ozone amount is available at this point. The ration- 
ale’ behind 1 the use of the 1,000 mb tables and for referring the quanti- 
ties computed with the help of Eq. (2) as the coarse effective albedoes 
can be found in the, discussion under the heading "Step I" in §2.4. . 
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Step II: Obtain values of R^ 100Q at A = 0.3125, 0.3175, 0.3312, 

C 

and 0.3398 pm after a linear extrapolation of the R innn vs. A curve 

A y XUUU 

obtained from values of R 0 . 3800, 1000 and R 0. 3600, 1000* 

_ c 

Step III: Follow Step III of §2.4 to obtain values of 1000 an< * 

Q 

£qqq for x = 1 only if 6 q < 79.6°, and for x = 2 for all values 
of 6 0 . 


Step IV: This Step is exactly identical to Step IV in §2.4. 

Step V: Follow Step V of §2.4 for obtaining values of the improved 
effective albedo (R £ iqoo^ at ^ ~ °’3398, and 0.3600 pm. Then, obtain 
values of R £ ^qoo at ^ = ^.3125, 0.3175, and 0.3312 ym using a 
quadratic extrapolation based on the values of R d j^qqq at 1 = 0.3398, 

0.3600, and 0.3800 pm. (Note that R 0 . 3800, 1000 = R 0. 3800, 1000^ * 

Step VI: Repeat Step III above using values of R £ ^QOO ^ nstea< ^ 

Q 

of those for ^QOO* then have values of the quantities 

I I 

Q i nnn and S for x = 1, and 2. 

x,1000 x,1000 

Step VII: Repeat Steps I through VI above by using the 400 mb 

tables instead of the 1,000 mb tables. We will then have values of 

c I ' c 

R A 400 and ®a 400 at add s£x wavelengths, and of the quantities 

and S x,400’ aS wel1 aS fi x,400 and S x,400 £or x = X » and 2 * 

Step VIII: Determine values of the effective albedo R by making 

A 

use of the following equation: 


, R 1 ' 

,1000 T K A ,400 ' 
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Step IX: Set ft to zero if 

^ x,y 


> 0.030 atm-cm for a 
implies that that par- 


T c 

| ft - ft , 
x,y x,y ( 

given x,y combination. A zero value of ft 

ticular value of ft F is unavailable. At that time, the corresponding 

x,y 

value of the slope (S^ ) is set to - 100. 

x,y / 


Step X: This final step .is for obtaining a value for the quantity 

ft , the best ozone estimate, from the computed values of ft^ for 
e x,y 

x = 1, and 2, and for y = 1,000, and 400 mb. The following paths are 

X 

taken depending upon the number of ft ' s carrying values greater than 

x,y 


zero: 


(i) All four values of ft 1 are available. Select the first 

x,y 

wavelength pair only if S 1}1000 > S 2>1000 and S lj400 > S 2j4Q0 
If not, select the second wavelength pair. 

I I 

(ii) If only ft^ 100Q and 4Q q are available, then select the 
first wavelength pair. 

I I 

(iii) If only ft 2 - 100Q and ft 2 40Q are available, then select the 

second wavelength pair. • - * - 

(iv) If the conditions listed under (i) , (ii) , or (iii) above are 

not satisfied, proceed as follows after assigning a zero value 

i 

. to ft : 
e 

i 

F ° r 5 0. 3398 s °- 2 - 


ft = Greater of 
e 


I I 

n i,1000’ ^2,1000 


( 13 ) 


For R 0. 3398 " °* 8 ’ 
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ft = Greater of 
e 


ft 1 ft 1 

w l,400* 2,400 


( 14 ) 


The zero value assigned to ft^ ’ remains .unchanged for 0.2 < 

^0.3398 < °* 8 ' ” ~ 

(v) If one of the conditions listed under (i) , (ii) , or (iii) above 
is satisfied, determine the mean effective albedo R-^ ^ for 
the wavelength pair selected as follows: 


R X ±j Aj = °- 5 


R X ± + R Xj 


(15) 


X. = 0.3312 urn, and X. = 0.3125 pm if the first wavelength is 
i 3 

selected; and X^ = 0.3398 pm and A_. = 0.3175 pm if the sec- 
ond wavelength pair is selected. Then, obtain a value for ft & 
as follows: 


®e ^x,1000 


for Rx ±, x i " 0,2 » 


(16) 


^e ^x, 400 


for Ri . \ ■ - °* 8 > 
a i, a 3 


(17) 


and 


for 0.2 < R^ < 0.8, 


ft = 
e 


°- 8-R Xi Ai i S Xi X-j ~°' 2 i . 

A>- -J- o'*" -i l J O 

0.6 5 x, 1000 0.6 x,400 


(18) 


6.2 Effective Albedo: From the outline of the procedure used for the 

estimation of total ozone with the SBUV/TOMS configuration given in the 
preceding section, we find that any real- or pseudo- dependence of the 
effective albedo on wavelength will be brought forward by this analysis. 
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In §5.2, we presented results to show that the BUV total-ozone estima- 
tion procedure can introduce an artificial dependence of R on 0 q and 

ft. , and that values of R so computed can fall outside the normal 
xn 

range (0.0- 1.0). In this section, we will present several typical ex- 
amples to show that a pseudo dependence of R^ vs. X can be exhibited 
in the results of analysis made with the total-ozone estimation proce- 
dure of §6.1, even when the simulated measurements- are for an atmospheric 
model resting on a Lambert surface whose reflectivity R is assumed to 
be independent of wavelength. 


Let us consider the atmospheric model No. 2 (see Table III) resting 

on a Lambert surface with a spectrally independent .reflectivity (R) of 

0.1. This model has a surface pressure (Po) of 1,000 mb, and a total 

ozone content (£2. ) of 0.200 atm-cm. We will further restrict our dis- 
xn 

cussion to the case when the sun is at 70° from the local zenith, i.e., 
9 0 = 70°. For this particular subcase, we find that the values of the 
quantity R? . nnn for X — 0.3125, 0.3175, 0.3312, 0.3398, 0.3600, and 
0.3800 pm to be given by 0.102, 0.102, 0.102, 0.101, 0.101, and 0.100, 
respectively. A weak spectral dependence of the coarse albedo on X in 
this particular case is due to the use of tables for 0.400 atm-cm total 


ozone (Step I of §6.1). This weak spectral dependence of the coarse 

albedo on X is not reflected in the values of the -improved effective 

C 

albedos which are calculated using the value of 1000“ We- f ind that 
R* innn = 0.100 for six values of X mentioned above. On the other 
hand, use of the 400 mb tables leads to the following sets of values for 
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R A,400 5 and R 


? , nn : R? /on = 0.334, 0.326, 0.306, 0.293, 0.263, 0.233 

A, 400 A, 400. 


and R x 40Q = 0.356, 0.345, 0.315, 0.298, 0.263,*, . and 0.233 for A = 0.3125, 
0.3175, 0.3312, 0.3398, 0.3600, and 0.3800 pm, respectively. This pseudo 
spectral dependence of R^ on A is due "to the use of- tables for a 

surface pressure different from that of the model under investigation. 
From the results presented above, we see that the effective albedo given 
by Eq. (12) will also exhibit an artificial spectral dependence in this 
particular case. (For the six values of A in the increasing order of 
A, R^ = 0.228, 0.222, 0.208, 0.199,-0.181, and 0.167, respectively.) 

For the atmospheric model No. 2 resting on a -Lambert surface with 
a spectrally independent* R of 0.8 and. illuminated by the sun at 70° 

r 

* 

from the local zenith, the above-mentioned pseudo-spectral dependence 
of R on A is somewhat smaller in magnitude than the R = 0.1 sub- 

A ■ - 

case- discussed in the .preceding paragraph. " For this latter case, values 
of R at A = 0.3125, 0.3175, 0.3312, 0.3398, 0.3600, and 0.3800 pm 

A 

are given by 0.798, 6.797, 0.795, 0.793, 0.791j and 0.789, respectively. 


We will next consider the atmospheric model No. 13 with a spectral- 
ly independent value of R .of 0.1, and 0 q = 70°. This model has a 
surface pressure of 400 mb, but a total ozone content (fL^) of 0.200 

atm-cm when the ozone contents of the deleted layers are accounted for 

* 

(Table III). The values of R for this particular subcase which is 

A 

similar to the subcase, discussed in the pen-ultimate paragraph except 
for the difference in the surface pressure, are found to increase with 
an increase in A.- R. = - 0.195, -0.173, - 0.121, - 0.092, - 0.038,. and 

A 
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-0.002 for X = 0.3125, 0:3175, 0.3312, 0.3398, 0.3600, and 0.3800 pm, 
respectively. 

6.3 Results for the R-Independent-of-\ Cases: In 1 this section, we pro- 

pose to compare values of the best ozone estimate (52 ) as obtained from 

6 

the SITC4 procedure described in §5.1 (analysis of the simulated BUV 
measurements at five wavelengths) , and from the SITDl procedure described 
in §6.1 (analysis of the simulated SBUV/TOMS measurements at six wave- 
lengths). This comparison will be carried out for the atmospheric models 
No. 2 through 11, and 13 through 22 (Table III) resting on a Lambert sur- 
face whose reflectivity R, is assumed to be independent of wavelength. 

A 

As mentioned earlier, the quantity 52 was evaluated for 10 different ' 
directions of illumination of an atmospheric model from above, and for 8 
different values of the Lambert reflectivity of the surface underneath 
the model. 

Values of the deviation [A = 1,000 * (52 -52. )] for the atmospheric 

e in 

models. No. 6 and 17 with an average amount of total ozone (52. = 0.400 

xn 

atm-cm) are given in Table XX and XXI, respectively, for. 10 different 
values of the parameter 9 0 ,. and 8 different values of the reflectivity 
(R. independent of X) . Values of A as obtained with the different 

A 

procedures (viz, SITC4 and SITDl) are given in adjacent columns. 

For the atmospheric model with 1,000 mb surface pressure (Table XX), 
values of the deviation A as obtained with the SITC4 and SITDl are 
practically the same for low (viz, R = 0.0, and 0.1), and high (viz. 
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Table XX. Values of the deviation [a = 1,.000 x (ft -ft. )] for the atmo- 

e in 

spheric model No. 6 as obtained using the SITC4 and SITDI 
procedures, ft = 0.400 atm-cm; P 0 = 1 9 000 mb. 

Ground Reflectivity: (= R) independent of X. 
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Table XXI. Values of the deviation '[ A = 1,000 x (ft --ft. )] for the atmo- 

e m 

spheric model No. 17 as obtained using the SITC4 and SITDl 
procedures, ft. = 0.400 atm-cm; P 0 = 400 mb. 
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R = 0.8, and 1.0) values of the parameter R. • For the intermediate 

values of the Lambert reflectivity of the surface underlying the model, 

the deviation for a given 0 q, R combination as obtained with the 

SITDl procedure, is significantly greater than that obtained with the 

SITC4 procedure. It may be noted that it is possible to evaluate 

.for four subcases (0 q = 90°, R = 0.1, 0.2, 0.3, and 0.4) with the SITDl 

procedure, but not with SITC4. The mean values of where the means 

are taken over all 8 values of R, and over all 10 values of 0 q for 

which an Q is available, is 0.406 atm- cm and 0.409 atm-cm for the 
e 

SITC4 and SITDl procedures, respectively. Similar analysis of re- 
sults for the models with 1,000 mb surface pressure suggests that the 
SITDl procedure tends to over-estimate the total ozone of an atmospheric 
column by a significant amount when its results are compared with those 
obtained with the SITC4 procedure. 

For the atmospheric model with 400 mb surface pressure (Table XXI), 
we find that the use’ of the SITDl procedure results in the recovery of 
meaningful values for the quantity for several subcases which were 

previously declared "undeterminable" by the SITC4 procedure. For the 
models with low surface pressures, the SITDl procedure has a tendency 
of over-estimating the total ozone amounts • at low values of R, but of 
under-estimating them at intermediate values of R. 

Mean values of the best ozone estimate as obtained with the SITC4 

and SITDl procedures and where a mean for a given atmospheric model is 

taken over all available values of Q (10 values of 0n, and 8 values 

e 
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of R) , are given in Table XXII for the models No. 2 through 11, and 
13 through 22. Results presented in this table support our findings of 
the preceding two paragraphs. Mean values of £2^ as obtained with the 
SITDl procedure are significantly greater (smaller) than those obtained 
with the SITC4 procedure for the models with 1,000 (400) mb surface 
pressure. Comparison of the grand means taken, over a maximum of 160 
subcases (2 values of Pg x 10 values of 0 O x 8 values of R; results 
presented in columns 4 and 7 of Table XXII) suggests ' that, on the aver- 
age, the SITDl procedure tends to provide results of somewhat better 
quality than the other procedure especially if £2_^ is small-. How- 
ever, this statement should not be construed to imply any superiority 
of the SITDl procedure over the SITC4 procedure. A more realistic in- 
terpretation of the results presented in Tables XX, XXI, and XXII will 

be as follows: Errors in the values of the best ozone estimate (£2 ) as. 

e 

obtained with the SITC4 and SITDl procedures, are of the same magnitude 
for all practical purposes. The distribution of this error amongst 
values of £2^ for various subcases of a given model, depends signifi- - 
cantly on the 'procedure used. 

6.4 Results for the Cases with R dependent upon In §5.3, we 

present ecf - r e s' ui't's f o r the analysis of the simulated BUV measurements 
for the atmospheric models resting on a Lambert surface whose reflec- 
tivity OR.) is varied with wavelength. These results were presented 
for three distinct types of simulated reflection' characteristics. The 
first .(second) type of characteristics represents surfaces whose Lpnbert 
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Table XXII. Mean values (mean over a maximum of 80 subcases; 10 values 
of 0q» and 8 values of R) of the best ozone estimate for 
atmospheric models with a given surface pressure (Pq in 
mb), and a given amount of ozone (ft . in atm- cm) . 


fi in 



Mean ft 

e 

(atm- cm) 



PROCEDURE USED: 

S1TC4 

PROCEDURE USED: 

SITDl 

atm- cm 


Pq = 400 mb 

mean of cols. 
2 and 3 

P 0 = 1000 mb 

Pq = 400 mb 

mean of cols. 
5 and 6 

0.200 

0.203 

0.207 

0.205 

0.204 

0.202 

0.203 

0.250 

0.254 

0.256 

0.255- 

0.256 . 

0.252 

0.254 

0.300 

O'. 305 

0.304 

0.305 

0.307 

0.301 

0.304 

0.350 

0.356 

0.352 , 

0.354 

0.358 . 

0.350 

0.354 

0.400 

0.406 

0.401 

0.404 

0.409 

0.399 

0.403 

0.450 

0.458 

0.450 

0.454 

0.460 

0.448 

0.454 

0.500 

0.509 

0.499 

0.504 

0.512 

0.496 

0.504 

0.550 

0.560 

0.548 

0.554 

0.563 

0.545 

0.554 

0.600 

0.611 

0.598 

0.604 

0.613 

0.593 

0.603 

0.650 

0.650 

0.646 

0.648 

- 

0.640 

- 
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reflectivity increases (decreases) with an increase of wavelength, and 

l 

is represented by the curves A, B, and C (D, E, and F) of Fig. 8. The 
third type of characteristics represents surfaces whose R vs. X curve 

A 

has either a maximum (curve G) , or a minimum (curve H) in the middle of 
the spectral region of interest, viz, 0.3100- 0.3800 pm. These results 
were obtained using two different procedures named SITC6 and SITC8. We 
concluded that, on the average, results obtained with the SITC8 proce- 
dure are of somewhat better quality than those obtained with the SITC6 
procedure. 

In order to investigate the effect of the spectrally-dependent Lam- 
bert reflectivity on the best ozone estimate value (ft £ ) obtained with 
the SITDl procedure specially developed for the analysis of the simu- 
lated SBUV/TOMS observations (§6.1), -the program of the SITDl procedure 
was modified to handle cases for which the Lambert reflectivity of the 
surface underlying the atmospheric model can be an arbitrary function 
of wavelength. This modified version of SITDl will be called SXTD2. 
Atmospheric models No. 2 through 11, and 13 through 22 (Table III) were 
again used to obtain values of ft and R at six wavelengths for 10 

6 A 

different values of 0 q listed earlier, and 8 different R^ vs. X 
curves of Fig. 8. 

For brevity, we will again restrict our discussion to a mean value 

of ft where this mean is taken over a maximum of 20 subcases for a 
e 

given value of the input ozone amount (ft_^) , but for 10 values of 0 q, 
and 2 values of the surface pressure (Pq). It should be again pointed 
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out that the analysis of results presented in this subsection is, by no 

\ 

l 

means, complete. Other additional aspects which should be looked into, 

are the maximum, mean, and R.M. S. deviation- of ft around its mean val- 

e 

ue. The deviation [1000 x (ft g - ft ) ] will be found to depend very sig- 
nificantly on various input parameters such as ft. , 0 q, Fo» and R, vs. X. 

in a 

characteristics . 

From the results presented in Table XXIII for the models resting oh 
a Lambert surface whose reflectivity increases with an increase of' wave- 
length, it can be seen that the availability of the additional measure- 

# * 

ment at 0.3600 ym (i.e., the use of SITD2 procedure in the place of 

SITC8) results in a very significant increase- in the number of subcases 

for which a value' of ft can be evaluated. Mean values of ft obtained 

e - e 

with the SITD2 procedure are nearer to the corresponding values of ft^ n 
when compared to those obtained with the SITC8 procedure. This improve- 
ment' is in the right direction but probably not of the expected magni- 
tude. Since. R^ is a linear function of X and as the new estimation 
procedure (SITD2) simulates this trend very accurately, one would expect • 

ft. and the corresponding mean ft for the SITD2 procedure to be exact- 
in e 

ly equal. We found that the absolute difference between the individual 

values of ft and ft. not to exceed 0.003 atm-cm for all models with 
e in 

P'o = 1,000 mb, and 0 O < 86.7°. For F 0 = 1>000 mb models, the quantity 

ift -ft. I exceeded 0.003 atm-cm limit only six times, and exceeded 0.005 
'em 1 

atm-cm limit only twice. It was never observed to exceed 0.009 atm-cm 

value. Thus, most of the differences between the mean ft of SITD2 and 
’ e 
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Table XXIII. Number of subcases for which the best' ozone estimate (ft ) 

s e 

• l 

could be evaluated out of a maximum of 20 subcases (10 

values' of 8 qj and 2 values of Pq) for a given input ozone 

amount (ft. ) and the corresponding mean ft for these 
xn e 

subcases as obtained with the SITC8 and SXTD2 procedures. 

Ground Reflectivity: R^ is assumed to increase linearly 

with an increase of A (curves A, B, and C of Fig. 8).- 


Number of subcases for. which ft could be .evaluated. 

e 


^in 

atm- cm 

Ref. Curye A. 

Ref. Curve B 

Ref . Curve C 

SITC8 

SITD2 

S1TC8 

SITD2 

SITC8 

SITD2 

0.200 

« 

20 

17 

20 

15 

19 

0.250 

wBm ■■/M 

20 

15 

20 • 

13 

20 

0.300 

15 

20 . 

1.2 

20 

8 

20 

0.350 

10 

20 

9 

20 

8 

20 

0-. 400 

9 

20 

8 

20 

5 

20 

0.450 

9 

20 

8 

20 

0 

20 

0.500 

8 

19 

7 

19 

0 

19 

0.550 

7 

18 

6 

18 

0 

18 

0.600 

7 

•16 

5 

16 

0 

16 


Corresponding value of mean 

ft in atm-cm. 

„ e 



0.211 

m 

0.214 

0.205 

0.221 

0.208 

1 

0.262 

gjpp 

0.264 

0.256 

0.272 

0.260 

0.300 

0.312 

ins 

0.311 

' 0.305 

0;315 

-0.309 


0.362. 

0.353 

0.357 

0.354 

0.367 

0.358 

0.400 

0.405 

.0.402 

0.408 

, 0.405 

■ 0.417 

0.409 

. . 0.45.0 

0.455 

0.452 

0.458 

0.454 

- 

0.459 

II 

0.505 

0.501 

0.509 

. 0.503 

- 

0.508 

- 0.550 . 

0.555 

0.550 

0.559 

0. 552 

- 

0.557 

HEBSB 

0.606 

0.600 

0.608 

0.602 

- 

0.606- 
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corresponding • ft_^ in Table XXIII are due to significant differences 

\ 

between the individual values of ft and ft. for the models with 400 

. e xn 

mb surface pressure. The reason for the over-estimation of ft for- 

£ 

the 400 mb models is that the corresponding R^ values are very small 

(< 0.2), and hence the estimation procedure favors the value of ft^ j_qqq 

over that of ft'*' /nri . It should be added that a subcase which is dis— 
x,400 

carded by the SITC8 procedure but is accepted by the SITD2 procedure, 
generally carries a large value for the difference, ft & - ft . As for 
example, for the model No. 13 resting on a Lambert surface having the 
.reflectivity curve C of Fig. 8, and illuminated by the sun at local zen- 
ith, ft -ft. = 0.037 atm-cm. 
e in 

Results presented in Table XXIV for the models resting on a Lambert 
surface whose reflectivity decreases .with an increase of the wavelength, 
also show a very significant increase in the number of cases for which 
ft^ could be evaluated with the SITD2 procedure. The differences be- 
tween the mean ft e and the corresponding value of ft^ are very small 
for this second type of surfaces compared to those for the first type 
of surfaces. For the atmospheric models resting on Lambert surfaces 
with -their spectral- reflectivity decreasing with an increase of wave- 
length, we- found the total ozone content of the column to be over-esti- 
mated for the models with 1,000 mb surface pressure, and. to be under- 
estimated for those with 400 mb surface pressure. The magnitude of the 
over-estimation and under-estimation are practically equal for a given 
ft^, especially for the models resting on a surface with its R^ vs. A 
variation represented by the curve F of Fig. 8 
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Table XXIV. Number of subcases for which .the best ozone estimate (ft ) 

\ 

could be evaluated out of a maximum of 20 subcases (10 
values of 0g, and 2 values of P 0 ) for a given input ozone 
amount (ft^) » and the corresponding mean ft^ for these 
subcases as obtained with the SXTC8 and SITD2 procedures. 


Ground Reflectivity: is assumed to decrease linearly 

with an increase of X (curves D, E, and F of Fig. 8). 


Number of subcases for which ft could be evaluated. 

e 


ft. 

in 

atm-cm 

Ref . Curve D 

Ref. Curve E 

Ref. Curve F 

SITC8 

SITD2 

SITC8 

SITD2 

SITC8 

SITD2 

0.200 

20 

20 

20 

20 ' 

18 

20 

0.250 

20 

20 

20 

20 

18 

20 

0.300 

20 

20 

18 

20 

17 

20 

0.350 

19 

20 

18 

. 20 

13 

20 

0.400 

18 

20 

16 

20 

4 

20 

0.450 

18 

20 

14 

20 

0 

20 

0.500 

16 

20 

. 12 

20 . 

0 

20 

0.550 

14 

20 

12 

20 

0 

20 

0.600 

13 

19 

10 

19 

0 

19 


Corresponding value of mean ft g in atm-cm. 


0.200 

0.196 

mmm 

■SE9 

mrm 

0.180 ' 

■ 0.200 

0.250 

0.247 

8119 


■sill 

' 0.231 

0.251 

0.300. 

0.297 

mmm 

tn 


0.278 

0.300 • 

0.350 

0.347 

0.353 

0.343 

0.352 

0.324 

0.350 

0.400 

0.397 

0.403 

0.392 

0.403 

0.381 

0.401 

0.450 

0.446 

0.454 

0.442 

0.453 

- 

0.451 


0.497 

0.504 

0.491 

0.503 

- 

0.501 

0.550 

0.547 

0.555 

0.541 

0.554 

- 

0.550 


0.597 

0.604 

0.591 

0.603 

- 

0.600 
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When the R^ vs. A curves of the Lambert surface underlying an 

' \ - 

atmospheric model exhibit a strong maximum or a minimum in the middle of 

the spectral region of interest, use of the SITD2 procedure results in a 

significant increase (decrease) in the number of subcases for which 

could be evaluated provided is large (small). However, an increase 

in the number of subcases for which the S2 is recovered is also gener- 

£ 

ally associated with a very significant under— estimation of mean 52^ for 

models with surface having its R^ vs. X variation represented by the 

curve H of Fig. 8. As for example, for = 0.500 atm-cm, the value 

of the mean fi = 0. 424 atm-cm. This under-estimation was generally 
e 

found to be very large at large values of 0q. As for example, for the 
atmospheric model No. 18 resting on a Lambert surface with the curve H 
type reflectivity, we found the devt-cvtion (A) to be — 58, -44, — 28, — 63, 
- 53, - 50, - 54, - 70, and -131 units for the values of 6 0 given by 
6 0 = 45°, 60°, 70°, 75.6°, 79.6°, 82. 5° , ' 84. 7° , 86.7°, and 90°, respec- 
tively. Thus, we may conclude that the use of SITD2 procedure leads to 
better values of compared to those obtained with SITC8 procedure - 

provided the spectral reflectivity of the surface underlying the atmo- 
sphere increases or decreases linearly with the wavelength. For surfaces 
with other .types of characteristics, results are of mixed quality, in- 
deed. 

6.5 Recormendatvons : A procedure, for the estimation of total ozone in 

an atmospheric column using measurements at six wavelengths (SBUV/TOMS 
configuration) was outlined in §6.1. This procedure referred to as the 
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Table XXV. Number of subcases for which the best ozone • estimate (52 ) 

e 

could be evaluated out of a maximum of 20 subcases (10 
values of 6 q, and 2 values of Pq) for a given input. ozone 
amount (52^) > aad the corresponding mean 52^ for these sub- 
cases as obtained with the SITC8 and SITD2 procedures. 

Ground Reflectivity: Variations of R^ vs. A. used for 

these results are as shown by the curves G and H of Fig. 8. 


Number of '.subcases for which 52 • could be evaluated. 

e 


n in 

atm- cm 

Ref. Curve G 

Ref . Curve H 

• SITC8 

SITD2 

SITC8 

SITD2 

0.200 

15 

4 

15 

8 

0.250 

7 

. 4 

14 

9 

0.300 

5 

3 ’ 

.9 

9 

0.350 

0 • 

5 

6 

10 

' 0.400 

0 

5 

4 

9- 

0.450 

0 

6 

0 

10 

0.500 

0 

. 7 

■ 0 

8 

0.550 

0 

9 

0 

11 

0.600 

0 

5 

0 

8 

Corresponding values of mean 52^ in 

atm- cm. 

■ 0.200 

0.231 

0.208 

0.173 

0.181 

0.250 

0.275 

0.259 

0.226 

0.233 ' 

0.300 

. 0.325 

0.306 

0.275 

0.280 

' O'. 350 

— 

0.358 

0.329 

0.328 

0.400 

- 

0.404 

0.376 

0.375 

0.450 

- 

0.450 

- 

0.418 

0.500 

— 

0.519 

- 

0.424 

0 . 550 

— 

0-.576 

- 

0.502 

0.600 

- 

0.613 

- 

0.549 
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SITD1 procedure (SITD2 procedure if the simulated measurements are for 
the atmospheric models resting on a Lambert surface with spectrally de- 
pendent reflectivity) contains recommendations I, III, and IV (§5.6) for 
■a better analysis of the BUV measurements (i.e., the procedure SITC4 of 
§5.1) -which are made at five wavelengths only. Besides, the additional 
measurement at 0.3600 pm is used in the SITDl procedure for obtaining 

values of the coarse and -improved effective atbedoes at shorter wave— 

Z 

lengths by extrapolation of the R^ ^ vs. X curves on the longer wave- 
length region. Results obtained with the analysis of simulated measure- 
ments for the five, and six' wavelength configurations, were compared in 
§6.2, §6.3, and §6.4. The following recommendations are based on the 
aforementioned analysis. 

Recommendation I: The SITDl procedure should be modified for com- 
puting values of at X = 0.3125, 0.3175, and 0.3312 pm using a 

x ,y 

linear extrapolation procedure in the place of the 'current quadraic ex- 
trapolation procedure (Step' V of '§6.1). Values of ft obtained with 

0 

such a modified procedure should be compared with the corresponding val- 
ues of obtained with the original procedure for all types of models - 

discussed -in §6.3 and §6.4. This additional test can assist in improv- 
ing the poor quality of values for the- atmospheric models resting 

on surfaces with their R^ vs. X characteristics represented by curve 
H of F-lg. 8 (see last paragraph of §6.4). 

Recommendation XI: For the SBUV/TOMS configurations, plans are 

* 

made for obtaining the best ozone estimate values from the analysis of 
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measurements of intensity along directions making some .angle with the 
local nadir direction. Sufficient evidence is presented in the preced- - 
ing sections to demonstrate the unexpected (and probably unpredictable) 
nature of the errors in due to .absence of any information on the 

surface pressure, and on the spectral dependence of the. surface reflec- 
tivity at the place and time of observation. We therefore recommend 
that' the analysis of the type outlined in the preceding sections, should 
also be carried for several additional directions of observation. 

Re commendation III: Appropriate versions of the Recommendations 

IX, V, VI, and VII given in §5.6 should also be included for the analy- 
sis of actual SBUV/TOMS observations for the estimation of total ozone 
in an atmospheric column. 

Recommendation IV: It is. recommended that an. extensive analysis of 

selected SBUV/TOMS data for which approximate coincidences of the highly 
reliable gvound-tvuth measurements of total ozone are available', be car- 
ried out. Attempts should be made to encompass as many different situa- 
tions as possible for this selected study. The analysis of the entire 
SBUV/TOMS data must await results of such a detailed study. 


VII. CONCLUSION' 

The 1 procedure for estimation of, total ozone in an atmospheric col- 
umn’ from the measurements- of the ultraviolet radiation back-scattered 
along the local nadir direction by the earth-atmosphere system, makes 
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use of a quantity called N which is proportional to -the logarithm of 

t 

the ratio of intensities at two. wavelengths . One of these wavelengths 
is located in the spectral region of moderate absorption by ozone, while 
the other is located in the region of weak absorption. This quantity 
given by Eq. (4), is a function of the wavelength pair used, of the 
amount of ozone in the atmospheric column, zenith angle of the sun, re- 
flectivity of the underlying surface, and of. the surface pressure at the 
bottom of the atmospheric column defined by the ground or the top of the 
cloud layer. Variations of as a function of the various parameters 

listed above,, and computed for models of the terrestrial atmosphere free 
of aerosols, water droplets, and snow crystals are discussed in §4. 

The total-ozone estimation procedure currently being used at NASA/ 
GSFC for the analysis of BUV measurements was modified for the analysis 
of simulated measurements of intensity. This modified' total-ozone esti- 
mation procedure is described in §2.4. The original as well as the 
modified procedures attempt to estimate total ozone content of an atmo- 
spheric column (called the best ozone estimate) from the measurements of 
intensities at wavelength pairs for which tables 'of the quantity 
are available, from a knowledge of the solar zenith angle at the place 
and time of observation, and from the knowledge of the so-called effec- 
tive albedo based on measurements of intensities in the adjacent spectral 
regions with little or no absorption by ozone. No- information ‘is avail- 
able about the surface pressure, and actual reflectivities of the under- 
lying surface at various wavelengths of interest. Comparison of the 
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numerical results obtained from the analysis of simulated measurements 
at five wavelengths (BUV configuration) using the - original procedure 
and several modifications to it. based on the findings of §4, is made in- 
15. 1 through §5.5. This comparison then culminates into a set of recom- 
mendations (§5.6) which are too numerous to be summarized here. 

The experience gained during the investigation reported in §5 is 
then used to set up an estimation procedure for total ozone with the 
'SBUV/TOMS configuration where one additional measurement of intensity 
is also available (§6.1). Results obtained with the final version of 
the procedure for the BUV' configuration, and with the procedure for the 
SBUV/TOMS configuration are compared in §6.2,. §6.3, and §6.4. Recom- 
mendations which can assist in improving’ the quality of ozone values 
obtained with the procedure described in §6.1, are then outlined in §6.5. 

Analysis presented in this particular "report is for the so-called 
Rayleigh-, or molecular-scattering models of the terrestrial atmosphere. 
In other words, these models are assumed to be free .of any particles 
comparable to, or larger than, the wavelengths of the electromagnetic 
radiation specified for measurements. Analysis of numerical results 
for atmospheric models containing spherical polydispersions of known 
characteristics, will be the subject matter of the final report with 
the same- title as' this one. 
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APPENDIX A 


A. 1 Description of the Program: In §2.4, §5.1, §5.3, §6.1, and §6.4, 

we mentioned several programs (SITCl, • SITC2, SITC3, SITC4, SITC6, SITC8, 
SITDl, and SITD2) for estimating total ozone from the simulated measure- 
ments of the intensity of the scattered radiation emerging at the top 
of an atmospheric model along the local nadir direction. These proce- 
dures are arrived at after carrying out minor changes, (from the program- 
ming considerations) in a basic program, i.e. , SITCl. We will there- 
fore describe only one program (viz, SITC4) , in great, detail, in this 
section. The listing of the FORTRAN statements, and printout for a test 
run with simulated measurements for one atmospheric model are ..given in 
§A. 2 and A. 3, respectively. The subroutine OMESLO which is being 
called by all programs, will be described first. In Table A. I, we have 
given mathematical equivalence as appearing in the text and/or defini- 
tions of important variables appearing in the subroutine OMESLO, and 
the calling program SITC4. 

Subroutine OMESLO: The function of this subroutine is to return 

Z Z 

values of the quantities ft and S for a given value of Sn, and 

x,y x,y • u 

for the stated values of the subscripts x and y, as well as the super- 
script Z. It has 11 arguments and contains 66 statements. 

The first three arguments (viz, EIC, TIC, and SBAR) are for speci- 
fying locations for the quantities I (X, ft. , Pq , R. = 0. 0, 0g) , T (X, . 

c in A c 

ft. , P 0 , 9g) , and S(X, ft. , Pq) defined in §2.2. These are essentially 
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Table A. I Variables appearing in the subroutine OMESLO, and program 
SITC4 and their mathematical equivalence or definition. 

Variable Equivalence or Definition 

A Variable for reading in values of I c (A, Pq, R^ = 0.0, 

0, 0 o ) for a given model defined by the parameters £2 
and P 0 . These values are available on the disc (Data 
File No. 16) for 18 values of 0, for 10 values of 
0Q 9 and for 6 different wavelengths. See Hi and H2 
of §2.2 for the definition. 

ALBD vector called by the subroutine OMESLO. Z = C 

x,y 

or I (see Steps. I, II, and V of §2.4). y = 1,000 mb 
or 400 mb. 

ALDA • Wavelength X . 

B Variable for reading in values of T c (A, £2, Pq, 0, 0q) 

for a given model. See Hi and 1F2 of §2.2, and also 
the definition of the quantity A above. 

BEST Best ozone estimate, £1 ; Eqs. (9), (10), and (11). 

e 

C Variable for reading in values of S(A, £2, Pg) fo,r a 

given model. See HI and H2 of §2.2, and also the 
definition of the quantity A above. 

DEV £2 - £2 

e 

C 

EFFALC R ; coarse effective albedo at the wavelength X 

A , y 

obtained by using y mb pressure tables; see Step I 
and II of §2.4. 
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Table A.I 3 cont’d. 


Variable 

Equivalence or Definition 

EFFALI 

R* ; improved effective albedo ; see EFFALC. 

x>y 

EIG 

I (X, & , P 0 , R = 0.0, 6q) for 6 values of X, 20 

c in A 

basic atmospheric models [the first (second) 10 val- 
ues of the second subscript of EIC for 10 increasing 

values of ft. , but for Pn = 1,000 (400) mb], and 
m u 

10 values of 0g (the last subscript of EIC). 

EIM 

I m (X,0 o ); §2.4. 

XGRP 

Parameter representing the ground pressure value for 
the basic tables; IGRP = 1, and 2 for Pq = 1,000, 
and 400 mb, respectively. 

ISLT 

Parameter selecting the wavelength (or wavelength 
pair) to be used for the next set of computations. 

LAMl, LAM2 

Parameters for selecting individual wavelengths of a 


pair; see Eq. (4) and the statement No. SITC 11. 

NMOD 

Model number for which values of ft and R are to 

e 

be obtained. 

NTHETO 

Serial number of the subscript of the parameter 0q 
(THETO) ; see the statement No. SITC 10. 

OMEGA 

z 

Total ozone amount, ft for Z = C, or I; x = 1, or 

x,y 

2, and y = 1,000, or 400 mb. 

OMEGC 

r 

ft ; see Step III of §2.4. 

x,y 

OMEGI 

ft"*" ; see Step VI of §2.4. 

x,y’ 
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Table A.I, cont’d. 


Variable 

Equivalence or Definition 

OMEGIN 

n. ; see the second paragraph following Eq. (1). 
xn 

RAT 

given by Eq. (4); for a given wavelength pair, 
and given values of the parameters P 0 and 0 q. 

RATM 

N : see Eq. (5). 
m’ 

REF 

R; Lambert reflectivity of the surface underlying the 
atmospheric model. 

SBAR 

S(A, P 0 ) for 6 values of A, and 20 basic atmo- 

spheric models (see definition of EIC) . 

SLOPE 

S Z ‘ ; see Eq. (6) . 
x,y 

SLPC ■ ' • 

S C ; see Step III of §2.4. 
x,y 

SLPI 

S^ ; see Step VI of §2.4. 
x,y 

STALB 

£ R summed over all values of 9 q for a given value 
of R (quantity REF above); see Eq. (12). 

TEMA, TEMB, TEMC 

Temporary variables. 

THETO 

0 O ; solar zenith angle. 

TIC 

1^(1, Pq, 0q) for 6 .values of A, 20 basic, atmo- 

spheric models (see definition of EIC) , and 10 values 


Of 0Q. 

TITB ■ 

Title information. 
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the basic tables of these quantities for six different values of wave- 
lengths used in the SBUV/TOMS configuration (only five for the BUV con- 
figuration), ten atmospheric models with P 0 = 1,000 mb and the next ten 
models with Pq = 400 mb, and for ten different values of the solar zen- 
ith angle. The input ozone amounts for these models are those listed in 
Table III; models No. 1 and 12 are not to be included. 


The next two arguments (viz, LAMl and LAM2) specify wavelengths A^ 

and A. for which the measurements of I (A, 8 0 ) are to be analyzed- for 
1 m 

Z Z 

obtaining values of ft and S . The sixth argument (NTHET0) is for 

x»y x,y 

specifying one of the ten preselected (statement No. SITC 10) values of 

the solar zenith- angle, 0 q. The seventh argument (ALBD) carries values 

for the quantities and R^_. to be used in computations of N c 

given by Eq.- (4). The eighth argument. IGRP = 1 -if 1,000 mb tables are 

to be used, and = 2 for the use of 400 mb tables. Values of I m (A, 9 q) 

for six wavelengths, but for a given value of 8q are available in the 

Z Z 

locations specified by EIM, while - computed values of ft and S are 

x,y x,y . 

returned in locations marked OMEGA and SLOPE, respectively. 


Statements No. 0MESL 3 through 41 contain comments, a format state- 
ment, a dimension statement, and a data (ft. ) initialization statement. 

xn 

The quantity N (A., A.; 0 q) given by Eq. (5) is computed at the end of 

in l j 

successful execution of statement No. 0MESL 42. 


The DO LOOP 200 contained between the statements No. OMESL 43 
through 56 is for computations of values of the quantity N^ for all 
ten values of fi^, but for the given values of the parameters A^, -A^. , 
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R X-> r X'» 0 O» and The quantity N is set to - 999.99 (OMESL 55) 

J *- 

under an unlikely condition of finding a very small (< 10 -1 °) value for 
the computed intensity, I c (OMESL 49). An appropriate message (OMESL 53 
and 54) is also printed out at the terminal. 


■ Statements contained in the DO LOOP 300 (OMESL 57 through 62) are 

Z Z 

then executed for computing values of (OMEGA) , and S x (SLOPE) 

C 

by making use of the procedure outlined in §2.3. The quantity is 

computed by backward extrapolation if a given value of N^ is less than 

the first computed value of N which is for £2. = 0.200 atm-cm. The 

cm 

quantities OMEGA and SLOPE are set to 0.0, and - 100.0, respectively 


when a given value of N is greater than the last computed value of 

N which is for ft. = 0.650 atm-cm (OMESL 63 and 64). 
c • xn 


Program SITC4: This program contains a total of 235 statements. . 

1 

Statements No. SITC 1 through 15 contain comment cards, dimension state- 
ments, and data initialization statement. Further information about the 
DEFINE FILE 16 statement (SITC 13) can be found on page No. 55 of the 
Technical Report: I of this contract. The first 22 records of this 
data set contains values of I , T^, S for 22 atmospheric models listed 
in Table III. 


The input to this program is the values of 1^, T^, and S for the 
22 atmospheric models listed in Table III, and also for the atmospheric 
model (NM0D, see SITC 61) for which' values of the best ozone estimate 
(Q & ) and effective albedo (R) are to be computed for 10 values of 0 O 
(SITC 10), and 8 values of R (SITC 12). A given record on the data- set 
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No. 16 contains 14 words of the title information describing the atmo- 
spheric model for which values of and are available for 18 val- 
ues of 0 -and 10 values of 0 O , and of 1^, and S for 6 wavelengths. 

The function of the part of the program contained between the state- 
ments SITC 40 through 59 is to read in data for 20 atmospheric models 
with non-zero ozone content (Table III), and to store the part of these, 
data required for future use in appropriate locations. It should be 
noted that only values of I c and T £ along the nadir direction [i.e., 

6 = 0° ; or A('1,N,L) and B(1„N, L)] are stored; values of these quanti- 
ties in the remaining 17 directions are discarded. At the end of a 
successful execution of- this part of the program, values of 1^, T^, and 
3 for the .L-th . wavelength,' for the J-th atmospheric model, and for 
the N-th value -of- 0 O (not 'for S) can be found in the location 
EIC(L, J1,N), TIC(L, Jl, N) ,• and SBAR(L, Jl) , respectively, where. J1 = 

J - 1 ■ for 1,000 mb models, and Jl = J-2 for 400 mb models. 

The next input to the program (value of the parameter NMOD; the 
atmospheric model number for which the printout is desired) is in a con- 
versational mode (at a terminal) under VM/CMS. Execution of the state- 
ment No. SITC 60 leads to the appearance of the. following message at the 
terminal: 

MODEL NUMBER PLEASE? 

? 

A value for NMOD can then be typed and entered. This value is then read 
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in as the value of the quantity NMOD during execution of statement No. 
61. 


The program is terminated if the value of the parameter NMOD read 
in is less than, or equal to, zero (SITC 62). For values of NMOD great- 
er than zero, the quantities and R are evaluated and printed out 
for all ten values of 9 0 and a given R during the execution of the 
450 DO LOOP (SITC 71 through 227), and then this- process is repeated 
for the remaining values of R (500 DO LOOP; SITC 67 through 230). In- 
formation about the atmospheric model under study is printed out before 
entering the DO LOOP 500 (SITC 64 through 66). 


For' a given atmospheric model NMOD, and a given value of reflecti- 
vity [REF (I) ] as well as of the solar zenith angle [THETO(N) ] , values of 
1^(1, 6 q) for all six values of A are computed during the execution of 
the statements contained between the statements numbered SITC 75 and 77, 


by making use of Eq. (1). Computed value of 1^(0.3800, 0 q) is then 

Q 

used for computations of the quantity R^ ^ggg iqOO t Ec l s * (2) and (3)] 
using 1,000 mb (IGRP = 1) tables, [Note that 0.3800 ym radiation is un- 
affected by ozone absorption. Thus, values of EIC(6, IT,N), TIC(6; IT,N), 
and SBAR(6, IT) are the same for all values of the subscript IT in the 
range 1 through 10.], and Eqs. (2) and (3).’ Statements in the 158 DO 


LOOP (SITC 87 through 89) are then used to assign the value of 
C C 

R n oon/i i nnn to R^ , at the remaining five wavelengths (Step II of 

U « jOUU ) 1UUU AjlUUU 

§2.4). 
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The section -of the program contained between, the statements No. 

C - 

SITC 97 and 108 is then used for computing values of £2^ 1000 an< ^ 

C 

x 1000 ^ 0r x = and 2 [computations for the first pair (x = l) are- 

carried out only if 0 O 5 79.6°]. The quantity Q 


'x, 1000 


is set to 


zero, and the quantity is set to -100 for 'cases for which it 

c 

is not possible to evaluate ^qqo (Step III of §2.4). Computations 


of SI 


x, 1000 


and S 


x , 1000 


are performed by calling the subroutine OMESLO, 


and after using the values of R^ j^qqq (EFFALC) . 


If n 


xaooo 


= 0 for x = 1 and 2, i.e., none of the ozone values 


is available at stage SITC 109, it is not possible to proceed with com- 
putations of ; rJ j1000 ,- s £ j1000 3.000 using Steps IV and V of 


§2.4. Under such circumstances, we assign appropriate values to these 


quantities . and. transfer control to statement No.. SITC 169. If one or 

c 

both values of £ 2 ^ ^qqq are available, one of the pairs (i.e., param- 
eter ISLT) is selected for computations of Rq 3393 ^000 during execu- 
tion of the statements No. SITC .118, 119, and 120. 


The part of the program contained between statements No. SITC 127 
through 150 is for computations of the impTOVed effective albedo repre- 
sented by the symbol R^ ^qqq using the procedure outlined in Step V 
of §2.4. The parameter JT is assigned a value between 1 and 10 (11 
and 20) if the set of computations are to be performed with 1,000 (400) 

mb tables. Values of and S 1 innA are then obtained during 

x ,1000 x ,1000 

the execution of statements No. SITC 157 through 168. This- latter part 
of the program is identical to that contained between statements No. 
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SITC 97 through 108, except that values of EFFALI are used in the place 
of those for the quantity EFFALC. 

A second pass through the part of the program contained between 
SITC 84 and SITC 169 is then made for computations of various quantities 
listed in Step VII of §2.4. Since this second pass is for computations 
with the 400 mb tables, parameters IGRP and IT are set to 2 and to 11, 
respectively (SITC 170 through 173). 

Statements SITC 178 through 185 are used for satisfying the- con- 
straint described under SITC3 in §5.1. Values of ft and R are then 

e 

obtained following Step IX of §2.4 by executing statements No. SITC 186 
through 221. 

For computations for a given value of the Lambert reflectivity of' 
the underlying surface (variable I of the DO LOOP 500), a summation of 
R for all ten values of 0 q is accumulated at the location STALB(I). 
Mean values of R for each of the eight values of R are printed out 
after completing all work for the atmospheric model, NMOD. 
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A. 2 Listing of the FORTRAN Statements 
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noon 


SITC 

PROGRAM FOR ESTIMATING TOTAL OZONE FROM SIMULATED OBSERVATIONS SITC 

IN THE NADIR DIRECTION FOR THE NIMBUS - 4 CONFIGURATION# . • - SITC 

SITC 

REAL '* 4 TITB{14),AU8,10,6) f 8CI8, 10, 6i ,C(6j ,REF ( 8 ) , EIC ( 6, 20,10) , SITC 

1 TIC(6,20,10i ,SBAR(6,20) ,ALDA(6) ,THETOf 10) , EFFALC(6,2) , SITC 

2 EFFALIC6»2) ,OMEGIN C 1GI ,OMEGC< 2,2) ,OMEGIf2,2) , SLPC i 2,2) ,SLPII2,2) SITC 

3 » EIMC 6 ) SITC 

REAL * 4 STALBC 8) • SITC' 

DATA THETO /O #0 , 45# Of 60# 0, 70.0 , 75# 6 f 79. 6 , 82# 5 , 84#7» 86.7 ,90.0/ SITC 

DATA ALDA /O. 3125,0 . 31T5 ,0.3312,0 #3398 , 0 #3600, 0. 3800/ SITC 

DATA REF /O.O, 0# 1 , 0# 2, 0.3 , Q*4, 0.6 , 0 = 8, 1 , 0/ SITC 

DEFINE FILE 16 ( 200, 8720, L, IFF) SITC 

DATA OMEGIN /0= 200, G=25Q ,0=300 ,0=350 ,0 .400, 0=450, 0= 500, 0, 550, SITC 

1 0.600,0=650/ SITC 

10 FORMAT (1H1 ) SITC 

15 FORMAT CT3 , * MODEL NUMBER PLEASE ?®i SITC 

20 FORMAT {1H1/T30, * ATMOSPHERIC MODEL NUMBER : °,I5,T80, ' SITC 

1 ‘PROGRAM USED : SITC4V//T10, SITC 

2 'SURFACE PRESSURE (MB) $ « ,F10=2,T55 ? 'TOTAL OZONE AMOUNT ( ATM-CMJSITC 

3 : *,F10.3////T40,‘D U S T DATA ®///T9»'TYPE OF', SITC 

4 T24, * SIZE DISTRIBUTION' ,T48,' PART OF THE REFRACTIVE INDEX*, SITC 

5 T86, 'TOT AL'/T8, 'AEROSOLS' ,T30, 'FUNCTION® »T52, 'REAL* , SITC 

6 T62, 6 IMAGINARY*, T85, 6 AMOUNT'/ i SITC 

25 FORMAT {/T4, « STRATOSPHERIC 8 ,8X ,3A4,4X, 0P2F15 .3 , 1 PI E20. 2 ) SITC 

30 FORMAT I/T4, ' TROPOSPHERIC '» 9X » 3A4, 4X,0P2F15 . 3, XP1 E20»2/ ) SITC 

35 FORMAT ( /T2, ' THETA* ,T14, » US ING 1,000 MB SURFACE-PRESSURE TABLES', SITC 

1 T71, 'USING 400 MB SURFACE-PRESSURE TABLES' , T120, ' BEST' / SITC 

2 T2, 8 SUB* ,T14, 'COARSE VALUES* , T42 , 9 IMPROVED VALUES', SITC 

3 T72, 'COARSE VALUES® ,T97, * IMPROVED VALUES * ,T120, 'OZONE'/ SITC 

4 T2» 'ZERO ® ,T10, *ALB EDO 6 , T1 9 , * OMEG. A * ,T27, 0 GMEG = B® , T37 ,* ALBEDO* , SITC 

5 T46, ' OMEG® A® ,T54, ' OMEG. B * ,T66, * ALBEDO' ,T75, '0MEG=A',T83, SITC 

6 'OMEG.B® ,T93, * ALBEDO 0 ,T102» *0MEG.A e ,TI 10 , * OMEGo B* ,T120, 'ESTIo* , SITC 

7 T128 , ' DEV. '/ ) SITC 

40 FORMAT |/T10,'MQDEL NUMBER :*,I5, SITC 

1 T45, 'LAMBERT REFLECTIVITY OF THE UNDERLYING SURFACE = SITC 

2 F10.3/ ) - SITC 

45 FORMAT £F6= 2, 2X,3F8 =3 , 3X, 3F8=3 ,5X,3F8=3 ,3X? 3F8 = 3 , F9 =3,F7=3) SITC 

IW =8 - SITC 

JMIN = 2 SITC 

UMAX = II SITC 

110 CONTINUE SITC 

DO 130 J = JMIN,JMAX SITC 

READ (16'J) TI TB , A, B ,C SITC 

J1 = J - 1 SITC 

IF < J *GT . . 11 ) Jl = J - 2 SITC 

DO 120 L = 1,6 SITC 

DO 115 N = 1,10 SITC 

EIC(L»J1,N) = A (1 ,N ,Li SITC 

TIC(L,J1,N) = B(1,N,L) SITC 

115 CONTINUE SITC 

SBAR<L,J1) = C ( L) SITC 

120 CONTINUE SITC 

130 CONTINUE SITC 

IF ( JMAX .GT. 11 ) GO TO 135 SITC 


1 

2 

3 

4 

5 

6 
7 
3 
9 

10 

11 

12 

13 

14 

15 

16 
17 
13 

19 

20 
21 
22 

23 

24 

25 

26 
27 
23 

29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 
43 

49 

50 

51 

52 

53 

54 

55 
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ooonnoo anon non 


135 

140 


150 


155 


158 


159 


160 


JMIN - 13 
UMAX = 22 
GO TO 110 
CONTINUE 
WRITE {6*15} 

READ {5**} NMOD 

IF ( NMOD oLEo 0) GO TO 1000 

READ ( 16.® NM0DJ TITB,A,B,C 

WRITE UW,20l NMOD? TITBCU ,TITB(2) 

WRITE ( I W , 25 ) TITB(9),TXT8C10J S TITBC 11 } ,TITB { 5 J , TITB{ 6) » TI T8( 3) 

WRITE { IW»30i TITB{ 12) ,TITB{ 131 ,TITB(X4) ,TITBC?) S TITB(8) ,TITB(4) 

DO 500 I = 1,8 

STALBU) = 0.0 

WRITE { IW ,401 NMOD* REF{ I i 

WRITE- ( I W ,35 ) 

DO 450 N = 1,10 

COMPUTE SIMULATED. MEASUREMENTS OF THE INTENSITY., o . . 

DO 150 L = 1*6 

EIM(L) * Ai 1 sN* L) + Bd,N,L) * CREFCIJ/CloO - C«L) * REF(I)J) 
CONTINUE 

COMPUTE VALUES OF THE COARSE EFFECTIVE ALBEDO USING 1,000 MB 
(IGRP=1) TABLES, OR 400 MB UGRP=2i TABLES. 

I GRP = 1 
IT = 1 
CONTINUE 

TEMA = E I M ( 6 -) - EIC(6,IT»N) 

EFFALC(6, IGRP) = TEMA/ (TIC (6 , IT, N ) + -SBARC6, IT l * TEMA) 

DO 158 J = 1,4 

EFFALC( J, IGRP) = EFFALC(6, IGRP ) 

CONTINUE 

ESTIMATE COARSE TOTAL-OZONE AMOUNTS FOR BOTH WAVELENGTH PAIRS, 
AND DECIDE. THE PAIR FOR DETERMINING IMPROVED EFECTIVE ALBEDO. 


00 NOT USE THE FIRST PAIR IF THETA SUB ZERO IS GREATER THAN 79.6 

IF { N . LE. 61 GO TO 159 
OMEGCd, IGRP) = 0.0 
SLPC( 1,1 GRP ) = -100.0 
GO TO 160 
CONTINUE. 

CALL OMESLO ( El C, TI C , S BAR , 1 , 3 , N , EFFALC T IGRP , El M, OMEGA , SLOPE 1 ' 

OMEGCd, IGRP) = OMEGA 
SLPCdjIGRP) = SLOPE 
CONTINUE 

CALL OMESLO ( El C , TI C, SBAR , 2 ,4, N, EFFALC , IGRP , El M , OMEGA , SLOPE) 

OMEGC (2, IGRP) = OMEGA 
SLPC( 2, IGRP) = SLOPE 

I F (OMEGCC 1 , IGRP ) .GT. 0.0 .OR. OMEGC ( 2 , IGRP ) .GT. 0.0) GO TO 165 
OMEGI ( 1 , IGRP ) = 0.0 


SITC 

56 

SITC 

57 

SITC 

58 

SITC 

59 

SITC 

60 

SITC 

61 

SITC 

62 

SITC 

63 

SITC 

64 

SITC 

65 

SITC 

66 

SITC 

67 

SITC 

68 

SITC 

69 

SITC 

70 

SITC 

71 

SITC 

72 

SITC 

73 

SITC 

74 

SITC 

75 

SITC 

76 

SITC 

77 

SITC 

78 

SITC 

79 

SITC 

80 

SITC 

81 

SITC 

82 

SITC 

83 

SITC 

84 

SITC 

85 

SITC 

86 

SITC 

87 

SITC 

88 

SITC 

89 

SITC 

90 

SITC 

91 

SITC 

92 

SITC 

93 

SITC 

94 

.SITC 

95 

SITC 

96 

SITC 

97 

SITC 

98 

SITC 

99 

SITC 

100 

SITC 

101 

SITC 

102 

SITC 

103 

SITC 

104 

SITC 

105 

SITC 

106 

SITC 

107 

SITC 

108 

SITC 

109 

SITC 

110 
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oooooo oooooo 


OMEGI (2? IGRP) * 0.0 
SLPIU, IGRP) = -lOOcO 
SLPI(2,IGRP) = -100.0 
DO 162 J = 1 ,4 

EFFAL'H J, IGRP) = EFFALC ( J , IGRP ) 

162 CONTINUE 
GO TO 190 
165 CONTINUE 
I 5LT = 2 

IF t SLPCC1, IGRP) .GT. SLPC(2,IGRP) ) ISLT = 1 

» 

COMPUTE VALUES OF THE IMPROVED EFFECTIVE ALBEDO, 


INTERPOLATE BETWEEN TWO OZONE VALUES. 

TEMA = OMEGC ( ISLT , IGRP } 

DO 170 J = 2,10 

IF ( TEMA .GT. 0MEGINCJ5 ) GO TO 170 
JT * ( IGRP— 1 > * 10 + J - 1 
GO TO 175 
170 CONTINUE 

' JT = { IGRP-1) *10+10 

175 CONTINUE 

TEMB = EIMC4) - EIC(4,JT,NI 

TEMB = TEMB/ITICC4, JT,N) + SBAR14,JT) * TEMB) 

JT = JT + 1 

IF ( IGRP .EG. 2 ) GO TO 176 
IF { JT o GT. 10 ) JT « 10 
GO TO 177 

176 IF ( JT .GT. 20 ) JT = 20 

177 CONTINUE 

TEMC = EIM(4) - EIC(4,JT,N1 

TEMC = TEMC/(TIC(4,JT,N) + SBAR(4,JT) * TEMC) 

JT = JT — 1 

IF t JT .GT. 10 ) JT * JT - 10 

EFFAL 1(4? IGRP) = TEMB + 20.0 * (TEMC - TEMB) * (TEMA 
DO 180 J = 1,3 

■ EFFALI(J,IGRP) = EFFALI (4, IGRP) 

180 CONTINUE 

ESTIMATE IMPROVED TOTAL OZONE AMOUNTS FOR BOTH WAVELENGTH PAIRS. 

DO NOT USE THE- FIRST PAIR IF THETA SUB ZERO IS. GREATER THAN 79.6 

IF ( N .LE. 6 ) GO TO 185 
0MEGM1* IGRP) - 0.0 
SLPI(ltlGRP) = -100.0 
GO TO 187 
185 CONTINUE 

CALL OMESLO ( El C ,TI C , S BAR , 1 ,3 , N , EFFALI , I GRP , El M, OMEGA, SLOPir ) 

OMEGI ( 1, IGRP ) = OMEGA 
SLPI ( 1 , IGRP) = SLOPE 
187 CONTINUE 


OMEGI NCJT) 


SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
) SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
.SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 
SITC 


111 

112 

113 

114 

115 

116 

117 

118 
'119 
120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 
156- 

157 

158 

159 

160 
161 
162 

163 

164 

165 
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non 


CALL OMESLO ( E I C ,TI C ,S BAR » 2 , 4, N, EFFALI , IGRP , E I M » OMEGA , SLORC) 

SITC 

166 

OMEGI ( 2 9 1 GRP ) = OMEGA 


SITC 

167 

SLP 1(2, IGRP 1 = SLOPE 


-SITC 

168 

190 CONTINUE 


SITC 

169 

IF ( IGRP -EG. 2 ) GO TO 200 


SITC 

170 

IGRP = 2 


SITC 

171 

IT = 11 


SITC 

172 

GO TO 155 


SITC 

173 



SITC 

174 

SELECT THE BEST ESTIMATE . . 


SITC 

175 



SITC 

176 

200 CONTINUE 


SITC 

177 

DO 210 K = 1,2 


SITC 

178 

DO 205 J = 1,2 


SITC 

179 

TEMA = AB S ( OMEGI { J , K ) - 0MEGC(J,K13 


SITC 

180 

IF ( TEMA »LT. 0.030 ) GO TO 205 


SITC 

181 

OMEGI ( J , K ) = 0.0 


SITC 

182 

SLPI (JsKi = -100.0 


SITC 

183 

205 CONTINUE 


SITC 

184 

210 CONTINUE 


SITC 

185 

IPAIRA = 0 


SITC 

186 

IPAIRB = 0 


SITC 

187 

IF ( OMEGI (1,1) oGT. 0,0 .AND. 0MEGI(1,2) .GT. 0.0) 

IPAIRA = I 

SITC 

188 

IF I OMEGI (2,13 oGT. 0.0 .AND. 0MEGI(2,2) .GT. 0.0) 

IPAIRB = 1 

SITC 

189 

IF UPAIRA .EQ. 1 .AND. IPAIRB .EQ. 1) GO TO 300 


SITC 

190 

IF '(IPAIRA -EQ- 1) GO TO 230 


SITC 

191 

IF (IPAIRB .EQ. 1 ) GO TO 270 


SITC 

192 

BEST = 0.0 


SITC 

193 

TEMA = 0.5 * (EFFAL 1(4,1) + EFFALI(4,2)) 


SITC 

194 

IF ( TEMA .GT. 0.2 ) GO TO 215 


SITC 

195 

BEST = AMAX1 < OMEGI < 1,1) tOMEGI (2, 1) ) 


SITC 

196 

GO TO 380 


SITC 

197 

215 IF ( TEMA .LT. 0.8 ) GO TO 380 


SITC 

198 

BEST = AMAX1 ( OMEGI { 1, 2) , OMEGI ( 2, 2 ) ) 


SITC 

199 

GO TO 380 


SITC 

200 

230 ISLT = 1 


SITC 

201 

GO TO 330 


SITC 

202 

•270 ISLT = 2 


SITC 

203 

GO TO 330 


SITC 

204 

300 ISLT = 2 


SITC 

205 

IF (SLPK 1,1) .GT. SLPI (2*1) .AND. SLPI{1,2) .GT. 

SLPI l 2,211 

SITC 

206 

1 ISLT = 1 


SITC 

207 

330 CONTINUE 


SITC 

208 

TEMA = 0.5 * (EFFAL 1(4,1) + EFFALI(4,2)> 


SITC 

209 

IF ( TEMA .GT. 0.2 J GO TO 340 


SITC 

210 

BEST = OMEGI( ISLT ,1 ) 


SITC 

211 

GO TO 380 


SITC 

212 

340 IF ( TEMA .GE. 0.8 ) GO TO 360 


SITC 

213 

BEST = (1.0/0. 6 ) * ( (0 . 8 — TEMA) * OMEGI ( I SLT ,1 ) + 


SITC 

214 

1 (TEMA-0.2) * OMEGI ( ISLT, 2} ) 


SITC 

215 

GO TO 380 


SITC 

216 

360 BEST = OMEGI ( ISLT ,2) 


SITC 

217 



SITC 

219 

PRINT OUT-. 


SITC 

219 



SITC 

220 



. 380 CONTINUE 

STALBCI) * STALBU* + TEMA 
DEV = BEST - TITB(2 ) . 

WRITE C I W ? 45 ) THETOCN) s EFFALCC 4, 1 K C OMEGC ( K r l 3 ,K=1,2) , 

1 EFFALXC4, 1) ,{OMEGHK, i),K=I,2) S EFFALC(4* 2) , C 0MEGC(K,2) ,K=I,2) * 

2 EFFALI(4j 2} , { OMEGI « K ? 2 ) , K=l, 2 J » BEST, DEV 
450 CONTINUE 

STALBUJ = 0.1 * STALBCI) 

IF ( I oEQe 2 .OR. I .EQ. 5 .OR. I .EQ. 8 ) WRITE (IW,10) 

500 CONTINUE 

WRITE UW?502) ( STALBC I ) , 1=1 , 8 ) 

502 FORMAT C/T4, 'EIGHT MEAN IMPROVED ALBEDOES s B ,8F8.3) 

GO TO 140 
1000 RETURN 
END • 


SITC 221 
SITC 222 
SITC 223 
SITC 224 
SITC 225 
SITC 226 
SITC 227 
SITC 228 
SITC 229 
SITC 230 
SITC 231 
SITC 232 
SITC 233 
SITC 234 
SITC 235 


122 



noooooooonoonoooooooooonoooooooo 


SUBROUTINE OMESLO { EIC ,TIC, SBAR, LAMl » LAM2, NTHETO, 

1 ALBD, IGRP, EIM, OMEGA, SLOPE) 

SUBROUTINE FOR ESTIMATING TOTAL OZONE AND THE SLOPE 

OF THE N - OMEGA CURVE IN THE REGION OF THE ESTIMATED OZONE--. 

EIC, TIC, AND SBAR AT 6 WAVELENGTHS, VIZ- 0 . 3 125, 0, 3175 , 
0.3312,0.3398,0.3600, AND 0.3800 MICROMETERS. 

E IC ,TIC, AND SBAR FOR 10 DIFFERENT OZONE AMOUNTS (0.200(0.050) 
0.650) FOR 1,000 MB SURFACE-PRESSURE MODELS, AND THEN FOR 
400 MB SURFACE-PRESSURE MODELS® . . . 

EIC AND TIC FOR 10 DIFFERENT VALUES OF THETA SUB ZERO... 

L AMI : SERIAL NUMBER OF THE SHORTER WAVELENGTH. 

LAM2 : SERIAL NUMBER OF THE LONGER WAVELENGTH® 

NTHETO s SERIAL NUMBER OF THETA SUB ZERO. 

ALBD s MEASURED ALBEDOES AT 6 WAVELENGTHS, AND FOR TWO GROUND 
PRESSURES... 

I GRP = 1 : USE TABLES CORRESPONDING TO 1,000 .MB® 

IGRP = 2 s USE TABLES CORRESPONDING TO 400 MB. 

EIM ; MEASURED INTENSITIES AT 6 WAVELENGTHS. 

OMEGA s ESTIMATED OZONE AMOUNT? SET TO 0.0 IF NO ESTIMATION ' 
IS POSSIBLE. 

SLOPE : SLOPE OF THE N - OMEGA CURVE IN THE REGION OF THE 

ESTIMATED TOTAL OZONE? SET TO -100.0 IF OMEGA - 0.0. 

20 FORMAT (/T5,«N0T POSSIBLE TO DETERMINE N -VALUE FOR 8 / 

1 T5, 'IGRP, NTHETO, LAM1, LAM2,TEMA,TEMB *• ,415, 1P2E12® 4/ 

2 T5, * EFFECTIVE ALBEDO AT LAM1 AND LAM2 -»,2F12.5$ 

REAL =5= 4 EIC(6,20,10),TICC6,20,10) , SBAR (6, 20) , ALBD (6, 2) , 

1 EIM(6),RAT(10),OMEGIN(10) 

DATA OMEGIN /O. 200, 0.250, 0.300,0.350,0.400, 0.450, 0® 500, 

1 0.550,0.600,0.650/ 

RATM = 100.0 * AL0G10 ( El M( LAM2 i/E IMCLAM1 ! } 

DO 200 J = 1,10 

JJ = (IGRP - 1) * 10 + J 

TEMA = EIC ( LAM 1 ,JJ» NTHETO) + TIC(LAM1, JJ , NTHETO I * 

1 ( ALBD( LAMl , IGRP )/( 1 .0 - SBARC LAM1 , J J ) *ALBD ( LAM1 , IGRP )) ) 

TEMB = EIC(LAM2,JJ, NTHETO) + TIC(LAM2,JJ, NTHETO) * 

1 ( ALBD(LAM2, IGRP)/( 1.0 - SB AR{ LAM2 , J J) *ALBD( LAM2 , IGRP ) ) ) 

IF (TEMA .LE. 1.0E-10 .OR. TEMB .LE= 1.0E-10) GO TO 150 
RAT ( J ) = 100.0 * AL0G10(TEMB/TEMA) 

GO TO 200 
150 CONTINUE 

WRITE (6, 20) IGRP, NTHETO, LAMl, LAM2, TEMA, TEMB, ALBD ( LAMl , IGRP ) , 

1 ALBD (LAM2, IGRP ) 

" RAT ( J ) = -9999.99 


OMESL 

1 

OMESL 

2 

OMESL 

3 

OMESL 

4 

OMESL 

5 

OMESL 

6 

OMESL 

7 

OMESL 

8 

OMESL 

9 

OMESL 

10 

OMESL 

11 

OMESL 

12 

OMESL 

13 

OMESL 

14 

OMESL 

15 

OMESL 

16 

OMESL 

17 

OMESL 

18 

OMESL 

19 

OMESL 

20 

OMESL 

21 

OMESL 

22 

OMESL 

23 

OMESL 

24 

OMESL 

25 

OMESL 

26 

OMESL 

27 

OMESL 

28 

OMESL 

29 

OMESL 

30 

OMESL 

31- 

OMESL 

32 

OMESL 

33 

OMESL 

34 

OMESL 

35 

OMESL 

36 

OMESL 

37 

OMESL 

38 

OMESL 

39 

OMESL 

40 

OMESL 

41 

OMESL 

42 

OMESL 

43 

OMESL 

44 

OMESL 

45 

OMESL 

46 

OMESL 

47 

OMESL 

48 

OMESL 

49 

OMESL 

50 

OMESL 

51 

OMESL 

52 

OMESL 

53 

OMESL 

54 

OMESL 

55 
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200 CONTINUE 

DO 300 J = l s 9 

IF { RATM «GTo RATCJ+l) ) GO TO 300 
SLOPE = {RATU+II - RAT(J>) * 20-0 
OMEGA = OMEGINI J } -5- (RATM - RAT ( J J ) /SLOPE 
GO TO 350 
300 CONTINUE 

OMEGA = 0*0 
SLOPE * -XOOoO 
350 RETURN 
END 


OMESL 56 
OMESL 57 
OMESL 58 
OMESL 59 
OMESL 60 
OMESL 61 
OMESL 62 
OMESL 63 
OMESL 64 
OMESL 65 
OMESL 66 
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A. 3 Output from a Test Run: 


125 



126 


ATMOSPHERIC MODEL NUMBER s 5 

SURFACE PRESSURE (MB) i 1000*00 TOTAL OZONE AMOUNT (ATM- 




•DUST 


DATA 


TYPE OF- 

SIZE DISTRIBUTION 

PART OF THE 

REFRACTIVE 

AEROSOLS 

FUNCTION 

REAL 

IMAGINARY 

STRATOSPHERIC 


loOOO 

GoO 

TROPOSPHERIC 


loOOO 

OqO 



MODEL NUMBER s 

5 

- 

LAMBERT 

REFLECTIVITY 

OF THE UNDERLY] 

THETA 

USING 

1,000 

MB SURFACE- 

-PRESSURE -TABLES 

• 

USING 400 

SUB 

COARSE VALUES 

IMPROVED VALUES 

COARSE V/ 

ZERO 

ALBEDO 

OMEGcA 

OH EG* B ' 

ALBEDO 

GMEG.A 

OMEG.B 

ALBEDO 

OMEG.? 

0.0 

0.0 

-0.350 

0.350 

-0.000 

0.350 

0.350 

0. 105 

0.34 

45.00 

0.0 

0.350 

0.350 

-0.000 

0.350 

0.350 

0.120 

0.35 

60.00 

0.0 

0.350 

0.350 

-0.000 

0.350 

0.350 

0.140 

0.3& 

70.00 

0.0 

0.350 

0.350 

-0.000 

0.350 

0.350 

0.167 

0.3 & 

75.60 

0.0 

0.350 

0.350 

-0.000 

0.350 

0.350 

0.187 

0. 37S 

79.60 

0.0 

0.350 

0.350 

-0.000 

0.350 

0.350 

0.201 

0.37; 

82.50. 

0.0 

0.0 

0.350 

-0.000 

0.0 

0.350 

0.205 

0.0 

84.70 

0.0 

0.0 

0.350 

-0.000 

0.0 

' 0.350 

0.193 

0.0 

86.70 

0.0 

0.0 

0.350 

-0.000 

0.0 

0.350 

0.180 

0.0 

90.00 

0.0 

0.0 

0.350 

-0.000 

0.0 

0.350 

0.147 

0.0 


MODEL NUMBER : 

- 5 


LAMBERT 

REFLECTIVITY 

OF THE UNDERLY[ 

THETA 

USING 

1,000 

MB SURFACE-PRESSURE TABLES 


USING 400 

SUB 

COARSE VALUES 

IMPROVED VALUES 

COARSE Vi 

ZERO 

ALBEDO 

OMEG*A 

OMEG.B 

ALBEDO 

QMEG.A 

OMEG.B 

ALBEDO 

GMEG.5 

. 0.0 

0.100 

0.350 

0.350 

0.100 

0.350 

0.350 

0.183 

0.34] 

45. 00 

0.100 

0.350 

0.350 

0.100 

0.350 

0.350 

0.194 

0.35: 

60.00 

0.100 

0.350 

0.350 

0.100 

0.350 

• 0.350 

0.211 

0.35? 

70.00 

0.100 

0.35C 

0.350 

0.100 

0.350 

0.350 

0.233 

0.36: 

75.60 

0.100 

0.350 

0.350 

0.100 

■ 0.350 

0.350 

0.250 

0.36? 

79.60 

0.100 

0.350 

0.350 

0.100 

0.350 

0.350 

0.263 

0.37: 

82.50 

0.100 

0.0 

0.350 

0.100 

0.0 

0.350 

0.267 

0.0 

84.70 

0.100 

0.0 

0.350 

0.100 

0.0 

0.350 

0.262 

0.0 

86.70 

0.100 

0.0 ■ 

' 0.350 

0.100 

0.0 

0.350 

0.247 

0.0 

90.00 

0.100 

0.0 

0.350 

0.100 

0.0 

0.350 

0.217 

0.0 


/$ 4 - 



5 


PROGRAM USED s SXTC4 


Ll OZONE AMOUNT (ATM-CMl s 0.350 


ft 

= REFRACTIVE INDEX TOTAL 

IMAGINARY AMOUNT 


OoO Q.Q 

0.0 OoO 


fcTIVITY OF THE UNDERLYING SURFACE = 0.0 



USING 400 MB 

SURFACE- 

PRESSURE 

TABLES 


BEST 



COARSE VALUES 

IMPROVED VALUES 

OZONE 


|oB 

ALBEDO 

OMEG.A OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ESTI. 

DEV® 

bo 

0. 105 

0.346 

0.305 

0.166 

0.324 

0.275 

0.350 

0.000 

po 

0.120 

0.353 

0. 319 

0.189 

0.332 

0.291 

0.350 

0.000 

|50 

0.140 

0.360 

0.333 

0.216 

0.339 

0.309 

0.350 

0.000 

50 

0.167 

0-366 

0.345 

0.244 

0.346 

0.326 

0.350 

0.000 

50 

0.187 

0.372 

0.352 

0.259 

0.350 

0.336 

0.350 

0.000 

50 

0.201 

0.379 

0.356 

0.260 

0.354 

0.342 

0.350 

0.000 

50 

0.205 

0.0 

0.356 

0.252 

0.0 

0.345 

0.350 

0.000 

50 

.0. 198 

0.0 

0.356 

0.239 

0.0 

0.344 

0.350 

0.000 

50 

0.180 

0.0 

0.360 

0.226 

0.0 

0.340 

0.350 

0.000 

50 

0.147 

0.0 

0.450 

0.241 

0.0 

0.0 

0.350 

0.000 

CTIVITY 

i 

OF THE UNDERLYING 

SURFACE 

= 0. 

100 

- 

- 


. ' 

USING 400 MB 

SURFACErPRESSURE 

TABLES 


BEST 



COARSE VALUES 

IMPROVED VALUES 

OZONE 


.8 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG. B 

ESTI. 

DEV . 

50 

0.183 

0.348 

0. 318 

0.233 

0.336 

0.302 

0.350 

0.000 

50 

0.194 

0.353 

0.327 

0.252 

0.340 

0.311 

0.350 

0.000 

50 

0.211 

0.357 

0.337 

0.275 

0.345 

0.322 

0.350 

0.000 

50 

0.233 

0.363 

0.347 

0.299 

0.349 

0.333 

0.350 

0.000 

50 

0.250 

0.367 

0.352 

0.312 

0.352 

0.341 

0.350 

0.000 

50 

0.263 

0.373 

0.355 

0.313 

0.355 

0.345 

0.350 

-0.000 

■50 

0.267 

0.0 

0.356 

0.307 

0.0 

0.347 

0.350 

-0.000 

50 

0.262 

0.0 

0.356 

0.295 

0.0 

0.347 

0.350 

0.000 

50 

0.247 

0.0 

0.358 

0.284 

0.0 

0.344 

0.350 

0.000 

50 

0.217 

0.0 

0.405 

0.290 

0.0 

0.0 

0.350 

0.000 


CL 
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MODEL NUMBER i 

5 


LAMBERT 

REFLECTIVITY 

OF THE UNDERLYI! 

THETA 

USING 

1,000 

MB SURFACE- 

•PRESSURE 

TABLES 


USING 400 ( 

SUB 

COARSE VALUES 

IMPROVED VALUES 

COARSE VA 

ZERO 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

0.0 

0.200 

0.350 

0.350 

0.200 

0.350 

0.350 

0.264 

0.353 

45.00 

0.200 

0.350 

0.350 

0.200 

0.350 

0.350 

0.271 

,0.355 

60.00 

0.200 

0.350 

0.350 

0.200 

0.350 

0.350 

0.284 

0-358 

70.00 

0.200 

0.350 

0.350 

0.200 

0.350 

0.350 

0.302 

0.361 

75.60 

0.200 

0.350 

0.350 

0.200 

0.350 

0.350 

0.316 

0.364 

79.60 

0.200 

0.350 

0.350 

0.200 

0.350 

0.350 

0.327 

0.369 

82. 50 

0.200 

0.0 

0.350 

0.200 

0.0 

0.350 

0.332 

0.0 

84.70 

0.200 

0.0 

0.350 

0.200 

0.0 

0.350 

0.328 

0.0 

86.70 

0.200 

0.0 

0.350 

0.200 

0.0 

0.350 

0.317 

0.0 

90.00 

0.200 

0.0 

0.350 

0.200 

0.0 

0.350 

0.290 

0.0 


MODEL NUMBER : 

5 


LAMBERT 

REFLECTIVITY 

OF THE UNDERLYIt 

THETA 

USING 

1,000 

MB SURFACE- 

-PRESSURE 

TABLES 


USING 400 f 

SUB 

COARSE VALUES 

I MPROVEO VALUES 

COARSE VAi 

ZERO 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

0.0 

0.300 

0.350 

0.350 

0.300. 

0.350 

0.350 

0.347 

0.353 

45.00. 

0.300 

0.350 

0.350 

0.300 

0.350 

0.350 

0.352 

0.358 

60.00 

0.300 

0.350 

-* 0.350 

0.300 

0.350 

0.350 

0.361 

0.359 

70. 00 

0.300 

0.350 

0.350 

0.300 

0.350 

0.350 

0.374 

0.361’ 

75.60 

0.300 

0.350 

0.350 

0.300 

0.350 

0.350 

0.385 

0.363 

79.60 

0.300 

0.350 

0.350 

0.300 

0.350 

0.35.0 

0.394 

0.366 

82.50 

0.300 

0.0 

0.350 

0.300 

0.0 

0.350 

0.399 

0.0 

84.70 

0.300 

0.0 

0.350 

0.300 

0.0 

0.350 

0.398 

0.0 

86.70 

0.300 

0.0 

0. 350 

0.300 

0.0 

0.350 

0.389 

0.0 

90. 00 

0.300 

0.0 

0.350 

0.300 

0.0 

0.350 

0.365 

0.0 1 


MODEL NUMBER : 

5 


LAMBERT 

REFLECTIVITY 

OF THE UNDERLYIN 

THETA 

USING 

1,000 

MB SURFACE- 

-PRESSURE 

TABLES 


USING 400 M 

SUB 

COARSE VALUES 

IMPROVED VALUES 

- ' COARSE VAL 

ZERO 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

0.0 • 

0.400 

0.350 

0.350 

0.400 

0.350 

0.350 

0.434 

0.363 

45.00 

0.400 

0.350 

0.350 

0.400 

0.350 

0.350 

0.435 

0.362 

60.00 

0.400 

0.350 

0.350 

0.400 

0.350 

0.350 

0.440 

0.361 

70.00 

0.400 

0.350 

0.350 

0.400 

0.350 

0.350 

0.448 

0.361 

75.60 

0.400 

0.350 

0.350- 

0.400 

0.350 

0.350 

0.457 

0.362 

79.60 

0.400 

0.350 

0.350 

0.400 

0.350 

0.350 

0.464 

0.364 

-82.50 

0.400 

0.0 

0. 350 

0.40Q 

0.0 

0.350 

0.469 

0.0 

84.70 

0.400 

0.0 

0.350 

0.400 

' 0.0 

0.350 

0.469 

0.0 

86.70 

0.400 

0.0 

0.350 

0.400 

0.0 

0.350 

0.464 

0.0 

90.00 

0.400 

0.0 

0.350 

0.400 

0.0 

0.350 

0.444 

0.0 





ItIVITY OF THE UNDERLYING SURFACE = 0.200 



USING 400 MB 

SURFACE-PRESSURE 

TABLES 


BEST 


!' 

COARSE VALUES 

IMPROVED VALUES 

OZONE 


|b 

ALBEDO 

OMEG.A I 

□MEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ESTI . 

DEV. 

1° 

0.264 

0.353 

0.332 

0.303 

0.346 

0.323 

0.350 

-0.000 

|o 

0® 271 

0.355 

G. 336 

0.318 

0.348 

0.327 

0.350 

- 0.000 

so 

0.284 

0.358 

0.342 

0.337 

0.350 

0.333 

0.350 

-0.000 

so 

0.302 

0.361 

0.349 

0.357 

0.352 

' 0.340 

0.350 

0.000 

io 

0.316 

0.364 

0.353 

0.368 

0.353 

0.345 

0.350 

0.000 

io 

0.327 

0.369 

0.355 

0.370 

0.356 

0.348 

0.350 

-0.000 

1 5 0 

0.332 

OeO 

0.356 

0.365 

0.0 

0.349 

0.350 

-0.000 

>0 

0.328 

0.0 

0.356 

0.355 

0.0 

0.349 

0.350 

”0.000 

>0 

0.317 

0.0 

0. 357 

0.345 

0.0 

0.347 

0.350 

-0.000 

50 

0.290 

0.0 

0.388 

0.348 

0.0 

OoQ 

0.0 

-0.350 

pIVITY 

OF THE UNDERLYING 

SURFACE 

= 0. 

.300 





USING 400 MB 

SURFACE- 

-PRESSURE 

TABLES 


BEST 



COARSE VALUES 

IMPROVED VALUES 

OZONE 


,B 

ALBEDO 

OMEG.A OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ESTI. 

DEV. 

>0 

0.347 

0.358 

0.344 

0.377 

. 0.354 

0.339 

0.351 

0.001 

50 

0.352 

0.358 

0.345 

0.388 

0.354 

0.340 

0.351 

0.001 

50 

0.361 

0.359 

0.348 

0.402 

0.354 

0.342 

0.351 

0.001 

50 

0.374 

' 0.361 

0. 352 

0.418 

0.355 

0.346 

0.351 

0.001 

50 

0.385 

0.363 

0.354 

0.428 

0.355 

0.349 

0.351 

0.001 

50 

0.394 

0.366 

0.356 

0.429 

0.357 

0.351 

0.350 

0.000 

50 

0.399 

0.0 

0.356 

0.426 

•0.0 

0.352 

0.350 

0.000 

SO 

0.398 

0.0 

0.356 

0.419 

0.0 

0.351 

0.350 

0.000 

50 

0.389 

0.0 

0.357 

0.410 

0.0 

0.350 

0.350 

-0.000 

50 

0.365 

0.0 

0.380 

0.411 

. 0.0 

0.0 

0.0 

-0.350 

:tivity 

OF THE 1 

JNDERLYING 

SURFACE 

= 0 , 

.400 






USING 400 MB 

SURFACE- 

•PRESSURE 

TABLES 


BEST 



COARSE VALUES 

IMPROVED VALUES 

OZONE 


.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ESTI. 

DEV. 

50 

0.434 

0.363 

0.354 

0.455 

0.361 

0.351 

0.354 

0.004 

50 

0.435 

0.362 

0.353 

0.461 

0.359 

0.350 

0.354 

0.004 

F° 

0.440 

0.361 

0.354 

0.471 

0.358 

0.350 

0.353 

0.003 

50 

0.448 

0.361 

0.355 

0.483 

0.357 

0.351 

0.353 

0.003 

50 

0.457 

0.362 

0.356 

0.491 

0.357 

0.352 

0.353 

0.003 

50 

0.464 

0.364 

0.356 

0.493 

0.358 

0.353 

0.351 

0.001 

50 

0.469 

0 = 0 

0.357 

0.491 

0.0 

0.353 

0.351 

0.001 

50 

0.469 

0.0 

0.357 

0.486 

0.0 

0.354 

0.351 

0.001 

50 

0-464 

0.0 

0.357 

0.479 

0.0 

0.353 

0.351 

0.001 

50 

0.444 

0.0 

0.375 

0.480 

0.0 

0.348 

0.349 

-0.001 
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MODEL NUMBER 


5 


LAMBERT REFLECTIVITY OF THE UNDERLYIhl 


THETA 

USING 

1,000 

MB SURFACE- 

•PRESSURE 

TABLES 


USING 400 V 

SUB 

COARSE VALUES 

IMPROVED VALUES 

COARSE VAL 

ZERO 

ALBEOO 

□MEG. A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A. 

0.0 

0.600 

0.350 

0.350 

0.600 

0.350 

Qj>350 

0.617 

0.371 

45.00 

0.600 

0.350 

0.350 

0.600 

0.350 

0.350 

0.610 

0.369 

60.00 

0.600 

0.350 

0.350 

0.600 

0.350 

0.350 

0.607 

0.366 

70.00 

0.600 

0.350 

0.350 

0.600 

0.350 

0.350 

0.607 

0.363 

75.60 

0.600 

0.350 

0.350 

0.600 

0.350 

0.350 

0.609 

0.362 

79.60 

0.600 

0.350 

0.350 

0.600 

0.350 

0.350 

0.612 

0.363 

82.50 

0.600 

0.0 

0.350 

0.600 

0.0 

0.350 

0.617 

0.0 

84.70 

0.600 

0.0 

0.350 

0.600 

0.0 

0.350 

0.622 

0.0 

86.70 

0.600 

0.0 

0.350 

0.600 

0.0 

0.350 

0.624 

0.0 

90.00 

0.600 

0.0 

0.350 

0.600 

0.0 

0.350 

0.610 

0.0 


MODEL NUMBER S- 

5 


LAMBERT 

REFLECTIVITY 

OF THE UNDERLYIN 


THETA 

USING 

1,000 

MB SURFACE- 

-PRESSURE 

TABLES 


. USING 400 M 

SUB 

COARSE 

; VALUES 

IMPROVED VALUES 

COARSE VAlJ 

ZERO 

ALBEDO OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 1 

0.0 

0. 800 

0.350 

0.350 

0.800 

0.350 

0.350 

0.313 

0.377 

45.00 

0-800 

0.350 

0.350 

0.800 

0.350 

0.350 

0.800 

0.374 

60.00 

0.800 

0.350 

0.350 

0.800 

0.350 

0.350 

0.788 

0.370 

70.00 

0.800 

0.350 

0.350 

0.800 

0.350 

0.350 

0.779 

0 • 366 

75.60 

0.800 

0.350 

0.350 

0.800 

0.350 

0.350 

0.774 

0.364 

79.60 

0.800 

0.350 

0.350 

0.800 

0.350 

0.350 

0.774 

0.363 

82.50 

0.800 

0.0 

0.350 

0.800 

0.0 

0.350 

0.779 

0.0 

84.70 

0.800 

0.0 

0.350 

0.800 

OoG 

0.350 

0.738 

0.0 

86.70 

0.800 

0.0 

0.350 

0.800 

0.0 

0.350 

0.797 

0.0 

90.00 

0.800 

0-0 

0.350 

0-800 

0.0 

0.350 

0.789 

0.0 


MODEL NUMBER : 

5 


LAMBERT 

REFLECTIVITY 

OF THE UNDERLYIN 

THETA 

USING 

1 ,000 

MB SURFACE- 

-PRESSURE 

TABLES 


USING 400 M 

SUB 

COARSE 

: VALUES 

IMPROVED VALUES 

-- COARSE VALt 

ZERO 

ALBEDO OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ALBEDO 

OMEG.A 


0.0 

1.000 

0.350 

0.350 

1.000 

0.350 

0.350 

1.026 

0.382 

45.00 

1.000 

0.350 

0.350 

1.000 

0.350 

0.350 

1 o 006 

0.379 

60.00 

1.000 

0.350 

0.350 

1.000 

0.350 

0.350 

0.985 

0.374 

70.00 

1.000 

0.350 

0.350 

1.000 

0.350 

‘ 0.350 

. 0.966 

0.369 

75.60 

1.000 

0.350 

0.350 

1.000 

0.350 

0.350 

0.954 

0.366 

79.60 

1.000 

0.350 

0.350 

1.000 

0.350 

0.350 

0.951 

0.365 

82.50 

1.000 

0.0 

0.350 

1.000 

0.0 

0.350 

0.956 

0.0 

84.70 

1.000 

0.0 

0.350 

1.000 

0.0 

0.350 

0.970 

0.0 

86.70 

1.000 

0.0 

0.350 

1.000 

0.0 

0.350 

0.986 

0.0 

90.00 

1.000 

1 

0.0 

0.350 

1.000 

o.o 

0.350 

0.985 

\ 

0.0 

EIGHT 

1 

< 

MEAN IMPROVED ALBEDOES : 

0.114 

0.193 

0.273 

< 

■s 

0.355 0.439 

0.61 





0 = 600 


ECTIVITY OF THE UNDERLYING SURFACE = 



USING 400 MB 

SURFACE- 

PRESSURE 

TABLES 


BEST 



COARSE VALUES 

IMPROVED VALUES 

OZONE 


G.B 

ALBEDO 

OMEG.A. 

0MEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ESTI. 

DEV. 

350 

0=617 

0 = 371 

0 = 368 

0 = 624 

0.37C 

0.368 

0.364 

0 = 014 

350 

0=610 

0 = 369 

0=366 

0 = 621 

0 = 368 

0.365 

0 = 362 

0 = 012 

350 

0 = 607 

0 = 366 

0=363 

0 = 622 

0=365 

0 = 362 

0.360 

0=010 

350 

0=607 

0 = 363 

0=360 

0.625 

0.362 

0 = 359 

0.358 

0=008 

350 

0 = 609 

0 = 362 

0 = 359 

0.629 

0=360 

0 = 357 

0.357 

0=007 

350 

0 = 612 

0 = 363 

0=358 

0.631 

0.360 

0.357 

0.355 

0=005 

350 

0 = 617 

0.0 

0=358 

0.632 

OoO 

0=357 

0.355 

0=005 

350 

0 =622 

0=0 

0=358 

0 = 632 

0 = 0 

0.357 

0 = 355 

0=005 

350 

0 = 624 

0= Q 

0.360 

0 = 629 

0.0 

0 = 358 

0.356 

0=006 

350 

0=610 

0 = 0 

0.372 

0=631 

0 = 0 

0.361 

0.357 

0=007 


ECTIVITY 

□F THE UNDERLYING 

SURFACE 

- 0 = 

800 





USING 400 MB 

SURFACE- 

-PRESSURE 

TABLES 


BEST 



COARSE VALUES 

IMPROVED VALUES 

OZONE 


G.B 

ALBEDO 

OMEG.A OMEG.B 

ALBEDO 

OMEG.A 

OMEG.B 

ESTI. 

DEV. 

350 

0.813 

0 = 377 

0.378 

0.813 

0.377 

0 = 378 

0 = 377 

0=027 

350 

0.800 

0.374 

0.374 

0.801 

0.374 

0=374 

0.374 

0.024 

350 

0.788 

0.370 

0.370 

0 = 792 

0.370 

0.370 

• 0.370 

0=020 

350 

0.779 

0.366 

0 = 365 

0.786 

0.366 

0.365 

0.365 

0=015 

350 

0.774 

0.364 

0= 362 

0.786 

0.363 

0.362 

0.363 

0 = 013 

350 

0.774 

0.363 

0.360 

0.788 

0.362 

0.359 

0.359 

0.009 

350 

0 = 779 

0.0 

0=360 

0=793 

0.0 

0 = 359 

0.359 

0 = 009 

350 

0 = 788 

0.0 

0.361 

0.797 

0.0 

0 = 360 

0.360 

0.010 

350 

0.797 

0.0 

0.362 

0.799 

0 . 0 

0.362 

0.362 

0.012 

1350 

0.789 

0.0 

0.374 

0.803 

0.0 

0.368 

0.368 

0.018 


Activity of the underlying SURFACE = 1 = 000 



USING 

I 400 MB 

SURFACE- 

■PRESSURE 

TABLES 


BEST 


1 

COARSE VALUES 

IMPROVED VALUES 

OZONE 


JG.B 

ALBEDO OMEG.A OMEG.B 

ALBEDO 

QMEGo A 

OMEG.B 

ESTI. 

DEV. 

350 

1.026 

0.382 

0.383 

1.025 

0.382 

0.383 

0.382 

0.032 

350 

1.006 

0.379 

0.380 

1.003 

0.379 

0.380 

0.379 

0 = 029 

35 0 

0.985 

0.374 

0.375 

0.984 

0.374 

0.375 

0.374 

0.024 

350 

0.966 

0.369 

0.369 

0.969 

0.369 

0.369 

0.369 

0.019 

350 

0.954 

0.366 

0.365 

0.964 

0.365 

0.364 

0.365 

- 0.015 

350 

0.951 

0.365 

0.362 

0.968 

0.363 

0.361 

0.363 

0 = 013 

1350 

0.956 

0-0 

0.361 

0.975 

0.0 

0=360 

0.360 

0.010 

350 ' 

0.970 

0.0 

0.362 

0.984 

0.0 

0.361 

0.361 

0.011 

350 

0.986 

0.0 

0. 365 

0.992 

0.0 

0.364 

0.364 

0.014 

350 

0.985 

\ 

6.355 0.439 

0.0 

0.376 

0.998 

0.0 

0.372 

0.372 

0 = 022 

173 

0.614 

0.798 

0.993 






-a- 



